PubMed Health⌕ Search

Biomedical subjects

Nathan J DeYonker

Publications and source records attributed to Nathan J DeYonker.

6 recordsLinked to original sources

The correlation-consistent composite approach: application to the G3/99 test set.

The correlation-consistent composite approach (ccCA), an ab initio composite technique for computing atomic and molecular energies, recently has been shown to successfully reproduce experimental data for a number of systems. The ccCA is applied to the G3/99 test set, which includes 223 enthalpies of formation, 88 adiabatic ionization potentials, 58 adiabatic electron affinities, and 8 adiabatic proton affinities. Improvements on the original ccCA formalism include replacing the small basis set quadratic configuration interaction computation with a coupled cluster computation, employing a correction for scalar relativistic effects, utilizing the tight-d forms of the second-row correlation-consistent basis sets, and revisiting the basis set chosen for geometry optimization. With two types of complete basis set extrapolation of MP2 energies, ccCA results in an almost zero mean deviation for the G3/99 set (with a best value of -0.10 kcal mol(-1)), and a 0.96 kcal mol(-1) mean absolute deviation, which is equivalent to the accuracy of the G3X model chemistry. There are no optimized or empirical parameters included in the computation of ccCA energies. Except for a few systems to be discussed, ccCA performs as well as or better than Gn methods for most systems containing first-row atoms, while for systems containing second-row atoms, ccCA is an improvement over Gn model chemistries.

Journal Article↗

Accurate enthalpies of formation of alkali and alkaline earth metal oxides and hydroxides: assessment of the correlation consistent composite approach (ccCA).

Computing the enthalpies of formation for alkali metal and alkaline earth metal oxides (M(x)O) and hydroxides [M(OH)(n)] using the Gaussian-n (Gn) and Weismann-n (Wn) ab initio model chemistries is difficult due to an improper treatment of core-valence electron correlation effects. Using a new model chemistry called the correlation consistent Composite Approach (ccCA), enthalpies of formation were determined for eight different alkali/alkaline earth metal oxides and hydroxides. Unlike the Gn and Wn model chemistries, which must be modified to properly account for core-valence electron correlation, the standard implementations of the ccCA provide acceptable results, and all enthalpies of formation obtained with the ccCA are within the accepted range of recommended values.

Computer Simulation↗

The correlation consistent composite approach (ccCA): an alternative to the Gaussian-n methods.

An alternative to the Gaussian-n (G1, G2, and G3) composite methods of computing molecular energies is proposed and is named the "correlation consistent composite approach" (ccCA, ccCA-CBS-1, ccCA-CBS-2). This approach uses the correlation consistent polarized valence (cc-pVXZ) basis sets. The G2-1 test set of 48 enthalpies of formation (DeltaHf), 38 adiabatic ionization potentials (IPs), 25 adiabatic electron affinities (EAs), and 8 adiabatic proton affinities (PAs) are computed using this approach, as well as the DeltaHf values of 30 more systems. Equilibrium molecular geometries and vibrational frequencies are obtained using B3LYP density functional theory. When applying the ccCA-CBS method with the cc-pVXZ series of basis sets augmented with diffuse functions, mean absolute deviations within the G2-1 test set compared to experiment are 1.33 kcal mol(-1) for DeltaHf,0.81 kcal mol(-1) for IPs, 1.02 kcal mol(-1) for EAs, and 1.51 kcal mol(-1) for PAs, without including the "high-level correction" (HLC) contained in the original Gn methods. Whereas the HLC originated in the Gaussian-1 method as an isogyric correction, it evolved into a fitted parameter that minimized the error of the composite methods, eliminating its physical meaning. Recomputing the G1 and G3 enthalpies of formation without the HLC reveals a systematic trend where most DeltaHf values are significantly higher than experimental values. By extrapolating electronic energies to the complete basis set (CBS) limit and adding G3-like corrections for the core-valence and infinite-order electron correlation effects, ccCA-CBS-2 often underestimates the experimental DeltaHf, especially for larger systems. This is desired as inclusion of relativistic and atomic spin-orbit effects subsequently improves theoretical DeltaHf values to give a 0.81 kcal mol(-1) mean absolute deviation with ccCA-CBS-2. The ccCA-CBS method is a viable "black box" method that can be used on systems with at least 10-15 heavy atoms.

Journal Article↗

The ground and two lowest-lying singlet excited electronic states of copper hydroxide (CuOH).

Various ab initio methods, including self-consistent field (SCF), configuration interaction, coupled cluster (CC), and complete-active-space SCF (CASSCF), have been employed to study the electronic structure of copper hydroxide (CuOH). Geometries, total energies, dipole moments, harmonic vibrational frequencies, and zero-point vibrational energies are reported for the linear 1Sigma+ and 1Pi stationary points, and for the bent ground-state X 1A', and excited-states 2 1A' and 1 1A". Six different basis sets have been used in the study, Wachters/DZP being the smallest and QZVPP being the largest. The ground- and excited-state bending modes present imaginary frequencies for the linear stationary points, indicating that bent structures are more favorable. The effects of relativity for CuOH are important and have been considered using the Douglas-Kroll approach with cc-pVTZ/cc-pVTZ_DK and cc-pVQZ/cc-pVQZ_DK basis sets. The bent ground and two lowest-lying singlet excited states of the CuOH molecule are indeed energetically more stable than the corresponding linear structures. The optimized geometrical parameters for the X 1A' and 1 1A" states agree fairly well with available experimental values. However, the 2 1A' structure and rotational constants are in poor agreement with experiment, and we suggest that the latter are in error. The predicted adiabatic excitation energies are also inconsistent with the experimental values of 45.5 kcal mol(-1) for the 2 1A' state and 52.6 kcal mol(-1) for the 1 1A" state. The theoretical CC and CASSCF methods show lower adiabatic excitation energies for the 1 1A" state (53.1 kcal mol(-1)) than those for the corresponding 2 1A' state (57.6 kcal mol(-1)), suggesting that the 1 1A" state might be the first singlet excited state while the 2 1A' state might be the second singlet excited state.

