Excited state second hyperpolarizability
-
- Posts: 5
- Joined: 02 Feb 2020, 17:18
- First name(s): Sidd
- Last name(s): Pandey
- Affiliation: Purdue Physics
- Country: United States
Excited state second hyperpolarizability
Hi,
I want to calculate the cubic response (second hyperpolarizability) of CO2, but on an excited state instead of the ground state. It seems that Dalton does its calculations with respect to a reference state, which is usually the ground state. Is there a way I can set an excited state as the reference state, to get the second hyperpolarizability on the excited state?
Thanks
I want to calculate the cubic response (second hyperpolarizability) of CO2, but on an excited state instead of the ground state. It seems that Dalton does its calculations with respect to a reference state, which is usually the ground state. Is there a way I can set an excited state as the reference state, to get the second hyperpolarizability on the excited state?
Thanks
- magnus
- Posts: 524
- Joined: 27 Jun 2013, 16:32
- First name(s): Jógvan Magnus
- Middle name(s): Haugaard
- Last name(s): Olsen
- Affiliation: Aarhus University
- Country: Denmark
Re: Excited state second hyperpolarizability
You can get excited-state polarizabilities from the double residue of the cubic response function but not hyperpolarizabilities.
-
- Posts: 270
- Joined: 27 Aug 2013, 16:42
- First name(s): Kenneth
- Last name(s): Ruud
- Affiliation: UiT The Arctic University of Norway
- Country: Norway
Re: Excited state second hyperpolarizability
CO2 is small enough that you may consider also doing an MCSCF calculation, selecting the excited state of interest as the reference state. You can then get the second hyperpolarizability by performing a MCSCF cubic response calculation.
- magnus
- Posts: 524
- Joined: 27 Jun 2013, 16:32
- First name(s): Jógvan Magnus
- Middle name(s): Haugaard
- Last name(s): Olsen
- Affiliation: Aarhus University
- Country: Denmark
Re: Excited state second hyperpolarizability
Yes, of course. Didn't consider MCSCF 

-
- Posts: 5
- Joined: 02 Feb 2020, 17:18
- First name(s): Sidd
- Last name(s): Pandey
- Affiliation: Purdue Physics
- Country: United States
Re: Excited state second hyperpolarizability
@kennethruud: How exactly do I "select the excited state of interest as the reference"?
-
- Posts: 270
- Joined: 27 Aug 2013, 16:42
- First name(s): Kenneth
- Last name(s): Ruud
- Affiliation: UiT The Arctic University of Norway
- Country: Norway
Re: Excited state second hyperpolarizability
In the
**WAVE FUNCTION
....
*CONFIGURATION
.SYMMETRY
2 <- Symmetry of the excited state
.STATE
2 <- Second state (i.e. first excited state) of symmetry 2.
See Dalton 2018 manual pages 294-299.
Kenneth
**WAVE FUNCTION
....
*CONFIGURATION
.SYMMETRY
2 <- Symmetry of the excited state
.STATE
2 <- Second state (i.e. first excited state) of symmetry 2.
See Dalton 2018 manual pages 294-299.
Kenneth
-
- Posts: 5
- Joined: 02 Feb 2020, 17:18
- First name(s): Sidd
- Last name(s): Pandey
- Affiliation: Purdue Physics
- Country: United States
Re: Excited state second hyperpolarizability
Thank you @kennethruud for the detailed explanation. However, when I try to specify the excited state to be calculated, I get an "illegal keyword" error. My .dal file is as follows:
The error I get is:
Keyword .STATE not recognized for *CONFIG
Code: Select all
**DALTON INPUT <- must start all input files
.RUN RESPONSE <- run integrals, wave function and response
**INTEGRALS
.DIPLEN. <- dipole length integrals for the chosen basis
**WAVE FUNCTIONS <- wave function input
.HF
.MP2 <- starting orbital for MCSCF
.MCSCF
*SCF INPUT <- HF input
.DOUBLY OCCUPIED
2 1 1 0 2 1 1 0
*CONFIGURATION INPUT
.SYMMETRY
1 <- symmetry of the excited state
.STATE
2 <- second (first excited) state of symmetry 1
.SPIN MULTIPLICITY
1
.INACTIVE
3 0 0 0 0 0 0 0 <- 1s orbitals on all atoms are inactive
.CAS SPACE
3 2 2 0 3 1 1 0 <- irrep. ordering: Ag B3u B2u B1g B1u B2g B3g Au
.ELECTRONS
16 <- no. of valence (2s & 2p) electrons in CO2
**RESPONSE
*CUBIC
.DIPLEN
**END OF DALTON INPUT
Keyword .STATE not recognized for *CONFIG
-
- Posts: 270
- Joined: 27 Aug 2013, 16:42
- First name(s): Kenneth
- Last name(s): Ruud
- Affiliation: UiT The Arctic University of Norway
- Country: Norway
Re: Excited state second hyperpolarizability
My error, .STATE should be under *OPTIMIZATION, i.e.
....
*OPTIMIZATION
.STATE
2 <- second (first excited) state of symmetry 1
*CONFIGURATION INPUT
.SYMMETRY
1 <- symmetry of the excited state
.SPIN MULTIPLICITY
1
......
....
*OPTIMIZATION
.STATE
2 <- second (first excited) state of symmetry 1
*CONFIGURATION INPUT
.SYMMETRY
1 <- symmetry of the excited state
.SPIN MULTIPLICITY
1
......
-
- Posts: 600
- Joined: 15 Oct 2013, 05:37
- First name(s): Peter
- Middle name(s): Robert
- Last name(s): Taylor
- Affiliation: Tianjin University
- Country: China
Re: Excited state second hyperpolarizability
I would be surprised if you need the STATE keyword, given that you are running with symmetry, because I am fairly certain that the first excited state of CO2 is a different symmetry to the ground state (I think perhaps a Piu state). This would then be the lowest state in symmetry 2 or 3. In that case, however, there is the possibility of Renner-Teller distortion of a degenerate electronic state of a linear molecule, and it is most likely that at the equilibrium geometry the molecule will be bent. Even if it is not (unlike the Jahn-Teller effect, where the symmetry is always lowered for a degenerate state of a nonlinear polyatomic molecule, the Renner-Teller effect may produce a linear equilibrium geometry) and the molecule remains linear in the excited state, you will have to deal with the issue of the electric field perturbation lifting the degeneracy of the electronic state, which will complicate analyzing the results of the response calculation (assuming you are interested in all components of the hyperpolarizability tensor).
Best regards
Pete
Best regards
Pete
Who is online
Users browsing this forum: No registered users and 8 guests