I. INTRODUCTION
In the course of a study aimed at obtaining analytic frequency-dependent hypermagnetizabilities and Cotton–Mouton constants of molecules using London atomic orbitals (LAO's),1 we came to discover the occurrence of an unfortunate mistake in the determination of the sign of the paramagnetic contribution to the anisotropy of the hypermagnetizability published in the studies of Refs. 2–4. In the latter, moreover, the sign of the hyperpolarizability contribution to Buckingham birefringence was wrong. Here we correct these mistakes and elaborate on the consequences for the discussions and comparisons made in the original papers.
II.
Table I replaces the original Table IV of Ref. 2. Since the paramagnetic and diamagnetic contributions are roughly of the same magnitudes and of opposite signs, the effect of the correction of the sign of is rather dramatic, with greatly reduced with respect to the erroneous estimates of Ref. 2. The effect on the Cotton-Mouton constant and on the anisotropy of the refractive index is marginal. With respect to the discussion of Cotton-Mouton effect (CME) in Sec. V D in Ref. 2, the need for significantly large basis sets, such as enhancement with double augmentation of the correlation consistent basis sets, is strengthened. Likewise, the originally reported agreement with the Hartree–Fock estimate in Ref. 5, where a lower quality basis was used, is no longer found. Indeed now the Hartree–Fock estimate for given by Augspurger and Dykstra in Ref. 5 is about five times larger than the value we compute. The effects induced by the change of sign of on the temperature dependence of the Cotton–Mouton constant are such that they do not affect the appearance of Fig. 2 in Ref. 2. In the temperature range of that figure, the contribution of is at most 0.4% of the total effect, instead of 2.2% as stated in Ref. 2. We add that the wide range of variations noted in the original paper for with the choice of the functional appears to be magnified.
. CME at . Anisotropy at . All data in a.u., with in . is defined as observed for an induction field of 1 T and a pressure of 1 atm.
Wf/Kernel . | Basis . | . | . | . | . | . | (cgs) . | . |
---|---|---|---|---|---|---|---|---|
HF-SCF | aug-cc-pVDZ | 483.9 | 160.5 | 2.83 | 1.70 | |||
HF-SCF | aug-cc-pVTZ | 519.1 | 159.5 | 3.01 | 1.81 | |||
HF-SCF | d-aug-cc-pVDZ | 539.2 | 150.1 | 2.97 | 1.79 | |||
LDA | 619.5 | 100.9 | 3.15 | 1.90 | ||||
LB94/LDA | 519.0 | 3.83 | 2.30 | |||||
BLYP | 655.3 | 154.5 | 3.11 | 1.87 | ||||
B3LYP/LDA | aug-cc-pVDZ | 491.1 | 173.5 | 2.59 | 1.56 | |||
B3LYP/BLYP | 462.6 | 156.8 | 2.63 | 1.58 | ||||
B3LYP | 584.1 | 148.9 | 3.06 | 1.85 | ||||
PBE | 620.4 | 101.1 | 3.15 | 1.90 | ||||
LDA | 622.9 | 106.5 | 3.18 | 1.91 | ||||
LB94/LDA | 516.0 | 3.6 | 3.57 | 2.15 | ||||
BLYP | 658.7 | 158.5 | 3.09 | 1.86 | ||||
B3LYP/LDA | aug-cc-pVTZ | 458.1 | 164.1 | 2.61 | 1.57 | |||
B3LYP/BLYP | 456.3 | 162.8 | 2.61 | 1.57 | ||||
B3LYP | 591.6 | 157.2 | 3.09 | 1.86 | ||||
PBE | 623.8 | 106.6 | 3.18 | 1.91 | ||||
CCSD | aug-cc-pVDZ | 527.2 | 165.1 | 2.76 | 1.66 | |||
d-aug-cc-pVDZ | 542.4 | 153.1 | 2.88 | 1.73 | ||||
a | ||||||||
Expt. | b | b | d | c d |
Wf/Kernel . | Basis . | . | . | . | . | . | (cgs) . | . |
---|---|---|---|---|---|---|---|---|
HF-SCF | aug-cc-pVDZ | 483.9 | 160.5 | 2.83 | 1.70 | |||
HF-SCF | aug-cc-pVTZ | 519.1 | 159.5 | 3.01 | 1.81 | |||
HF-SCF | d-aug-cc-pVDZ | 539.2 | 150.1 | 2.97 | 1.79 | |||
LDA | 619.5 | 100.9 | 3.15 | 1.90 | ||||
LB94/LDA | 519.0 | 3.83 | 2.30 | |||||
BLYP | 655.3 | 154.5 | 3.11 | 1.87 | ||||
B3LYP/LDA | aug-cc-pVDZ | 491.1 | 173.5 | 2.59 | 1.56 | |||
B3LYP/BLYP | 462.6 | 156.8 | 2.63 | 1.58 | ||||
B3LYP | 584.1 | 148.9 | 3.06 | 1.85 | ||||
PBE | 620.4 | 101.1 | 3.15 | 1.90 | ||||
LDA | 622.9 | 106.5 | 3.18 | 1.91 | ||||
LB94/LDA | 516.0 | 3.6 | 3.57 | 2.15 | ||||
BLYP | 658.7 | 158.5 | 3.09 | 1.86 | ||||
B3LYP/LDA | aug-cc-pVTZ | 458.1 | 164.1 | 2.61 | 1.57 | |||
B3LYP/BLYP | 456.3 | 162.8 | 2.61 | 1.57 | ||||
B3LYP | 591.6 | 157.2 | 3.09 | 1.86 | ||||
PBE | 623.8 | 106.6 | 3.18 | 1.91 | ||||
CCSD | aug-cc-pVDZ | 527.2 | 165.1 | 2.76 | 1.66 | |||
d-aug-cc-pVDZ | 542.4 | 153.1 | 2.88 | 1.73 | ||||
a | ||||||||
Expt. | b | b | d | c d |
III.
Tables II and III replace Tables IV and V, respectively, in Ref. 3. The whole paragraph at the end of the left column on p. 234314–7 in Ref. 3 should be replaced by the following:
Dynamic second electric-dipole hyperpolarizability and mixed electric-dipole hypersusceptibilities (see text for definitions). Atomic units, .
Wave function . | Basis . | . | . | . | . |
---|---|---|---|---|---|
HF-SCF | Sadlej | 3.24 | 4.57 | ||
LDA | 6.07 | 8.27 | |||
LB94/LDA | 4.19 | 5.59 | |||
BLYP | 6.61 | 8.99 | |||
B3LYP | 5.04 | 7.13 | |||
HF-SCF | aug-cc-pVDZ | 2.87 | 2.52 | ||
LDA | 5.39 | 6.04 | |||
LB94/LDA | 3.82 | 3.31 | |||
BLYP | 5.76 | 6.93 | |||
B3LYP | 4.64 | 5.05 | |||
HF-SCF | aug-cc-pVTZ | 3.25 | 4.71 | ||
B3LYP | 5.09 | 7.05 | |||
Expt. | a | b | c |
Wave function . | Basis . | . | . | . | . |
---|---|---|---|---|---|
HF-SCF | Sadlej | 3.24 | 4.57 | ||
LDA | 6.07 | 8.27 | |||
LB94/LDA | 4.19 | 5.59 | |||
BLYP | 6.61 | 8.99 | |||
B3LYP | 5.04 | 7.13 | |||
HF-SCF | aug-cc-pVDZ | 2.87 | 2.52 | ||
LDA | 5.39 | 6.04 | |||
LB94/LDA | 3.82 | 3.31 | |||
BLYP | 5.76 | 6.93 | |||
B3LYP | 4.64 | 5.05 | |||
HF-SCF | aug-cc-pVTZ | 3.25 | 4.71 | ||
B3LYP | 5.09 | 7.05 | |||
Expt. | a | b | c |
Linear birefringences of . . Atomic units except where noted.
Wave function . | Basis . | Kerr effecta . | CMEb . | Buckinghamc . | Jonesd . | ||||
---|---|---|---|---|---|---|---|---|---|
. | . | . | . | . | . | . | . | ||
HF-SCF | Sadlej | 1.91 | 7.66 | 0.44 | 2.46 | 11.31 | 11.52 | 3.54 | |
LDA | 2.98 | 11.98 | 0.23 | 1.31 | 10.44 | 20.85 | 6.41 | ||
LB94/LDA | 2.89 | 11.62 | 0.53 | 2.88 | 18.97 | 14.10 | 4.33 | ||
BLYP | 3.02 | 12.14 | 0.51 | 2.84 | 12.26 | 22.67 | 6.97 | ||
B3LYP | 2.69 | 10.81 | 0.61 | 3.42 | 12.05 | 17.99 | 5.53 | ||
HF-SCF | aug-cc-pVDZ | 1.95 | 7.85 | 0.23 | 1.29 | 11.57 | 6.36 | 1.96 | |
LDA | 3.01 | 12.08 | 1.73 | 9.73 | 10.34 | 15.23 | 4.68 | ||
LB94/LDA | 2.87 | 11.51 | 1.99 | 11.17 | 18.62 | 8.35 | 2.56 | ||
BLYP | 3.05 | 12.26 | 1.55 | 8.72 | 12.21 | 17.49 | 5.37 | ||
B3LYP | 2.72 | 10.95 | 1.34 | 7.54 | 12.06 | 12.73 | 3.91 | ||
HF-SCF | aug-cc-pVTZ | 1.93 | 7.76 | 1.91 | 10.72 | 11.68 | 11.88 | 3.65 | |
B3LYP | 2.69 | 10.79 | 2.09 | 11.72 | 11.99 | 17.78 | 5.46 | ||
Derived from Experiment | e | f | f | g | |||||
h | h |
Wave function . | Basis . | Kerr effecta . | CMEb . | Buckinghamc . | Jonesd . | ||||
---|---|---|---|---|---|---|---|---|---|
. | . | . | . | . | . | . | . | ||
HF-SCF | Sadlej | 1.91 | 7.66 | 0.44 | 2.46 | 11.31 | 11.52 | 3.54 | |
LDA | 2.98 | 11.98 | 0.23 | 1.31 | 10.44 | 20.85 | 6.41 | ||
LB94/LDA | 2.89 | 11.62 | 0.53 | 2.88 | 18.97 | 14.10 | 4.33 | ||
BLYP | 3.02 | 12.14 | 0.51 | 2.84 | 12.26 | 22.67 | 6.97 | ||
B3LYP | 2.69 | 10.81 | 0.61 | 3.42 | 12.05 | 17.99 | 5.53 | ||
HF-SCF | aug-cc-pVDZ | 1.95 | 7.85 | 0.23 | 1.29 | 11.57 | 6.36 | 1.96 | |
LDA | 3.01 | 12.08 | 1.73 | 9.73 | 10.34 | 15.23 | 4.68 | ||
LB94/LDA | 2.87 | 11.51 | 1.99 | 11.17 | 18.62 | 8.35 | 2.56 | ||
BLYP | 3.05 | 12.26 | 1.55 | 8.72 | 12.21 | 17.49 | 5.37 | ||
B3LYP | 2.72 | 10.95 | 1.34 | 7.54 | 12.06 | 12.73 | 3.91 | ||
HF-SCF | aug-cc-pVTZ | 1.93 | 7.76 | 1.91 | 10.72 | 11.68 | 11.88 | 3.65 | |
B3LYP | 2.69 | 10.79 | 2.09 | 11.72 | 11.99 | 17.78 | 5.46 | ||
Derived from Experiment | e | f | f | g | |||||
h | h |
given in SI units of . Pressure of 1 bar, and electric field strength of ;
given in cgs units; . defined for an induction field of 1 T, pressure of 1 atm, and ;
, , and ;
, , , and ;
from Ref. 16. We have extrapolated their fitted temperature-dependence linear equation to to calculate ;
from Ref. 13 with , . We have extrapolated their data to calculate to ;
from Ref. 14. We have extrapolated their fitted temperature-dependence linear equation to to calculate ;
measured at and across the temperature range of 304.1–453.5 K [Ref. 12]. We have extrapolated their data to 273.15 K to calculate .
“The value of is negative and very close to the center of the distribution of the experimental data (, as measured at ). The effect of electron correlation depends strongly on the functional, leading to a reduction of as much as for LB94/LDA, and to increases of up to 25% for the other functionals. With the aug-cc-pVTZ basis set, the increase is for B3LYP.”
No changes are needed for the rest of the Discussion in Ref. 3, and Fig. 3 thereof also remains unchanged.
IV.
Table IV replaces Table IV of Ref. 4. As a first consequence of the changes, in particular due to the reversal of the sign of , the revised value of the traceless quadrupole moment of originally measured in Ref. 6 is (instead of the previous incorrect value of ). The entry in the next to the last row in Table III of Ref. 4 must therefore be updated. Note then that, contrary to what we stated originally, the revision of the “apparent” quadrupole moment given in Ref. 6, made taking into account the nonvanishing contribution, brings the Buckingham-birefringence-derived experimental value closer to our ab initio theoretical best estimate . The inclusion of zero-point vibrational average further improves the comparison.
CME and BE for at and . Atomic units, with in cgs units of . is the birefringence defined for an induction field of 1 T and a pressure of 1 atm according to Ref. 17. for Buckingham birefringence is given for a pressure of and EFG of .
Wave function . | Basis . | CME . | Buckingham . | ||||
---|---|---|---|---|---|---|---|
. | . | . | . | . | . | ||
HF-SCF | aug-cc-pVDZ | 6.17 | 3.72 | 26.3 | 3.74 | ||
d-aug-cc-pVDZ | 3.46 | 2.09 | 33.6 | 3.81 | |||
aug-cc-pVTZ | 0.8 | 5.55 | 3.34 | 30.7 | 3.77 | ||
d-aug-cc-pVTZ | 3.9 | 5.35 | 3.22 | 34.0 | 3.64 | ||
aug-cc-pVQZ | 2.5 | 5.38 | 3.24 | 32.9 | 3.65 | ||
LDA | aug-cc-pVDZ | 13.70 | 8.25 | 45.7 | 4.31 | ||
d-aug-cc-pVDZ | 11.86 | 7.15 | 58.2 | 4.38 | |||
aug-cc-pVTZ | 5.3 | 12.82 | 7.72 | 53.3 | 4.47 | ||
d-aug-cc-pVTZ | 9.1 | 12.37 | 7.45 | 59.8 | 4.29 | ||
aug-cc-pVQZ | 7.0 | 12.46 | 7.50 | 57.1 | 4.32 | ||
B3LYP | aug-cc-pVDZ | 10.79 | 6.50 | 40.6 | 4.35 | ||
d-aug-cc-pVDZ | 8.78 | 5.29 | 51.8 | 4.39 | |||
aug-cc-pVTZ | 4.2 | 9.86 | 5.94 | 47.3 | 4.44 | ||
d-aug-cc-pVTZ | 7.6 | 9.44 | 5.69 | 53.2 | 4.26 | ||
aug-cc-pVQZ | 5.8 | 9.65 | 5.81 | 50.5 | 4.30 | ||
KT1 | aug-cc-pVDZ | 11.72 | 7.06 | 51.4 | 4.10 | ||
d-aug-cc-pVDZ | 9.15 | 5.51 | 64.5 | 4.14 | |||
aug-cc-pVTZ | 7.2 | 10.40 | 6.26 | 60.0 | 4.22 | ||
d-aug-cc-pVTZ | 11.0 | 10.13 | 6.10 | 66.9 | 4.06 | ||
aug-cc-pVQZ | 8.8 | 10.33 | 6.22 | 64.0 | 4.10 | ||
CCSD | aug-cc-pVDZ | 10.40 | 6.26 | 30.0 | 4.46 | ||
d-aug-cc-pVDZ | 8.49 | 5.12 | 47.6 | 4.47 | |||
aug-cc-pVTZ | 1.9 | 9.35 | 5.63 | 41.9 | 4.38 |
Wave function . | Basis . | CME . | Buckingham . | ||||
---|---|---|---|---|---|---|---|
. | . | . | . | . | . | ||
HF-SCF | aug-cc-pVDZ | 6.17 | 3.72 | 26.3 | 3.74 | ||
d-aug-cc-pVDZ | 3.46 | 2.09 | 33.6 | 3.81 | |||
aug-cc-pVTZ | 0.8 | 5.55 | 3.34 | 30.7 | 3.77 | ||
d-aug-cc-pVTZ | 3.9 | 5.35 | 3.22 | 34.0 | 3.64 | ||
aug-cc-pVQZ | 2.5 | 5.38 | 3.24 | 32.9 | 3.65 | ||
LDA | aug-cc-pVDZ | 13.70 | 8.25 | 45.7 | 4.31 | ||
d-aug-cc-pVDZ | 11.86 | 7.15 | 58.2 | 4.38 | |||
aug-cc-pVTZ | 5.3 | 12.82 | 7.72 | 53.3 | 4.47 | ||
d-aug-cc-pVTZ | 9.1 | 12.37 | 7.45 | 59.8 | 4.29 | ||
aug-cc-pVQZ | 7.0 | 12.46 | 7.50 | 57.1 | 4.32 | ||
B3LYP | aug-cc-pVDZ | 10.79 | 6.50 | 40.6 | 4.35 | ||
d-aug-cc-pVDZ | 8.78 | 5.29 | 51.8 | 4.39 | |||
aug-cc-pVTZ | 4.2 | 9.86 | 5.94 | 47.3 | 4.44 | ||
d-aug-cc-pVTZ | 7.6 | 9.44 | 5.69 | 53.2 | 4.26 | ||
aug-cc-pVQZ | 5.8 | 9.65 | 5.81 | 50.5 | 4.30 | ||
KT1 | aug-cc-pVDZ | 11.72 | 7.06 | 51.4 | 4.10 | ||
d-aug-cc-pVDZ | 9.15 | 5.51 | 64.5 | 4.14 | |||
aug-cc-pVTZ | 7.2 | 10.40 | 6.26 | 60.0 | 4.22 | ||
d-aug-cc-pVTZ | 11.0 | 10.13 | 6.10 | 66.9 | 4.06 | ||
aug-cc-pVQZ | 8.8 | 10.33 | 6.22 | 64.0 | 4.10 | ||
CCSD | aug-cc-pVDZ | 10.40 | 6.26 | 30.0 | 4.46 | ||
d-aug-cc-pVDZ | 8.49 | 5.12 | 47.6 | 4.47 | |||
aug-cc-pVTZ | 1.9 | 9.35 | 5.63 | 41.9 | 4.38 |
B3LYP GIAO magnetizability anisotropies are employed in place of the nonGIAO CCSD results.
The temperature-independent contribution to the CME is still significant, but yet on the average, when browsing through the results in Table IV, is only about 3% to 4% of (instead of the 20% given in Ref. 4).
The effect of electron correlation on is far more dramatic than seen in the original paper mostly due to the fact that the near cancellation of the paramagnetic and diamagnetic contribution yields an anisotropy of the hypermagnetizability that is far smaller than originally computed, and therefore far more sensitive to changes in the electron correlation treatment. Note that basis sets of double zeta quality yield the sign of opposite to that obtained with more extended (triple and quadruple zeta) basis sets. As a consequence of the reduced importance of the contribution of to and , we reduce our estimates given in Ref. 4 to , and , respectively. With the change of sign in also, the prediction of the Buckingham Effect (BE) constant and BE birefringence of changes, although only slightly: and , respectively. While the change in sign of modifies remarkably our revised value of the quadrupole moment of (see above), the revised value of at , given in the original paper as , does not change.
V.
CME and BE for at and . Atomic units, with in cgs units of . is the birefringence defined for an induction field of 1 T and a pressure of 1 atm according to Ref. 17. for Buckingham birefringence is given for a pressure of and EFG of .
Wave function . | Basis . | CME . | Buckingham . | ||||
---|---|---|---|---|---|---|---|
. | . | . | . | . | . | ||
HF | aug-cc-pVDZ | 6.93 | 4.17 | 207.4 | 8.52 | ||
d-aug-cc-pVDZ | 5.05 | 3.05 | 278.9 | 7.58 | |||
aug-cc-pVTZ | 7.10 | 4.28 | 260.2 | 7.97 | |||
d-aug-cc-pVTZ | 7.31 | 4.40 | 285.0 | 7.63 | |||
LDA-DFT | aug-cc-pVDZ | 13.13 | 7.91 | 247.6 | 4.82 | ||
d-aug-cc-pVDZ | 11.42 | 6.88 | 344.6 | 4.29 | |||
aug-cc-pVTZ | 13.48 | 8.12 | 311.7 | 4.07 | |||
d-aug-cc-pVTZ | 19.9 | 13.70 | 8.25 | 352.7 | 3.77 | ||
B3LYP-DFT | aug-cc-pVDZ | 10.64 | 6.41 | 236.3 | 7.87 | ||
d-aug-cc-pVDZ | 8.94 | 5.39 | 328.6 | 7.13 | |||
aug-cc-pVTZ | 10.57 | 6.37 | 296.9 | 7.15 | |||
d-aug-cc-pVTZ | 11.5 | 10.75 | 6.47 | 336.8 | 6.69 | ||
KT1-DFT | aug-cc-pVDZ | 8.97 | 5.40 | 261.9 | 4.26 | ||
d-aug-cc-pVDZ | 7.29 | 4.39 | 360.6 | 3.61 | |||
aug-cc-pVTZ | 8.92 | 5.37 | 325.9 | 4.02 | |||
d-aug-cc-pVTZ | 11.3 | 9.47 | 5.70 | 369.7 | 3.69 | ||
CCSDa | aug-cc-pVDZ | 8.73 | 5.26 | 227.6 | 8.12 | ||
d-aug-cc-pVDZ | 6.45 | 3.89 | 310.1 | 7.02 | |||
aug-cc-pVTZ | 8.78 | 5.29 | 284.7 | 6.42 | |||
Extrapolated from expt. | b | b |
Wave function . | Basis . | CME . | Buckingham . | ||||
---|---|---|---|---|---|---|---|
. | . | . | . | . | . | ||
HF | aug-cc-pVDZ | 6.93 | 4.17 | 207.4 | 8.52 | ||
d-aug-cc-pVDZ | 5.05 | 3.05 | 278.9 | 7.58 | |||
aug-cc-pVTZ | 7.10 | 4.28 | 260.2 | 7.97 | |||
d-aug-cc-pVTZ | 7.31 | 4.40 | 285.0 | 7.63 | |||
LDA-DFT | aug-cc-pVDZ | 13.13 | 7.91 | 247.6 | 4.82 | ||
d-aug-cc-pVDZ | 11.42 | 6.88 | 344.6 | 4.29 | |||
aug-cc-pVTZ | 13.48 | 8.12 | 311.7 | 4.07 | |||
d-aug-cc-pVTZ | 19.9 | 13.70 | 8.25 | 352.7 | 3.77 | ||
B3LYP-DFT | aug-cc-pVDZ | 10.64 | 6.41 | 236.3 | 7.87 | ||
d-aug-cc-pVDZ | 8.94 | 5.39 | 328.6 | 7.13 | |||
aug-cc-pVTZ | 10.57 | 6.37 | 296.9 | 7.15 | |||
d-aug-cc-pVTZ | 11.5 | 10.75 | 6.47 | 336.8 | 6.69 | ||
KT1-DFT | aug-cc-pVDZ | 8.97 | 5.40 | 261.9 | 4.26 | ||
d-aug-cc-pVDZ | 7.29 | 4.39 | 360.6 | 3.61 | |||
aug-cc-pVTZ | 8.92 | 5.37 | 325.9 | 4.02 | |||
d-aug-cc-pVTZ | 11.3 | 9.47 | 5.70 | 369.7 | 3.69 | ||
CCSDa | aug-cc-pVDZ | 8.73 | 5.26 | 227.6 | 8.12 | ||
d-aug-cc-pVDZ | 6.45 | 3.89 | 310.1 | 7.02 | |||
aug-cc-pVTZ | 8.78 | 5.29 | 284.7 | 6.42 | |||
Extrapolated from expt. | b | b |
B3LYP-DFT GIAO magnetizability anisotropies are employed in place of the nonGIAO CCSD results.
Data in Refs. 7 and 8 fitted in this work assuming linear regression with a value of equal to our B3LYP-DFT/daug-cc-pVTZ “best value” of
The temperature dependence of the CME of as computed here, compared to the experiment. The labels “aTZ” and “d-aTZ” stand for aug-cc-pVTZ and daug-cc-pVTZ, respectively.
The temperature dependence of the CME of as computed here, compared to the experiment. The labels “aTZ” and “d-aTZ” stand for aug-cc-pVTZ and daug-cc-pVTZ, respectively.
The temperature-independent CME contribution of is indeed similar to that of , as stated in Ref. 4, meaning that on the average, it is of the , never exceeding 10%. The value of 20% given originally is therefore overestimated.
The revision of the magnetizability anisotropy value of given by Lamb and Ritchie,7,8 revision made by employing the new best estimate for the hypermagnetizability anisotropy, the B3LYP-DFT/d-aug-cc-pVTZ value in Table V , now confirms not surprisingly the validity of the assumptions made by Lamb and Ritchie7,8 when they neglected the temperature-independent contribution to the CME. The first two sentences of the second paragraph on p. 114307–10 of Ref. 4 now read:
“Ritchie and Lamb in Refs. 7 and 8 neglected the temperature-independent contribution. We have rederived the value for by fitting their experimental data such that the line in Eq. (2) passes through our estimated intercept (B3LYP/daug-cc-pVTZ: ), assuming as Lamb and Ritchie a value for the anisotropy of the electric dipole polarizability of .9 The experimental estimate of of was confirmed.”
The last row of Table II of Ref. 4, with its old erroneously revised value of for of , should therefore be taken away together with the associated footnote.
The last two paragraphs of Sec. IV C 2 should be replaced by the following:
“We predict a Cotton–Mouton constant of with an associated birefringence of , under the conditions in Table V. With the exception of LDA-DFT, which overestimates the effect, we are inside the error bars of experiment and close to the center of the statistical distribution in particular with KT1-DFT and CCSD.
Finally, we predict the BE constant to be with an associated birefringence of . No experimental measurements are available for these constants.”
In the conclusions, Sec. V of Ref. 4, the last three paragraphs should be replaced by the following:
“The agreement with experimental data is satisfactory, particularly in view of the neglect of molecular vibrations. The temperature-independent contribution to the Cotton–Mouton birefringence is about 5% for both molecules, whereas the contributions to the BE are about 5% and 10% for and , respectively.
We have carried out a detailed and systematic investigation of the molecular quadrupole moment of both molecules, yielding and for and , respectively. For , this value is within one standard deviation of our revised experimental measurement, while for , our value supports the claims of Lamb and Ritchie7,8,10 that the measurement of Gierszal et al.11 is inaccurate.
Our best ab initio result for the magnetizability anisotropy of boron trichloride is in good agreement with the results of the measurements performed by Lamb and Ritchie,7,8 and excellent agreement is also observed between theory and experiment for the Cotton–Mouton constant of .”