Journal Article↗

Application of equation-of-motion coupled-cluster methods to low-lying singlet and triplet electronic states of HBO and BOH.

The equilibrium structures and physical properties of the X (1)sigma(+) linear electronic states, linear excited singlet and triplet electronic states of hydroboron monoxide (HBO) (A (1)sigma(-), B (1)delta, a (3)sigma(+), and b (3)delta) and boron hydroxide (BOH) (A (1)sigma(+), B (1)Pi, and b (3)Pi), and their bent counterparts (HBO a (3)A('), b (3)A("), A (1)A("), B (1)A(') and BOH X (1)A('), b (3)A('), c (3)A("), A (1)A('), B (1)A('), C (1)A(")) are investigated using excited electronic state ab initio equation-of-motion coupled-cluster (EOM-CC) methods. A new implementation of open-shell EOM-CC including iterative partial triple excitations (EOM-CC3) was tested. Coupled-cluster wave functions with single and double excitations (CCSD), single, double, and iterative partial triple excitations (CC3), and single, double, and full triple excitations (CCSDT) are employed with the correlation-consistent quadruple and quintuple zeta basis sets. The linear HBO X (1)sigma(+) state is predicted to lie 48.3 kcal mol(-1) (2.09 eV) lower in energy than the BOH X (1)sigma(+) linear stationary point at the CCSDT level of theory. The CCSDT BOH barrier to linearity is predicted to lie 3.7 kcal mol(-1) (0.16 eV). With a harmonic zero-point vibrational energy correction, the HBO X (1)sigma(+)-BOH X (1)A(') energy difference is 45.2 kcal mol(-1) (1.96 eV). The lowest triplet excited electronic state of HBO, a (3)A('), has a predicted excitation energy (T(e)) of 115 kcal mol(-1) (4.97 eV) from the HBO ground state minimum, while the lowest-bound BOH excited electronic state, b (3)A('), has a T(e) of 70.2 kcal mol(-1) (3.04 eV) with respect to BOH X (1)A('). The T(e) values predicted for the lowest singlet excited states are A (1)A(")<--X (1)sigma(+)=139 kcal mol(-1) (6.01 eV) for HBO and A (1)A(')<--X (1)A(')=102 kcal mol(-1) (4.42 eV) for BOH. Also for BOH, the triplet vertical transition energies are b (3)A(')<--X (1)A(')=71.4 kcal mol(-1) (3.10 eV) and c (3)A(")<--X (1)A(')=87.2 kcal mol(-1) (3.78 eV).

Journal Article↗

Low-lying electronic states of FeNC and FeCN: a theoretical journey into isomerization and quartet/sextet competition.

With several levels of multireference and restricted open-shell single-reference electronic structure theory, optimum structures, relative energetics, and spectroscopic properties of the low-lying (6)Delta, (6)Pi, (4)Delta, (4)Pi, and (4)Sigma(-) states of linear FeNC and FeCN have been investigated using five contracted Gaussian basis sets ranging from Fe[10s8p3d], C/N[4s2p1d] to Fe[6s8p6d3f2g1h], C/N[6s5p4d3f2g]. Based on multireference configuration interaction (MRCISD+Q) results with a correlation-consistent polarized valence quadruple-zeta (cc-pVQZ) basis set, appended with core correlation and relativistic corrections, we propose the relative energies: T(e)(FeNC), (6)Delta(0)<(6)Pi (2300 cm(-1))<(4)Delta (2700 cm(-1))<(4)Pi (4200 cm(-1))<(4)Sigma(-); and T(e)(FeCN), (6)Delta(0)<(6)Pi (1800 cm(-1))<(4)Delta (2500 cm(-1))<(4)Pi (2900 cm(-1))<(4)Sigma(-). The (4)Delta and (4)Pi states have massive multireference character, arising mostly from 11sigma-->12sigma promotions, whereas the sextet states are dominated by single electronic configurations. The single-reference CCSDT-3 (coupled cluster singles and doubles with iterative partial triples) method appears to significantly overshoot the stabilization of the quartet states provided by both static and dynamical correlation. The (4,6)Delta and (4,6)Pi states of both isomers are rather ionic, and all have dipole moments near 5 D. On the ground (6)Delta surface, FeNC is predicted to lie 0.6 kcal mol(-1) below FeCN, and the classical barrier for isocyanide/cyanide isomerization is about 6.5 kcal mol(-1). Our data support the recent spectroscopic characterization by Lei and Dagdigian [J. Chem. Phys. 114, 2137 (2000)] of linear (6)Delta FeNC as the first experimentally observed transition-metal monoisocyanide. Their assignments for the ground term symbol, isotopomeric rotational constants, and the Fe-N omega(3) stretching frequency are confirmed; however, we find rather different structural parameters for (6)Delta FeNC:r(e)(Fe-N)=1.940 A and r(N-C)=1.182 A at the cc-pVQZ MRCISD+Q level. Our results also reveal that the observed band of FeNC originating at 27 236 cm(-1) should have an analog in FeCN near 23 800 cm(-1) of almost equal intensity. Therefore, both thermodynamic stability and absorption intensity factors favor the eventual observation of FeCN via a (6)Pi<--(6)Delta transition in the near-UV.

Journal Article↗