Dear pw_forum,<br><div id=""><p>    I am trying to do a series of constrained ab initio molecular dynamics to obtain the the free energy profie from one intermediate to the next. In order to obtain the free energy of constrained molecular dynamics, I have to obtain the values of Lagrange multipliers for the constraint technique. Can someone tell me where the values of Lagrange multipilers are printed when the constrained ab initio MD is running? Thank you very much. Below is a sample of my input file, I try to constrain the three O-H distance of a hydronium ion in liquid water. There should be a serial of Lagrande multipliers related to the constraint, right? How should I obtain the values of  Lagrande multipliers when doing the calculation?</p><p><br></p><p>-------------------------------------------------------------------------------</p><p>-------------------------------------------------------------------------------</p><p>&CONTROL</p><p>  calculation='cp',</p><p>  dt=3.0D0,</p><p>  prefix='test',</p><p>  restart_mode='from_scratch',</p><p>  nstep=99000,</p><p>  iprint=1,</p><p>  isave=50,</p><p>  ndr=91,</p><p>  ndw=92,</p><p>  outdir='/home/temp',</p><p>  pseudo_dir="/home",</p><p>  tprnfor=.TRUE.,</p><p>  tstress=.FALSE.,</p><p>  disk_io='default',</p><p>  verbosity='default',</p><p>  max_seconds=86D+3</p><p>/</p><p>&SYSTEM</p><p>  ibrav=1,</p><p>  celldm(1)=24.802,</p><p>  nat=193,</p><p>  ntyp=2,</p><p>  nspin=1,</p><p>  ecutwfc=30.D0,</p><p>  ecutrho=240.D0,</p><p>  nr1b=28,</p><p>  nr2b=28,</p><p>  nr3b=28,</p><p>  tot_charge=0.0</p><p>  vdw_corr='DFT-D',</p><p>  london_s6=0.75</p><p>  london_rcut=200</p><p> ! degauss=0.002</p><p> ! occupations='smearing'</p><p> ! smearing='mp'</p><p>/</p><p>&ELECTRONS</p><p>  electron_dynamics='verlet',</p><p>  ! electron_temperature='nose',</p><p>  ! ekincw=0.035d0,</p><p>  ! fnosee=200.d0</p><p>  emass=400.d0,</p><p>  emass_cutoff=2.5d0,</p><p>  ortho_max=400,</p><p>/ </p><p>&IONS</p><p>  ion_dynamics='verlet',</p><p>    ion_radius(1)=1.0</p><p>    ion_radius(2)=1.0</p><p>  ion_temperature='nose',</p><p>  tempw=300.0,</p><p>  fnosep=40.0,</p><p>  nhptyp=2,</p><p>  nhpcl=4,</p><p>/</p><p>ATOMIC_SPECIES</p><p>  H  1.000  H.pbe-van_bm.UPF</p><p>  O  16.000 O.pbe-van_bm.UPF</p><p>ATOMIC_POSITIONS angstrom</p><p>193</p><p>generated  0.000000</p><p>H 10.0154  5.9257  7.6968</p><p>O  1.3610  5.7350  7.5670</p><p>H  0.7690  5.2550  8.2140</p><p>H  1.1890  5.3950  6.6430</p><p>O 11.0980  4.3090  3.2760</p><p>H 10.5150  3.9560  4.0060</p><p>H 10.5570  4.4260  2.4430</p><p>O  9.1760  6.1450  8.1060</p><p>H  8.4380  6.7670  7.8460</p><p>H  8.9630  5.7250  8.9880</p><p>O  4.0590 11.0180  5.7560</p><p>H  4.5850 10.9720  6.6050</p><p>H  3.8910 11.9750  5.5200</p><p>O  4.2620  6.3560  7.1420</p><p>H  3.3990  6.3890  7.6460</p><p>H  4.0740  6.3280  6.1600</p><p>O  8.0830  1.2070  3.5110</p><p>H  7.6000  1.5430  2.7030</p><p>H  8.8620  0.6480  3.2260</p><p>O  4.4710  3.6790  4.2970</p><p>H  3.8940  3.4070  3.5260</p><p>H  4.4660  2.9570  4.9880</p><p>O 12.1370  1.4190 12.1310</p><p>H 12.4300  0.4900 11.9000</p><p>H 11.2650  1.6140 11.6830</p><p>O  6.7530 10.4470 10.6850</p><p>H  7.6260 10.7580 11.0620</p><p>H  6.0020 10.8870 11.1770</p><p>O  4.5670  1.7530  6.3570</p><p>H  4.5730  1.7960  7.3560</p><p>H  5.4070  1.3140  6.0370</p><p>O  2.6960  0.9500 11.1090</p><p>H  2.1290  1.3620 10.3950</p><p>H  2.6150  1.4850 11.9500</p><p>O  9.0360 12.2130  7.4380</p><p>H  9.2020 11.9490  8.3880</p><p>H  9.2660 11.4490  6.8350</p><p>O  9.1710  9.9580  5.6590</p><p>H  8.2800  9.6800  6.0190</p><p>H  9.0400 10.5900  4.8960</p><p>O  0.3290  2.3180  2.3960</p><p>H -0.2660  3.0080  2.8090</p><p>H  0.0890  2.2050  1.4320</p><p>O  0.8050  8.4680  4.5200</p><p>H  1.0010  9.2390  5.1260</p><p>H -0.1820  8.4060  4.3670</p><p>O  4.7880 10.2540  1.4240</p><p>H  5.7450  9.9650  1.3980</p><p>H  4.3370  9.8430  2.2160</p><p>O  9.1290  3.2360  4.7200</p><p>H  8.3060  3.7730  4.9060</p><p>H  8.8810  2.4070  4.2200</p><p>O  0.2220  7.4780 11.5880</p><p>H  0.9800  7.8360 11.0410</p><p>H  0.4170  7.6170 12.5590</p><p>O  2.8060  4.5670 12.3120</p><p>H  1.8550  4.2850 12.1880</p><p>H  3.2320  4.7050 11.4180</p><p>O  7.0040  5.0180  4.3520</p><p>H  6.8100  5.9840  4.5200</p><p>H  6.1580  4.4900  4.4290</p><p>O  9.7160  8.0730 11.9000</p><p>H 10.6530  8.0320 11.5550</p><p>H  9.0760  7.9070 11.1490</p><p>O  5.0130  0.6990  3.2560</p><p>H  5.4990 -0.1600  3.4170</p><p>H  4.0890  0.5030  2.9270</p><p>O  7.9480  5.1020  1.3710</p><p>H  8.8180  4.8270  0.9620</p><p>H  7.9820  4.9550  2.3590</p><p>O  9.1970 11.6630 10.5210</p><p>H 10.1280 11.2980 10.5540</p><p>H  9.2130 12.6370 10.7480</p><p>O  0.5590 10.4940  1.2100</p><p>H  1.0660  9.6460  1.3640</p><p>H -0.4110 10.3410  1.3970</p><p>O  4.6840 11.9510 11.8510</p><p>H  4.5340 11.3000 12.5950</p><p>H  3.8060 12.3230 11.5500</p><p>O  3.5460  9.2240  3.8870</p><p>H  2.8180  8.6630  4.2810</p><p>H  3.7040 10.0220  4.4680</p><p>O  6.1300  2.9030  2.1740</p><p>H  5.5590  2.1400  2.4760</p><p>H  5.9940  3.6830  2.7850</p><p>O  8.4120  5.2620 10.5220</p><p>H  8.3510  6.2060 10.8460</p><p>H  7.6680  4.7250 10.9200</p><p>O  2.9560  3.0680  2.2970</p><p>H  2.9150  3.8150  1.6330</p><p>H  2.0300  2.7420  2.4890</p><p>O  1.2320 10.6220  6.2540</p><p>H  2.2230 10.7580  6.2290</p><p>H  1.0200  9.8230  6.8170</p><p>O  7.3500  7.7510 10.4780</p><p>H  7.2760  8.7420 10.3740</p><p>H  6.6780  7.4350 11.1470</p><p>O  0.4990  8.5380  7.7490</p><p>H -0.4910  8.4110  7.8020</p><p>H  0.9230  7.7080  7.3860</p><p>O  6.2320  7.9070  4.3870</p><p>H  5.3260  8.0760  3.9990</p><p>H  6.3820  8.5200  5.1630</p><p>O  0.8910  7.8290  1.8590</p><p>H  1.7870  7.3860  1.9040</p><p>H  0.6080  8.1010  2.7780</p><p>O 10.1370  9.4020  1.8420</p><p>H 10.1350  8.9860  0.9330</p><p>H 10.3340  8.7020  2.5290</p><p>O  3.5210  6.6410  1.6150</p><p>H  3.0830  5.9930  0.9910</p><p>H  4.3940  6.9380  1.2280</p><p>O 11.9420  4.4990  8.9590</p><p>H 11.6670  3.5910  8.6440</p><p>H 11.2730  5.1760  8.6530</p><p>O  1.3570  1.3190  8.5430</p><p>H  1.6370  1.8290  7.7300</p><p>H  0.4010  1.0420  8.4500</p><p>O  1.7160  9.4610 10.1300</p><p>H  2.5010  8.8780 10.3380</p><p>H  1.2920  9.1570  9.2760</p><p>O  6.6930  9.4320  6.4250</p><p>H  6.3310 10.1370  7.0350</p><p>H  6.7180  8.5560  6.9070</p><p>O  5.5820 11.1770  8.3150</p><p>H  6.1290 10.9440  9.1190</p><p>H  4.7040 11.5550  8.6080</p><p>O 10.5300  7.7940  4.2490</p><p>H 10.0580  8.2750  4.9880</p><p>H 10.2650  6.8290  4.2560</p><p>O  6.7750  7.1400  7.7570</p><p>H  6.7570  7.2610  8.7500</p><p>H  5.8860  6.8030  7.4460</p><p>O  3.0640 11.5930  8.8380</p><p>H  2.8250 11.1630  9.7090</p><p>H  2.4030 12.3130  8.6290</p><p>O 12.0460 11.1760 11.1600</p><p>H 12.6960 10.6090 10.6530</p><p>H 12.2320 11.1000 12.1400</p><p>O  7.0060  1.6040 12.1740</p><p>H  6.7520  2.1640 12.9620</p><p>H  6.1950  1.1360 11.8230</p><p>O  9.9150 11.8570  3.1040</p><p>H 10.0680 10.9590  2.6920</p><p>H 10.7500 12.1580  3.5660</p><p>O  1.8120  2.2450  6.0360</p><p>H  1.6070  3.0940  5.5480</p><p>H  2.7940  2.0610  5.9870</p><p>O  7.7370  2.7910  7.5270</p><p>H  8.6900  2.7870  7.2250</p><p>H  7.1550  2.4310  6.7980</p><p>O  0.1580  4.8520 11.5190</p><p>H  0.2110  5.8380 11.6760</p><p>H  0.0570  4.6750 10.5400</p><p>O  2.4540  0.0700  2.6630</p><p>H  1.9300 -0.5490  2.0790</p><p>H  1.9240  0.9030  2.8260</p><p>O  6.1400  3.4830 10.0360</p><p>H  6.7650  3.4160  9.2580</p><p>H  6.3060  2.7220 10.6640</p><p>O -0.1950 -0.0500  4.6980</p><p>H  0.3050 -0.5570  5.4000</p><p>H  0.4020  0.6510  4.3060</p><p>O  3.5980  7.5940 10.6820</p><p>H  3.4020  6.8100 10.0920</p><p>H  4.5490  7.5520 10.9870</p><p>O  7.2340  9.0690  1.8630</p><p>H  7.0040  8.8500  2.8110</p><p>H  8.2280  9.1000  1.7600</p><p>O 11.1450  1.6330  7.7110</p><p>H 11.6850  1.6740  6.8710</p><p>H 10.5250  0.8500  7.6710</p><p>O  3.7990  4.8780  9.8470</p><p>H  3.4560  4.7140  8.9220</p><p>H  4.6690  4.3990  9.9700</p><p>O  9.4720  2.0130 11.1900</p><p>H  9.4700  2.4190 10.2760</p><p>H  8.5340  1.9590 11.5330</p><p>O  6.8380  0.5160  5.7870</p><p>H  7.5440  0.4850  6.4940</p><p>H  7.2700  0.5990  4.8880</p><p>O  3.4930  6.2490  4.3110</p><p>H  3.4410  6.6120  3.3810</p><p>H  3.8670  5.3210  4.2860</p><p>O  1.0350  4.8600  5.2220</p><p>H  0.3720  4.6050  4.5180</p><p>H  1.5960  5.6190  4.8930</p><p>O  6.2370  6.9790  0.3590</p><p>H  6.7510  6.2580  0.8240</p><p>H  6.2690  7.8160  0.9040</p><p>O 10.3120  4.1090  0.5870</p><p>H 11.1220  4.3070  0.0340</p><p>H  9.7760  3.3890  0.1470</p><p><br></p><p><br></p><p>CONSTRAINTS</p><p> 3 1d-8</p><p> 'distance' 8 1</p><p> 'distance' 8 9</p><p> 'distance' 8 10</p><p><br></p><p>------------------------------------------------------------------------------------------</p><p>------------------------------------------------------------------------------------------  </p>
<p><br></p><p><br></p><p>Your sincerely,</p><p>Jinfan</p><p style="margin-right: 0px; margin-left: 0px; padding: 0px; line-height: 1.5; border: 0px; list-style: none; color: rgb(51, 51, 51); font-family: Verdana, 宋体, Helvetica, sans-serif; font-size: 14px;"><font face="Arial"><br></font></p><p style="margin-right: 0px; margin-left: 0px; padding: 0px; line-height: 1.5; border: 0px; list-style: none; color: rgb(51, 51, 51); font-family: Verdana, 宋体, Helvetica, sans-serif; font-size: 14px;"><font face="Arial">CHEN Jinfan</font></p><p style="margin-right: 0px; margin-left: 0px; padding: 0px; line-height: 1.5; border: 0px; list-style: none; color: rgb(51, 51, 51); font-family: Verdana, 宋体, Helvetica, sans-serif; font-size: 14px;"><font face="Arial">PhD candidate</font></p><p style="margin-right: 0px; margin-left: 0px; padding: 0px; line-height: 1.5; border: 0px; list-style: none; color: rgb(51, 51, 51); font-family: Verdana, 宋体, Helvetica, sans-serif; font-size: 14px;">School of Energy and Environment</p><p style="margin-right: 0px; margin-left: 0px; padding: 0px; line-height: 1.5; border: 0px; list-style: none; color: rgb(51, 51, 51); font-family: Verdana, 宋体, Helvetica, sans-serif; font-size: 14px;">City University of Hong Kong</p><p style="margin-right: 0px; margin-left: 0px; padding: 0px; line-height: 1.5; border: 0px; list-style: none; color: rgb(51, 51, 51); font-family: Verdana, 宋体, Helvetica, sans-serif; font-size: 14px;">Kowloon, Hong Kong SAR</p><p style="margin-right: 0px; margin-left: 0px; padding: 0px; line-height: 1.5; border: 0px; list-style: none; color: rgb(51, 51, 51); font-family: Verdana, 宋体, Helvetica, sans-serif; font-size: 14px;">Email: jinfachen2-c@my.cityu.edu.hk</p></div>