Supporting information for: Simulating Third-Order Nonlinear Optical. Properties Using Damped Cubic Response. Theory within Time-Dependent Density
|
|
- Brett French
- 5 years ago
- Views:
Transcription
1 Supporting information for: Simulating Third-Order Nonlinear Optical Properties Using Damped Cubic Response Theory within Time-Dependent Density Functional Theory Zhongwei Hu, Jochen Autschbach, and Lasse Jensen, Department of Chemistry, The Pennsylvania State University, 104 Chemistry Building, University Park, Pennsylvania , USA Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York , USA S1
2 1 Derivation of the relationship between ε αβ and ε αβ According to the damped second-order CPKS equation, the damped second-order Lagrangian multiplier matrices with arguments (+ω 1,+ω 2 ) and ( ω 1, ω 2 ) are given as ε αβ (+ω 1,+ω 2 ) = G αβ (+ω 1,+ω 2 ) +G α (+ω 1 )U β (+ω 2 )+G β (+ω 2 )U α +ω 1 ) U α (+ω 1 )ε β (+ω 2 ) U β (+ω 2 )ε α (+ω 1 ) (1) +ε 0 U αβ (+ω 1,+ω 2 ) U αβ (+ω 1,+ω 2 )ε 0 +(ω 1 +ω 2 +iγ)u αβ (+ω 1,+ω 2 ) and ε αβ ( ω 1, ω 2 ) = G αβ ( ω 1, ω 2 ) +G α ( ω 1 )U β ( ω 2 )+G β ( ω 2 )U α ( ω 1 ) U α ( ω 1 )ε β ( ω 2 ) U β ( ω 2 )ε α ( ω 1 ) (2) +ε 0 U αβ ( ω 1, ω 2 ) U αβ ( ω 1, ω 2 )ε 0 (ω 1 +ω 2 +iγ)u αβ ( ω 1, ω 2 ) respectively. Here we note that the expression for the negative argument has already transformed +iγ to iγ by construction. Therefore, taking the adjoint of Eq. (2) leads to the following equation that maintains the sign of iγ ε αβ ( ω 1, ω 2 ) = G αβ ( ω 1, ω 2 ) +U β ( ω 2 )G α ( ω 1 )+U α ( ω 1 )G β ( ω 2 ) ε β ( ω 2 )U α ( ω 1 ) ε α ( ω 1 )U β ( ω 2 ) (3) +U αβ ( ω 1, ω 2 )ε 0 ε 0 U αβ ( ω 1, ω 2 ) (ω 1 +ω 2 +iγ)u αβ ( ω 1, ω 2 ). S2
3 Using the facts ε 0 = ε 0, U a ( ω) = U a (+ω), G a ( ω) = G a (+ω), ε a ( ω) = ε a (+ω) and G ab ( ω 1, ω 2 ) = G ab (+ω 1,+ω 2 ), and substituting the following equations ε a (+ω) = G a (+ω)+ε 0 U a (+ω) U a (+ω)ε 0 +(ω +iγ)u a (+ω) G a (+ω) = ε a (+ω) ε 0 U a (+ω)+u a (+ω)ε 0 (ω +iγ)u a (+ω) (4) into Eq. (3), one can obtain ε αβ ( ω 1, ω 2 ) = G αβ (+ω 1,+ω 2 ) +G α (+ω 1 )U β (+ω 2 )+G β (+ω 2 )U α (+ω 1 ) U α (+ω 1 )ε β (+ω 2 ) U β (+ω 2 )ε α (+ω 1 ) [ ] U α (+ω 1 )U β (+ω 2 )+U β (+ω 2 )U α (+ω 1 ) ε 0 ] +ε [U 0 α (+ω 1 )U β (+ω 2 )+U β (+ω 2 )U α (+ω 1 ) [ ] +(ω 1 +ω 2 +2iΓ) U α (+ω 1 )U β (+ω 2 )+U β (+ω 2 )U α (+ω 1 ) (5) +U αβ ( ω 1, ω 2 )ε 0 ε 0 U αβ ( ω 1, ω 2 ) (ω 1 +ω 2 +iγ)u αβ ( ω 1, ω 2 ). This can be rewritten using W αβ (+ω 1,+ω 2 ) = U α (+ω 1 )U β (+ω 2 )+U β (+ω 2 )U α (+ω 1 ) as ε αβ ( ω 1, ω 2 ) = G αβ (+ω 1,+ω 2 ) +G α (+ω 1 )U β (+ω 2 )+G β (+ω 2 )U α (+ω 1 ) U α (+ω 1 )ε β (+ω 2 ) U β (+ω 2 )ε α (+ω 1 ) W αβ (+ω 1,+ω 2 )ε 0 +ε 0 W αβ (+ω 1,+ω 2 ) +(ω 1 +ω 2 +2iΓ)W αβ (+ω 1,+ω 2 ) [ ] + W αβ (+ω 1,+ω 2 ) U αβ (+ω 1,+ω 2 ) ε 0 ] ε [W 0 αβ (+ω 1,+ω 2 ) U αβ (+ω 1,+ω 2 ) [ ] (ω 1 +ω 2 +iγ) W αβ (+ω 1,+ω 2 ) U αβ (+ω 1,+ω 2 ), (6) S3
4 where the fact U αβ ( ω 1, ω 2 ) = W αβ (+ω 1,+ω 2 ) U αβ (+ω 1,+ω 2 ) is also adopted. Eq. (6) can be further simplified as ε αβ ( ω 1, ω 2 ) = G αβ (+ω 1,+ω 2 ) +G α (+ω 1 )U β (+ω 2 )+G β (+ω 2 )U α (+ω 1 ) U α (+ω 1 )ε b (+ω 2 ) U β (+ω 2 )ε α (+ω 1 ) +ε 0 U αβ (+ω 1,+ω 2 ) U αβ (+ω 1,+ω 2 )ε 0 (7) +(ω 1 +ω 2 +iγ)u αβ (+ω 1,+ω 2 ) +iγw αβ (+ω 1,+ω 2 ), which is equivalent to ε αβ ( ω 1, ω 2 ) = ε αβ (+ω 1,+ω 2 )+iγw αβ (+ω 1,+ω 2 ) (8) by substituting Eq. (1) into Eq. (7). S4
5 2 Derivation of the reduced IDRI second hyperpolarizability The full expression of IDRI can be obtained by substituting ω 1 = ω, ω 2 = ω, and ω 3 = ω into the general SOS equation for the second hyperpolarizability. To obtain the reduced form of it (γ TPA ) when ω is far from any one-photon resonances, we start by eliminating all pure one-photon terms in the full expression. This leads to a total of eight two-photon terms, where the dominant ones are given as 1 3 m 0 { n 0 p 0 0 µ α m m µ δ n n µ γ p p µ β 0 (ω m0 ω iγ)(ω n0 2ω iγ)(ω p0 ω iγ) 0 µ α m m µ + δ n n µ β p p µ γ 0 (ω m0 ω iγ)(ω n0 2ω iγ)(ω p0 ω iγ) 0 µ δ m m µ + α n n µ γ p p µ β 0 (ω m0 ω +iγ)(ω n0 2ω iγ)(ω p0 ω iγ) } 0 µ δ m m µ + α n n µ β p p µ γ 0, (ω m0 ω +iγ)(ω n0 2ω iγ)(ω p0 ω iγ) (9) equivalent to 1 3 n { S αδ 0n ( ω +iγ,ω iγ ) S βγ 0n ω n0 2ω iγ The expression above can be further rewritten as ( ω +iγ,ω +iγ ) }. (10) 1 3 n { S αδ 0n + ( ) ( ) ω iγ,ω iγ S βγ 0n ω +iγ,ω +iγ ω n0 2ω iγ [ S αδ 0n ( ω +iγ,ω iγ ) S αδ 0n ( ω iγ,ω iγ ) ] S βγ 0n ω n0 2ω iγ ( ω +iγ,ω +iγ ) }, (11) where the first part corresponds to the reduced IDRI second hyperpolarizability (γ TPA ), and the second part that represents a two-photon interference term, equivalent to 1 3 n { m [ 0 µ α m m µ δ n ω m0 (ω +iγ) 0 µα m m µ δ n ω m0 (ω iγ) S5 ] Sβγ 0n ( ω +iγ,ω +iγ ) ω n0 2ω iγ }, (12)
6 was also eliminated as being found contributing insignificantly to the two-photon poles when ω is far from any one-photon resonances. S6
7 3 Comparison between the SZ and STO-3G basis sets for the response properties of LiH at the TDDFT level of theory In Table S1, we represent the static α, β and γ values calculated at the TDDFT level of theory. All values in Table S1 are given in a.u. and only non-duplicate ones are shown based on Kleinman symmetry. Note that all static values were obtained without the lifetime, i.e., Γ = 0 a.u.. Table S1: Simulated static α, β and γ for LiH at the TDDFT level ADF Response Dalton Response Dalton SOS Property α xx (0;0) α zz (0;0) β zxx (0;0,0) β zzz (0;0,0) γ xxxx (0;0,0,0) γ xxyy (0;0,0,0) γ xxzz (0;0,0,0) γ zzzz (0;0,0,0) Within the TDDFT framework, all the properties calculated using response theory through ADF are in good agreement with those calculated using both response theory and the SOS approach through Dalton. This is an indication that, though identical minimal basis sets are not available, adopting SZ basis set for ADF while STO-3G basis set for Dalton is a reasonable approach for comparison purpose. Despite the same basis set (STO-3G) is ensured, discrepancies between those Dalton response and Dalton SOS results can still be seen for NLO properties at the TDDFT level of theory. S7
8 4 Comparison between response theory and the SOS approach for the response properties of LiH at the FCI level of theory In Table S2, we represent the static α, β and γ values calculated at the FCI level of theory. Including the damping factor, in Table S3 we also represent the resonant α and β values calculated at the FCI level of theory. All values in Tables S2 and S3 are given in a.u. and only non-duplicate ones are shown based on Kleinman symmetry. Note that all static values were obtained without the lifetime, i.e., Γ = 0 a.u., whereas the resonant α and β values were obtained using Γ = a.u at frequency equal to and a.u., respectively. The resonant γ values are not provided here because the FCI response theory approach to damped γ is currently not available in Dalton. Table S2: Simulated static α, β and γ for LiH at the FCI level Dalton Response Dalton SOS Property α xx (0;0) α zz (0;0) β zxx (0;0,0) β zzz (0;0,0) γ xxxx (0;0,0,0) γ xxyy (0;0,0,0) γ xxzz (0;0,0,0) γ zzzz (0;0,0,0) Table S3: Simulated resonant α and β for LiH at the FCI level Dalton Response Dalton SOS Property α xx ( ω;ω) α zz ( ω;ω) β xxz ( 2ω;ω,ω) β zxx ( 2ω;ω,ω) β zzz ( 2ω;ω,ω) S8
9 Within the FCI framework, all the properties calculated using response theory are nearly identical to those calculated using the SOS approach. This is expected for such an exact wavefunction method when the same basis (STO-3G) is ensured. S9
10 5 Comparison between response theory and the SOS approachfor thenlo processesoflih at thetddft level of theory EOPE: βeope( ω;ω,0) ( 10 4 a.u.) one photon pole βeope( ω;ω,0) ( 10 4 a.u.) one photon pole 6 6 (a) TDDFT Response (b) TDDFT SOS Figure S1: Simulated EOPE spectra for LiH using (a) response theory and (b) the SOS approach at the TDDFT level of theory. The vertical - - line indicates the one photon resonance due to the excited state. OR: S10
11 βor(0;ω, ω) ( 10 4 a.u.) one photon pole βor(0;ω, ω) ( 10 4 a.u.) one photon pole 6 6 (a) TDDFT Response (b) TDDFT SOS FigureS2: SimulatedORspectraforLiHusing(a)responsetheoryand(b)theSOSapproach at the TDDFT level of theory. The vertical - - line indicates the one photon resonance due to the excited state. EFIOR: γefior(0;ω, ω,0) ( 10 7 a.u.) one photon pole γefior(0;ω, ω,0) ( 10 7 a.u.) one photon pole (a) TDDFT Response (b) TDDFT SOS Figure S3: Simulated EFIOR spectra for LiH using (a) response theory and (b) the SOS approach at the TDDFT level of theory. The vertical - - line indicates the one photon resonance due to the excited state. EFISHG: S11
12 γefishg( 2ω;ω,ω,0) ( 10 6 a.u.) one photonpole two photonpole γefishg( 2ω;ω,ω,0) ( 10 6 a.u.) one photonpole two photonpole (a) TDDFT Response (b) TDDFT SOS Figure S4: Simulated EFISHG spectra for LiH using (a) response theory and (b) the SOS approach at the TDDFT level of theory. The vertical - - and -. lines indicate the oneand two-photon resonances due to the excited states, respectively. THG: For thepair ofbothβ OR andγ OKE spectra, we notethat the zero imaginary part is due to the orientational averaging rather than all zero tensor components. For the pair of γ THG spectra, the magnitude and shape of the second three-photon poles, as well as the first two-photon poles, slightly differ from each other. This is because the band extrema distance (the one within 1.5 to 2.0 ev) in Figure S5(b) is smaller than that in Figure S5(a), which causes a more pronounced doubly resonant effect that contributed by the overlap between two- and three-photon resonances. S12
13 γthg( 3ω;ω,ω,ω) ( 10 6 a.u.) one photon pole two photon pole three photon pole (a) TDDFT Response one photon pole two photon pole three photon pole (b) TDDFT SOS Figure S5: Simulated THG spectra for LiH using (a) response theory and (b) the SOS approach at the TDDFT level of theory. The vertical - -, -., and.. lines indicate the one-, two- and three-photon resonances due to the excited states, respectively. S13
14 6 Derivation of the spectral representation for the first hyperpolarizability We start with the first order transformation matrix (U) that can be obtained by solving the follow set of linear equations Uα (ω) U α ( ω) = A B B A ω V α V α, (13) wherev α isthedipolematrixinthemobasisandthespectral representationoftheresponse matrix reads ω B A X n n + ω n n ω Y N Y n 1 = 1 ω n +ω Y n X n Y n X n. (14) Then, the elementary form of the transition dipole moment can be written as µ α 0n,ia = V α ia (X n,ia +Y n,ia ) (15) and µ α n0,ia = V α ia(x n,ia +Y n,ia), (16) where i,j,... occupied (occ) orbitals, a,b,... virtual (virt) orbitals, and the Cartesian component is represented as α, β, γ... to avoid duplicate indices. This leads to the elementary form for the spectral representation of U as U α ia(+ω) = n X n,ia µ α n0,ia ω n ω + Y n,ia µα 0n,ia ω n +ω (17) S14
15 and U α ia ( ω) = n Y n,ia µ α n0,ia ω n ω + X n,ia µα 0n,ia ω n +ω. (18) Recall that the trace (Tr) operation in the 2n+1 expression for β can be rewritten as the sums over both occ and all (occ + virt) orbitals, S1 i.e., occ all Tr[nU α G β U γ ] = n UikG α β kl Uγ li (19) i kl and occ all Tr[nU α U γ G β ] = n Uik α Uγ kl Gβ li. (20) i kl Notice that Ref. S1 uses Tr[nU α U γ ε β ] in Eq. (20) but it is valid to replace ε with G here due to the zero diagonal block of U. Using the elementary form of U matrices and the first-order KS matrix G α ia (±ω) = Hα ia + j,b K HXC ia,jb Uα jb (±ω)+khxc ia,bj Uα bj ( ω), (21) Equations (19) and (20) can be rewritten as occ all n Uik( ω α σ )G β kl (+ω 1)U γ li (+ω occ all { 2) = n i [ m j,b kl X m,ik µ α m0,ik + Y ][ m,ik µα 0m,ik H β kl ω m +ω σ ω m ω + σ j,b i kl K HXC kl,jb ( Y Kkl,bj HXC n,jb µ β n0,jb + X n,jb µβ )][ 0n,jb ω n n ω 1 ω n +ω 1 p ( n X p,li µ γ p0,li X n,jb µ β n0,jb + Y n,jb µβ ) 0n,jb + ω n ω 1 ω n +ω 1 ω p ω 2 + Y p,li µγ 0p,li ω p +ω 2 ]} (22) S15
16 and occ all n Uik α ( ω σ)u γ kl (+ω 1)G β li (+ω occ all { 2) = n i [ m ( n kl X m,ik µ α m0,ik + Y ][ m,ik µα 0m,ik ω m +ω σ ω m ω σ X n,jb µ β n0,jb + Y n,jb µβ ) 0n,jb + ω n ω 2 ω n +ω 2 j,b p i kl X p,kl µ γ p0,kl ω p ω 1 + Y p,kl µγ 0p,kl ω p +ω 1 ][ H β li + j,b K HXC li,jb ( Y Kli,bj HXC n,jb µ β n0,jb + X n,jb µβ )]} 0n,jb, ω n n ω 2 ω n +ω 2 (23) respectively. We focus on the ALDA in response theory by extracting the terms containing thecouplingmatrixk HXC fromthe2n+1expressionofβ, inwhichtheone-photondominant ones are found as n p occ virt all exc. no. ij b kl mnp 0 [ Y m,ik µ α 0m,ik ω m ω σ K HXC kl,jb Xn,jbµ β n0,jb ω n ω 1 Xp,liµ γ p0,li ω p ω 2 Y m,ik µα 0m,ik ω m ω σ Xp,klµ γ p0,kl ω p ω 1 K HXC li,jb Xn,jbµ β n0,jb ω n ω 2 + Y m,ik µα 0m,ik ω m ω σ K HXC kl,bj Yn,jbµ β n0,jb ω n ω 1 Xp,liµ γ p0,li ω p ω 2 Y m,ik µα 0m,ik ω m ω σ Xp,klµ γ p0,kl ω p ω 1 K HXC li,bj Yn,jbµ β n0,jb ω n ω 2 ], (24) For the expression above, X n and Y n are the spectral representation of the response vectors obtained from the linear response equations, and p represents a summation over corresponding terms obtained by permuting ( 2ω,α), (+ω,β), and (+ω,γ). If one considers the S16
17 SHG case and takes the finite lifetime of the excited state into account, Eq. (24) becomes n p occ virt ij all b kl exc. no. mnp 0 [ Y m,ik µ α 0m,ik KHXC kl,jb X n,jb µ β n0,jb X p,liµ γ p0,li (ω m 2ω iγ)(ω n ω iγ)(ω p ω iγ) Y m,ik µα 0m,ik X p,klµ γ p0,kl KHXC li,jb X n,jb µ β n0,jb (ω m 2ω iγ)(ω p ω iγ)(ω n ω iγ) + Y m,ik µα 0m,ik KHXC kl,bj Y n,jb µ β n0,jb X p,liµ γ p0,li (ω m 2ω iγ)(ω n ω iγ)(ω p ω iγ) Y m,ik µα 0m,ik X p,klµ γ p0,kl KHXC li,bj Y n,jb µ β n0,jb (ω m 2ω iγ)(ω p ω iγ)(ω n ω iγ) ]. (25) Notice that, in addition to the part p Uα ( ω σ )[G β (+ω 1 ),U γ (+ω 2 )] being addressed above, the 2n+1 expression of β has another part that contains the second-order xc kernel, i.e., Tr[g ALDA xc D α ( ω σ )D β (+ω 1 )D γ (+ω 2 )]. (26) This term is also adiabatically approximated and potentially has a significant contribution to β in the vicinity of one-photon poles. However, as shown in Figure S6, we find it barely affects the shape or intensity of the one-photon structure in the SHG spectrum, and thus was not included within the derivation above. References (S1) Karna, S. P.; Dupuis, M. J. Comput. Chem. 1991, 12, S17
18 βshg( 2ω;ω,ω) ( 10 4 a.u.) No g xc No g xc Figure S6: Simulated SHG spectrum for LiH using response theory at the TDDFT level of theory. The dashed lines indicate the real and imaginary SHG values that were calculated by setting Eq. (26) equal to zero, respectively S18
Calculating frequency-dependent hyperpolarizabilities using time-dependent density functional theory
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 24 22 DECEMBER 1998 Calculating frequency-dependent hyperpolarizabilities using time-dependent density functional theory S. J. A. van Gisbergen Section Theoretical
More informationTDDFT in Chemistry and Biochemistry III
TDDFT in Chemistry and Biochemistry III Dmitrij Rappoport Department of Chemistry and Chemical Biology Harvard University TDDFT Winter School Benasque, January 2010 Dmitrij Rappoport (Harvard U.) TDDFT
More informationPolarizabilities. Patrick Norman
Polarizabilities Patrick Norman Division of Theoretical Chemistry and Biology School of Biotechnology KTH Royal Institute of Technology Stockholm, Sweden A perspective on nonresonant and resonant electronic
More informationSupporting Information
Supporting Information The Journal of Physical Chemistry A Determination of Binding Strengths of a Host-Guest Complex Using Resonance Raman Scattering Edward H. Witlicki, Stinne W. Hansen, Martin Christensen,
More informationDynamic second-order hyperpolarizabilities of Si 2 C and Si 3 C clusters. using coupled cluster singles-and-doubles response approach
Dynamic second-order hyperpolarizabilities of Si 2 C and Si 3 C clusters using coupled cluster singles-and-doubles response approach You-Zhao Lan* 1, Yun-Long Feng Zhejiang Key Laboratory for Reactive
More informationAccurate density functional calculations on frequency-dependent hyperpolarizabilities of small molecules
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 24 22 DECEMBER 1998 Accurate density functional calculations on frequency-dependent hyperpolarizabilities of small molecules S. J. A. van Gisbergen Section
More informationSUPPLEMENTARY INFORMATION
In the format provided by the authors and unedited. DOI: 10.1038/NPHYS3962 A global inversion-symmetry-broken phase inside the pseudogap region of YBa 2 Cu 3 O y Contents: S1. Simulated RA patterns for
More informationNonlinear optical response: New in silico Methods and Features
Nonlinear optical response: New in silico Methods and Features by Vladimir A. Goncharov Dissertation Submitted to the Faculty of the Graduate School of Vanderbilt University in partial fulfillment of the
More informationPure and zero-point vibrational corrections to molecular properties
Pure and zero-point vibrational corrections to molecular properties Kenneth Ruud UiT The Arctic University of Norway I UNIVERSITETET TROMSØ June 30 2015 Outline Why consider vibrational effects? General
More informationEJL R Sβ. sum. General objective: Reduce the complexity of the analysis by exploiting symmetry. Garth J. Simpson
e sum = EJL R Sβ General objective: Reduce the complexity of the analysis by exploiting symmetry. Specific Objectives: 1. The molecular symmetry matrix S. How to populate it.. Relationships between the
More informationCh. 2 NONLINEAR SUSCEPTIBILITIES
NONLINEAR OPTICS Ch. 2 NONLINEAR SUSCEPTIBILITIES Field notations Nonlinear susceptibility tensor : definition - 2 nd order NL susceptibility - 3 rd order NL susceptibility - n th order NL susceptibility
More informationCalculations of Linear and Nonlinear Optical Properties of H-Silsesquioxanes
J. Phys. Chem. B 2000, 104, 6737-6742 6737 Calculations of Linear and Nonlinear Optical Properties of H-Silsesquioxanes Wen-Dan Cheng,, Kai-Hua Xiang, Ravindra Pandey,*, and Udo C. Pernisz Department of
More information2.1 Experimental and theoretical studies
Chapter 2 NiO As stated before, the first-row transition-metal oxides are among the most interesting series of materials, exhibiting wide variations in physical properties related to electronic structure.
More informationHong Zhou, Guang-Xiang Liu,* Xiao-Feng Wang and Yan Wang * Supporting Information
Three cobalt(ii) coordination polymers based on V-shaped aromatic polycarboxylates and rigid bis(imidazole) ligand: Syntheses, crystal structures, physical properties and theoretical studies Hong Zhou,
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,5 8,.7 M Open access books available International authors and editors Downloads Our authors
More informationLinear response time-dependent density functional theory
Linear response time-dependent density functional theory Emmanuel Fromager Laboratoire de Chimie Quantique, Université de Strasbourg, France fromagere@unistra.fr Emmanuel Fromager (UdS) Cours RFCT, Strasbourg,
More informationThe frequency-dependent Sternheimer equation in TDDFT
The frequency-dependent Sternheimer equation in TDDFT A new look into an old equation Miguel A. L. Marques 1 Centre for Computational Physics, University of Coimbra, Portugal 2 LPMCN, Université Claude
More informationI. CSFs Are Used to Express the Full N-Electron Wavefunction
Chapter 11 One Must be Able to Evaluate the Matrix Elements Among Properly Symmetry Adapted N- Electron Configuration Functions for Any Operator, the Electronic Hamiltonian in Particular. The Slater-Condon
More informationSupporting Information. I. A refined two-state diabatic potential matrix
Signatures of a Conical Intersection in Adiabatic Dissociation on the Ground Electronic State Changjian Xie, Christopher L. Malbon, # David R. Yarkony, #,* Daiqian Xie,,%,* and Hua Guo,* Department of
More informationNonlinear Polarizabilities of Donor-Acceptor Substituted Conjugated Polyenes
11080 J. Phys. Chem. 1996, 100, 11080-11085 Nonlinear Polarizabilities of Donor-Acceptor Substituted Conjugated Polyenes Guanhua Chen and Shaul Mukamel* Department of Chemistry, UniVersity of Rochester,
More informationLecture 3. E α (ω) exp( iωt) dω = F 1 [E α ](t), (1 )
Nonlinear Optics 5A551 (200) Lecture notes Lecture Susceptibility tensors in the frequency domain The susceptibility tensors in the frequency domain arise when the electric field E α (t) of the light is
More informationFirst Order Hyperpolarizability and Homo-Lumo Analysis of L-Arginine Maleate (LArM) by Density Functional Theory Methods
51 Sciencia Acta Xaveriana An International Science Journal ISSN. 0976-1152 Volume 4 No. 2 pp. 51-58 September 2013 First Order Hyperpolarizability and Homo-Lumo Analysis of L-Arginine Maleate (LArM) by
More informationTime-Dependent Density-Functional Theory
Summer School on First Principles Calculations for Condensed Matter and Nanoscience August 21 September 3, 2005 Santa Barbara, California Time-Dependent Density-Functional Theory X. Gonze, Université Catholique
More informationAdvanced TDDFT II. II. Frequency-Dependent Kernels: Double Excitations. f xc
Advanced TDDFT II II. Frequency-Dependent Kernels: Double Excitations f xc Neepa T. Maitra Hunter College and the Graduate Center of the City University of New York Plan -- Double-Excitations in TDDFT
More information2.4 Properties of the nonlinear susceptibilities
2.4 Properties of the nonlinear susceptibilities 2.4.1 Physical fields are real 2.4.2 Permutation symmetry numbering 1 to n arbitrary use symmetric definition 1 2.4.3 Symmetry for lossless media two additional
More informationNonlinear optical phenomena on rough surfaces of metal thin films
PHYSICAL REVIEW B VOLUME 57, NUMBER 23 15 JUNE 1998-I Nonlinear optical phenomena on rough surfaces of metal thin films Evgeni Y. Poliakov, Vadim A. Markel,* and Vladimir M. Shalaev Department of Physics,
More information3: Many electrons. Orbital symmetries. l =2 1. m l
3: Many electrons Orbital symmetries Atomic orbitals are labelled according to the principal quantum number, n, and the orbital angular momentum quantum number, l. Electrons in a diatomic molecule experience
More informationInterference effects on the probe absorption in a driven three-level atomic system. by a coherent pumping field
Interference effects on the probe absorption in a driven three-level atomic system by a coherent pumping field V. Stancalie, O. Budriga, A. Mihailescu, V. Pais National Institute for Laser, Plasma and
More informationNWChem: Hartree-Fock, Density Functional Theory, Time-Dependent Density Functional Theory
NWChem: Hartree-Fock, Density Functional Theory, Time-Depent Density Functional Theory Hartree-Fock! Functionality! Input! Wavefunctions! Initial MO vectors! Direct and semidirect algorithms! Convergence,
More information3 Constitutive Relations: Macroscopic Properties of Matter
EECS 53 Lecture 3 c Kamal Sarabandi Fall 21 All rights reserved 3 Constitutive Relations: Macroscopic Properties of Matter As shown previously, out of the four Maxwell s equations only the Faraday s and
More informationThe interaction of light and matter
Outline The interaction of light and matter Denise Krol (Atom Optics) Photon physics 014 Lecture February 14, 014 1 / 3 Elementary processes Elementary processes 1 Elementary processes Einstein relations
More informationf xc Advanced TDDFT PDF Created with deskpdf PDF Writer - Trial :: II. Frequency-Dependent Kernels: Double Excitations
Advanced TDDFT II. Frequency-Dependent Kernels: Double Excitations f xc Neepa T. Maitra Hunter College and the Graduate Center of the City University of New York First, quick recall of how we get excitations
More informationTime-dependent density functional theory (TDDFT)
Advanced Workshop on High-Performance & High-Throughput Materials Simulations using Quantum ESPRESSO ICTP, Trieste, Italy, January 16 to 27, 2017 Time-dependent density functional theory (TDDFT) Ralph
More informationLinear-response excitations. Silvana Botti
from finite to extended systems 1 LSI, CNRS-CEA-École Polytechnique, Palaiseau, France 2 LPMCN, CNRS-Université Lyon 1, France 3 European Theoretical Spectroscopy Facility September 3, 2008 Benasque, TDDFT
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationCompound matrices and some classical inequalities
Compound matrices and some classical inequalities Tin-Yau Tam Mathematics & Statistics Auburn University Dec. 3, 04 We discuss some elegant proofs of several classical inequalities of matrices by using
More informationMagnetization-induced optical second- and third-harmonic generation in magnetic nanostructures
Laboratory of Nonlinear Optics of Nanostructures & Photonic Crystals Moscow State University Magnetization-induced optical second- and third-harmonic generation in magnetic nanostructures.shg.ru Т.V. Мurzina
More informationOptical and Photonic Glasses. Lecture 39. Non-Linear Optical Glasses III Metal Doped Nano-Glasses. Professor Rui Almeida
Optical and Photonic Glasses : Non-Linear Optical Glasses III Metal Doped Nano-Glasses Professor Rui Almeida International Materials Institute For New Functionality in Glass Lehigh University Metal-doped
More informationTime-dependent density functional theory (TDDFT)
04/05/16 Hands-on workshop and Humboldt-Kolleg: Density-Functional Theory and Beyond - Basic Principles and Modern Insights Isfahan University of Technology, Isfahan, Iran, May 2 to 13, 2016 Time-dependent
More informationNiS - An unusual self-doped, nearly compensated antiferromagnetic metal [Supplemental Material]
NiS - An unusual self-doped, nearly compensated antiferromagnetic metal [Supplemental Material] S. K. Panda, I. dasgupta, E. Şaşıoğlu, S. Blügel, and D. D. Sarma Partial DOS, Orbital projected band structure
More informationReview of Matrices. L A matrix is a rectangular array of numbers that combines with other such arrays according to specific rules.
Review of Matrices L A matrix is a rectangular array of numbers that combines with other such arrays according to specific rules. T The dimension of a matrix is given as rows x columns; i.e., m x n. Matrix
More informationPlasmonic eigenmodes in individual and bow-tie. graphene nanotriangles
Plasmonic eigenmodes in individual and bow-tie graphene nanotriangles Weihua Wang,, Thomas Christensen,, Antti-Pekka Jauho,, Kristian S. Thygesen,, Martijn Wubs,, and N. Asger Mortensen,, DTU Fotonik,
More informationComparison of frozen-density embedding and discrete reaction field solvent models for molecular propertiesw
PAPER www.rsc.org/pccp Physical Chemistry Chemical Physics Comparison of frozen-density embedding and discrete reaction field solvent models for molecular propertiesw Christoph R. Jacob,* a Johannes Neugebauer,
More informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Raman Centre for Atomic, Molecular and Optical
More informationSupporting information for Metamaterials with tailored nonlinear optical response
Supporting information for Metamaterials with tailored nonlinear optical response Hannu Husu, Roope Siikanen, Jouni Mäkitalo, Joonas Lehtolahti, Janne Laukkanen, Markku Kuittinen and Martti Kauranen Nonlinear
More informationChiroptical Spectroscopy
Chiroptical Spectroscopy Theory and Applications in Organic Chemistry Lecture 3: (Crash course in) Theory of optical activity Masters Level Class (181 041) Mondays, 8.15-9.45 am, NC 02/99 Wednesdays, 10.15-11.45
More information3.5 The random-phase approximation
154 3. LINEAR RESPONSE THEORY 3.5 The random-phase approximation Earlier in this chapter, we have demonstrated that many experimentally observable properties of solids can be expressed in terms of two-particle
More informationrequency generation spectroscopy Rahul N
requency generation spectroscopy Rahul N 2-11-2013 Sum frequency generation spectroscopy Sum frequency generation spectroscopy (SFG) is a technique used to analyze surfaces and interfaces. SFG was first
More informationQuantum Electrodynamical TDDFT: From basic theorems to approximate functionals
Quantum Electrodynamical TDDFT: From basic theorems to approximate functionals Ilya Tokatly NanoBio Spectroscopy Group - UPV/EHU San Sebastiàn - Spain IKERBASQUE, Basque Foundation for Science - Bilbao
More informationElectric properties of molecules
Electric properties of molecules For a molecule in a uniform electric fielde the Hamiltonian has the form: Ĥ(E) = Ĥ + E ˆµ x where we assume that the field is directed along the x axis and ˆµ x is the
More informationCorrelation spectroscopy
1 TWO-DIMENSIONAL SPECTROSCOPY Correlation spectroscopy What is two-dimensional spectroscopy? This is a method that will describe the underlying correlations between two spectral features. Our examination
More informationShort Course on Density Functional Theory and Applications VIII. Time-dependent DFT
Short Course on Density Functional Theory and Applications VIII. Time-dependent DFT Samuel B. Trickey Sept. 2008 Quantum Theory Project Dept. of Physics and Dept. of Chemistry trickey@qtp.ufl.edu Time-dependent
More informationDifferential of the Exponential Map
Differential of the Exponential Map Ethan Eade May 20, 207 Introduction This document computes ɛ0 log x + ɛ x ɛ where and log are the onential mapping and its inverse in a Lie group, and x and ɛ are elements
More informationThermal Emission in the Near Field from Polar Semiconductors and the Prospects for Energy Conversion
Thermal Emission in the Near Field from Polar Semiconductors and the Prospects for Energy Conversion R.J. Trew, K.W. Kim, V. Sokolov, and B.D Kong Electrical and Computer Engineering North Carolina State
More informationNonlinear Electrodynamics and Optics of Graphene
Nonlinear Electrodynamics and Optics of Graphene S. A. Mikhailov and N. A. Savostianova University of Augsburg, Institute of Physics, Universitätsstr. 1, 86159 Augsburg, Germany E-mail: sergey.mikhailov@physik.uni-augsburg.de
More informationIV. Calculations of X-ray Spectra in Real-space and Real-time. J. J. Rehr
TIMES Lecture Series SLAC-Stanford U March 2, 2017 IV. Calculations of X-ray Spectra in Real-space and Real-time J. J. Rehr Calculations of X-ray Spectra in Real-space and Real-time Goal: Real-space, real
More informationAdvanced TDDFT II. Double Excitations in TDDFT
Advanced TDDFT II. Double Excitations in TDDFT f xc Neepa T. Maitra Hunter College and the Graduate Center of the City University of New York First, quick recall of how we get excitations in TDDFT: Linear
More informationManipulator Dynamics (1) Read Chapter 6
Manipulator Dynamics (1) Read Capter 6 Wat is dynamics? Study te force (torque) required to cause te motion of robots just like engine power required to drive a automobile Most familiar formula: f = ma
More informationNIELINIOWA OPTYKA MOLEKULARNA
NIELINIOWA OPTYKA MOLEKULARNA chapter 1 by Stanisław Kielich translated by:tadeusz Bancewicz http://zon8.physd.amu.edu.pl/~tbancewi Poznan,luty 2008 ELEMENTS OF THE VECTOR AND TENSOR ANALYSIS Reference
More informationSupporting Materials
Supporting Materials Figure S1 Experimental Setup Page Figure S (a) (b) (c) Feynman Diagrams Page 3-6 Figure S3 D IR Spectra Page 7 Figure S4 Kinetic Model Page 8 Figure S5 Van t Hoff Plots Page 9 1 k
More informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
More informationA new generation of density-functional methods based on the adiabatic-connection fluctuation-dissipation theorem
A new generation of density-functional methods based on the adiabatic-connection fluctuation-dissipation theorem Andreas Görling, Patrick Bleiziffer, Daniel Schmidtel, and Andreas Heßelmann Lehrstuhl für
More informationarxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 26 Aug 2003
Circular dichroism simulated spectra of chiral gold nanoclusters: A dipole approximation arxiv:cond-mat/0308552v1 [cond-mat.mtrl-sci] 26 Aug 2003 Carlos E. Román-Velázquez, 1 Cecilia Noguez, 2, and Ignacio
More informationCalculations of X-ray Spectra in Real-space and Real-time
X-Ray Science in the 21st Century Calculations of X-ray Spectra in Real-space and Real-time J. J. Rehr, F. Vila, Y. Takimoto Department of Physics University of Washington Seattle, WA USA Time (s) KITP,
More informationChem 442 Review of Spectroscopy
Chem 44 Review of Spectroscopy General spectroscopy Wavelength (nm), frequency (s -1 ), wavenumber (cm -1 ) Frequency (s -1 ): n= c l Wavenumbers (cm -1 ): n =1 l Chart of photon energies and spectroscopies
More informationAssessment of range-separated time-dependent density-functional theory for calculating C 6 dispersion coefficients
1/10 Assessment of range-separated time-dependent density-functional theory for calculating C 6 dispersion coefficients Julien Toulouse 1,2, Elisa Rebolini 1, Tim Gould 3, John F. Dobson 3, Prasenjit Seal
More informationAdvanced Spectroscopy. Dr. P. Hunt Rm 167 (Chemistry) web-site:
Advanced Spectroscopy Dr. P. Hunt p.hunt@imperial.ac.uk Rm 167 (Chemistry) web-site: http://www.ch.ic.ac.uk/hunt Maths! Coordinate transformations rotations! example 18.1 p501 whole chapter on Matrices
More informationOptical and Photonic Glasses. Lecture 37. Non-Linear Optical Glasses I - Fundamentals. Professor Rui Almeida
Optical and Photonic Glasses : Non-Linear Optical Glasses I - Fundamentals Professor Rui Almeida International Materials Institute For New Functionality in Glass Lehigh University Non-linear optical glasses
More informationPredicting the Frequency Dispersion of Electronic Hyperpolarizabilities on the Basis of Absorption Data and Thomas-Kuhn Sum Rules
J. Phys. Chem. C 2010, 114, 2349 2359 2349 Predicting the Frequency Dispersion of Electronic Hyperpolarizabilities on the Basis of Absorption Data and Thomas-Kuhn Sum Rules Xiangqian Hu, Dequan Xiao, Shahar
More informationMolecular Vibrations C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology
Molecular Vibrations C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Why Estimate Molecular Vibrations? Simulation of vibrational spectrum (identification of molecules)
More informationSupporting information
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2014 Supporting information Resonant Raman Spectra of Molecules with Diradical Character:
More informationGarth J. Simpson. Department of Chemistry Purdue University. Garth J. Simpson
Objectives: 1. Discuss the benefits of coordinateindependent visualization of molecular tensors. 2. Describe sagittary representations for the resonant molecular tensor. 3. Introduce space-filling hyper-ellipsoid
More informationBasic concepts to start Mechanics of Materials
Basic concepts to start Mechanics of Materials Georges Cailletaud Centre des Matériaux Ecole des Mines de Paris/CNRS Notations Notations (maths) (1/2) A vector v (element of a vectorial space) can be seen
More informationPost Hartree-Fock: MP2 and RPA in CP2K
Post Hartree-Fock: MP2 and RPA in CP2K A tutorial Jan Wilhelm jan.wilhelm@chem.uzh.ch 4 September 2015 Further reading MP2 and RPA by Mauro Del Ben, Jürg Hutter, Joost VandeVondele Del Ben, M; Hutter,
More informationP A L A C E P IE R, S T. L E O N A R D S. R a n n o w, q u a r r y. W WALTER CR O TC H, Esq., Local Chairman. E. CO O PER EVANS, Esq.,.
? ( # [ ( 8? [ > 3 Q [ ««> » 9 Q { «33 Q> 8 \ \ 3 3 3> Q»«9 Q ««« 3 8 3 8 X \ [ 3 ( ( Z ( Z 3( 9 9 > < < > >? 8 98 ««3 ( 98 < # # Q 3 98? 98 > > 3 8 9 9 ««««> 3 «>
More informationSupporting information for: Second-order Nonlinear Optical Properties of. Stenhouse Photoswitches: Insights from. Density Functional Theory
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the wner Societies 218 Supporting information for: Second-order Nonlinear ptical Properties of Stenhouse
More informationarxiv:hep-ex/ v1 17 Sep 1999
Propagation of Errors for Matrix Inversion M. Lefebvre 1, R.K. Keeler 1, R. Sobie 1,2, J. White 1,3 arxiv:hep-ex/9909031v1 17 Sep 1999 Abstract A formula is given for the propagation of errors during matrix
More informationSCF calculation on HeH +
SCF calculation on HeH + Markus Meuwly Department of Chemistry, University of Basel, Basel, Switzerland Abstract This document describes the main steps involved in carrying out a SCF calculation on the
More informationWe will use C1 symmetry and assume an RHF reference throughout. Index classes to be used include,
1 Common Notation We will use C1 symmetry and assume an RHF reference throughout. Index classes to be used include, µ, ν, λ, σ: Primary basis functions, i.e., atomic orbitals (AOs) {φ µ ( r)}. The size
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
More informationThe frequency-dependent Sternheimer equation in TDDFT
The frequency-dependent Sternheimer equation in TDDFT A new look into an old equation Miguel A. L. Marques 1 Centre for Computational Physics, University of Coimbra, Portugal 2 European Theoretical Spectroscopy
More informationInternet Electronic Journal of Molecular Design
ISSN 1538 6414 Internet Electronic Journal of Molecular Design October 2004, Volume 3, Number 10, Pages 651 661 Editor: Ovidiu Ivanciuc Proceedings of the Internet Electronic Conference of Molecular Design,
More informationBSE and TDDFT at work
BSE and TDDFT at work Claudio Attaccalite http://abineel.grenoble.cnrs.fr/ CECAM Yambo School 2013 (Lausanne) Optical Absorption: Microscopic View Direct and indirect interactions between an e-h pair created
More informationAppendix to Lecture 2
PHYS 652: Astrophysics 1 Appendix to Lecture 2 An Alternative Lagrangian In class we used an alternative Lagrangian L = g γδ ẋ γ ẋ δ, instead of the traditional L = g γδ ẋ γ ẋ δ. Here is the justification
More informationA Brief Introduction to Tensors
A Brief Introduction to Tensors Jay R Walton Fall 2013 1 Preliminaries In general, a tensor is a multilinear transformation defined over an underlying finite dimensional vector space In this brief introduction,
More informationTime-dependent linear response theory: exact and approximate formulations. Emmanuel Fromager
June 215 ISTPC 215 summer school, Aussois, France Page 1 Time-dependent linear response theory: exact and approximate formulations Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
More informationarxiv: v1 [cond-mat.str-el] 10 Aug 2010
ROHF Theory Made Simple Takashi Tsuchimochi and Gustavo E. Scuseria Department of Chemistry, Rice University, Houston, TX 775-1892 Department of Physics and Astronomy, Rice University, Houston, TX 775-1892
More informationMultipole moments. Dipole moment. The second moment µ is more commonly called the dipole moment, of the charge. distribution and is a vector
Dipole moment Multipole moments The second moment µ is more commonly called the dipole moment, of the charge distribution and is a vector µ = µ x ˆx + µ y ŷ + µ z ẑ where the α component is given by µ
More informationSupplementary Figure 1 Interlayer exciton PL peak position and heterostructure twisting angle. a, Photoluminescence from the interlayer exciton for
Supplementary Figure 1 Interlayer exciton PL peak position and heterostructure twisting angle. a, Photoluminescence from the interlayer exciton for six WSe 2 -MoSe 2 heterostructures under cw laser excitation
More informationMOLECULAR SPECTROSCOPY
MOLECULAR SPECTROSCOPY First Edition Jeanne L. McHale University of Idaho PRENTICE HALL, Upper Saddle River, New Jersey 07458 CONTENTS PREFACE xiii 1 INTRODUCTION AND REVIEW 1 1.1 Historical Perspective
More informationAlkali metals show splitting of spectral lines in absence of magnetic field. s lines not split p, d lines split
Electron Spin Electron spin hypothesis Solution to H atom problem gave three quantum numbers, n,, m. These apply to all atoms. Experiments show not complete description. Something missing. Alkali metals
More informationThe arrow of time, black holes, and quantum mixing of large N Yang-Mills theories
The arrow of time, black holes, and quantum mixing of large N Yang-Mills theories Hong Liu Massachusetts Institute of Technology based on Guido Festuccia, HL, to appear The arrow of time and space-like
More informationLecture 7 From Dirac equation to Feynman diagramms. SS2011: Introduction to Nuclear and Particle Physics, Part 2 2
Lecture 7 From Dirac equation to Feynman diagramms SS2011: Introduction to Nuclear and Particle Physics, Part 2 2 1 Dirac equation* The Dirac equation - the wave-equation for free relativistic fermions
More informationJaroslav Hamrle. October 21, 2014
Magneto-optical Kerr effect (MOKE) Jaroslav Hamrle (jaroslav.hamrle@vsb.cz) October 21, 2014 Photon-photon spectroscopies (absorption) I: Type of investigations (polarized light x non-polarized light,
More informationGrading. Class attendance: (1 point/class) x 9 classes = 9 points maximum Homework: (10 points/hw) x 3 HW = 30 points maximum
Grading Class attendance: (1 point/class) x 9 classes = 9 points maximum Homework: (10 points/hw) x 3 HW = 30 points maximum Maximum total = 39 points Pass if total >= 20 points Fail if total
More information[ ( )] + ρ VIII. NONLINEAR OPTICS -- QUANTUM PICTURE: 45 THE INTERACTION OF RADIATION AND MATTER: QUANTUM THEORY PAGE 88
THE INTERACTION OF RADIATION AND MATTER: QUANTUM THEORY PAGE 88 VIII. NONLINEAR OPTICS -- QUANTUM PICTURE: 45 A QUANTUM MECHANICAL VIEW OF THE BASICS OF N ONLINEAR OPTICS 46 In what follows we draw on
More informationAdvanced Vitreous State The Physical Properties of Glass
Advanced Vitreous State The Physical Properties of Glass Active Optical Properties of Glass Lecture 21: Nonlinear Optics in Glass-Applications Denise Krol Department of Applied Science University of California,
More informationTime-dependent density functional theory : direct computation of excitation energies
CECAM Tutorial Electronic excitations and spectroscopies : Theory and Codes Lyon 007 Time-dependent density functional theory : direct computation of excitation energies X. Gonze Université Catholique
More informationAll-optical generation of surface plasmons in graphene
All-optical generation of surface plasmons in graphene T. J. Constant, 1, S. M. Hornett, 1 D. E. Chang, 2, and E. Hendry 1 1 Electromagnetic Materials Group, Department of Physics, College of Engineering,
More information