How to compute Two-photon absorption cross-sections with ADC(2)/ADC(3) methods?

Hi all,

I am currently using Qchem 6.0 and according to the manual it is possible to compute two-photon absorption cross-sections (TPAs) with ADC methods through the matrix inversion technique using the keyword ADC_PROP_TPA, however in the output I don’t find this information. If I include the keyword ADC_PROP_ES2ES I am able to recover TPAs but computed with the Sum over states expressions which is not the quantity I am interested in.

This is the input that I am using:

$molecule
0 1
O 0.000000 0.000000 -0.069903
H 0.000000 0.757532 0.518435
H 0.000000 -0.757532 0.518435
$end
$rem
METHOD adc(2)
BASIS aug-cc-pvdz
AUX_BASIS rimp2-aug-cc-pvdz
ADC_PROP_ES True
ADC_PROP_TPA True
EE_STATES [1,1,0,1]
$end
and this is the part of the output related to one of excited states:
Excited state 1 (singlet, B2) [converged]

Term symbol: 1 (1) B2 R^2 = 7.51604e-07

Total energy: -76.0072675789 a.u.
Excitation energy: 6.967051 eV

Osc. strength: 0.054586
Trans. dip. moment [a.u.]: [ -0.000000, -0.565506, -0.000000]
<i|r^2|0> [a.u.]: [ 0.000000, -0.000000, -0.000000]

Dip. moment [a.u.]: [ 0.000000, 0.000000, 0.690673]
Total dipole [Debye]: 1.755515
<r^2> [a.u.]: [ 15.147806, 10.478238, 13.978771]
Total <r^2> [a.u.]: 39.604816

V1^2 = 0.9292, V2^2 = 0.0708

Important amplitudes:
occ i occ j vir a vir b v
---------------------------------------------------------
1 (B2) A 4 (A1) A 0.6139
1 (B2) A 6 (A1) A 0.2067
---------------------------------------------------------

As you can see there is no trace of TPAs, do you know if I am doing anything wrong?

Best Regards,
Carmelo

I am not familiar with this methodology but I found the following sample input, whose output prints the two-photon abs. matrix without the sum-over-states message that appears in other jobs:

$comment
TPA cross-sections of H2O (direct)
$end

$molecule
0 1
O         0.000000    0.000000    0.146192
H        -0.748943    0.000000   -0.461081
H         0.748943    0.000000   -0.461081
$end

$rem
jobtype = sp
method = adc(2)
n_frozen_core = 0
basis = cc-pvdz
ee_singlets = [2,0,0,0]
adc_prop_tpa = true
$end

Output looks like:

  Excited state 1 (singlet, A1)                                    [converged]
  ----------------------------------------------------------------------------
  Term symbol:  2 (1) A1                                    R^2 =  1.96980e-08

  Total energy:                                            -75.8347960040 a.u.
  Excitation energy:                                              10.780969 eV

  Osc. strength:                                                      0.101837
  Trans. dip. moment [a.u.]:           [  -0.000000,    0.000000,    0.620933]
  <i|r^2|0> [a.u.]:                    [  -0.534828,   -0.104117,   -0.396852]
  Two-photon absorption cross-section [a.u.]:                       594.495795
  Two-photon absorption matrix [a.u.]:
                                       |   -4.896661   -0.000000   -0.000000 |
                                       |   -0.000000   -0.786115    0.000000 |
                                       |   -0.000000    0.000000   -6.795354 |

Hi jhelbert,

thank you very much for your useful answer! Indeed in this way it works, it seems that the problem is the employment of the resolution of identity which prevents the printing of these informations. Do you think it is a bug or maybe they haven’t implemented it yet?

Morevoer the same keywords don’t work for adc(3), which in principle should be possible. Did you find an example also for adc(3)?

Best
Carmelo

That is the only example that I found, in the samples directory that should ship with Q-Chem (sp_adc_tpa.in). Could be this is not yet implemented with RI or for ADC(3), please check the manual and if it’s not clear from there you could contact Q-Chem support.

Hi jhelbert, I am checking again this problem and I realized that even for adc2 there is something a bit weird for the two-photon absorption, in fact even in the output that you sent me the TPA is very large and not coherent with other values in the literature (see the paper of the implementation of this method on qchem Calculations of nonlinear response properties using the intermediate state representation and the algebraic-diagrammatic construction polarization propagator approach: Two-photon absorption spectra | The Journal of Chemical Physics | AIP Publishing). I am wondering if it is a conversion issue or a bug, do you know to whom I should talk about it?

Best,
Carmelo

Please email Q-Chem customer support.

My understanding form previously raising such queries is that TPA is only implemented (currently) for ADC(2), and (as mentioned above) not for Resolution of Identity calculations