Calculating parity violation energy difference between enantiomers at MP2 level using DIRAC package

[Kaushik Talukdar]

 Dear Users:

         Could somebody suggest an appropriate input format (with all the keywords/notations) for calculating parity violation energy difference between enantiomers at MP2 level using DIRAC package.

Thanking you in advance.

with regards,

Kaushik

[Visscher, L.]

Dear Kaushimk,

An example can be found in test/MP2_properties/mp2_gradient.inp I paste this below for convenience (with removal of for you irrelevant properties that are also calculated). Kept the somewhat increased SCF precision as you typically want something of higher precision. The implementation  and precise formalism used is discussed in 

J. van Stralen, L. Visscher, C. Larsen, H. J. Aa.Jensen, First-order MP2 molecular properties in a relativistic framework, Chem Phys. 311 (2005) 81–95. doi:10.1016/j.chemphys.2004.10.018.

best regards,

Luuk

**DIRAC
.TITLE
Analytical first-order MP2 properties
.WAVE F
.PROPER
.4INDEX
**WAVE FUNCTIONS
.SCF
.RELCCSD
*SCF
.FCKCNV
1.0D-9
**PROPERTIES
.PVC
**MOLTRA
.ACTIVE
energy -15.0 100.0 1.0
**RELCCSD
.GRADIENT
*CCENER
.NOSD
.NOSDT
*CCFOPR
.MP2G
.RELAXED
*END OF

--
You received this message because you are subscribed to the Google Groups "dirac-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to dirac-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/dirac-users/8c586728-9cd1-429d-9bfc-a496586b7975%40googlegroups.com.

[Kaushik Talukdar]

Dear Professor Visscher,

Thank you for your reply. I will use the given input to calculate relevant property and let you know if I find any difficulty.

Also, I shall go through the reference suggested by you.

Thanks and best regards,

Kaushik

To unsubscribe from this group and stop receiving emails from it, send an email to dirac...@googlegroups.com.

[Leonardo belpassi]

13th International Conference on Relativistic Effects in Heavy-Element
Chemistry and Physics (REHE2020), Assisi, Italy.

We kindly invite you to register and submit an abstract (deadline is the next 7th of February 2020)

to present your research at the ''13th International Conference on Relativistic Effects in

Heavy-Element Chemistry and Physics (REHE 2020)'', Assisi, Italy from 22
June to 26 June 2020.

The REHE conference series runs since 26 years now with typically about
100-120 participants, who are mostly working on heavy element physics
and chemistry. Typically, the atmosphere is quite informal but
scientifically effective.

Please visit the conference website at https://www.rehe2020.it/ for
further information.

We would like to remind you that the maximum number of delegates is

limited for this event (first come/first served) and we have a great location !!!!

The organizers

--

You received this message because you are subscribed to the Google Groups "dirac-users" group.

To unsubscribe from this group and stop receiving emails from it, send an email to dirac-users...@googlegroups.com.

[juan...@gmail.com]

Dear Dirac developers,

Is there a way to calculate the parity violation energy differences adding a corresponding term to the hamiltonian, previous to the SCF?
Would it be correct to add the term H_PV (http://www.diracprogram.org/doc/release-19/manual/properties.html?highlight=pvc#pvc) to the hamiltonian? Is it possible?

Thanks


Juan José

[Trond SAUE]

Hi Juan José,

you mean you would like to include the PV term in the SCF-calculation ? Technically, it should be possible, but its contribution is exceedingly small compared to the electronic energy and would therefore drown in numerical noise.

All the best,

   Trond

--

Trond Saue
Laboratoire de Chimie et Physique Quantiques
UMR 5626 CNRS --- Université Toulouse III-Paul Sabatier
118 route de Narbonne, F-31062 Toulouse, France

Phone : +33/561556361 Fax: +33/561556065
Mail : trond...@irsamc.ups-tlse.fr
Web : https://dirac.ups-tlse.fr/saue
DIRAC : http://www.diracprogram.org/
ESQC : http://www.esqc.org/

Book: Principles and Practices of Molecular Properties: Theory, Modeling, and Simulations

[Juan Jose Aucar]

Yes that's what I mean.

I'm plotting some vector fields using the VISUAL module and I'd like to see if I can get some differences this way.

How do I add the PV operator into the hamiltonian?

--

You received this message because you are subscribed to the Google Groups "dirac-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to dirac-users...@googlegroups.com.

To view this discussion on the web visit https://groups.google.com/d/msgid/dirac-users/be19f3ad-19af-1f6f-4553-000fbfa0c281%40irsamc.ups-tlse.fr.

[Trond SAUE]

Hi again,

start by running an SCF on your system + calculation of PV as an expectation value. In this run recover the file AOPROPER

pam .. -- get AOPROPER

In your build directory there is a program called labread.x Run

labread.x AOPROPER

and send the output to the list.

[Juan Jose Aucar]

Hi Trond,

Thank you very much for the answers.

These are the labels I was looking for (for my test system):

        13     ********  22Jan22   SY 1FTTF  PVCSe 01
        15     ********  22Jan22   SY 1FTTF  PVCO  02
        17     ********  22Jan22   SY 1FTTF  PVCI  03
        19     ********  22Jan22   SY 1FTTF  PVCCl 04


Now I should be able to add this terms to the hamiltonian as:

**HAMILTONIAN
.OPERATOR
 'PVCSe 01'
.OPERATOR
 'PVCO  02'
.OPERATOR
 'PVCI  03'
.OPERATOR
 'PVCCl 04'

but I'm missing the common factor, right?

The Fermi coupling constant
G_F = 2.22255×10^-14E_h a₀³
and the weak nuclear charge
Q_A^w=Z^ACp_V+N^ACn_V=Z^A(1−4sin²θ_W)−N^A
where DIRAC uses sin²θ_w=0.2319 so that
Q_A^W=0.0724 Z^A-N^A
should be taken into account so that for each term the COMFACTOR should be
2.22255×10^-14 * 0.0724 Z^A-N^A* (1/2√2)

Am I right?

Thanks again

--

You received this message because you are subscribed to the Google Groups "dirac-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to dirac-users...@googlegroups.com.

To view this discussion on the web visit https://groups.google.com/d/msgid/dirac-users/a95525c0-f6bd-9c24-cd3d-e6e195646eed%40irsamc.ups-tlse.fr.

[Juan Jose Aucar]

This are the terms that have been added



so it seems I'm doing something wrong... if under **HAMILTONIAN I introduce the operator's names without '' then DIRAC doesn't find the operators
I attach the outputs

[Trond Saue]

Hi Juan,

I do not have my laptop here, so I can not see the factors, but you can check that from the previous output.

All the best,

Trond

Sent from my iPad

On 23 Jan 2022, at 03:08, Juan Jose Aucar <juan...@gmail.com> wrote:



To view this discussion on the web visit https://groups.google.com/d/msgid/dirac-users/CAAG1oD_xha-HhHNfUmNEFRZmxt9Vra4U2DEF93rQ14m%2BeywJ3g%40mail.gmail.com.

[Juan Jose Aucar]

Hi Trond, 

Thanks again.
In the output (attached) I get the Q factor. 
To build the operator I also need the normalized nuclear charge densities and I can't find it in the one-electron operators availables.

The straightforward way would be to just enter the operators labels but I'm still getting the error described before (when I use the labels that I got from using labread.x it doesn't find them).

Thanks,
Juan José

To view this discussion on the web visit https://groups.google.com/d/msgid/dirac-users/F52F7BCE-D15F-46EA-82F8-A6E3322D1154%40irsamc.ups-tlse.fr.