共查询到20条相似文献,搜索用时 30 毫秒
1.
Xavier Prat-Resina Mireia Garcia-Viloca Gerald Monard Angels González-Lafont José M. Lluch Josep Maria Bofill Josep Maria Anglada 《Theoretical chemistry accounts》2002,107(3):147-153
We propose a methodology to locate stationary points on a quantum mechanical/molecular mechanical potential-energy surface.
This algorithm is based on a suitable approximation of an initial full Hessian matrix, either a modified Broyden–Fletcher–Goldfarg–Shanno
or a Powell update formula for the location of, respectively, a minimum or a transition state, and the so-called rational
function optimization. The latter avoids the Hessian matrix inversion required by a quasi-Newton–Raphson method. Some examples
are presented and analyzed.
Received: 16 July 2001 / Accepted: 9 October 2001 / Published online: 9 January 2002 相似文献
2.
Amy J. Austin Michael J. Frisch John A. Montgomery George A. Petersson 《Theoretical chemistry accounts》2002,107(3):180-186
An overlap criterion is defined that connects the identification of core orbitals in a molecular system, which can be problematic,
to that in isolated atoms, which is well defined. This approach has been tested on a variety of troublesome systems that have
been identified in the literature, including molecules containing third-row main-group elements, and is shown to remove errors
of up to 100 kcal/mol arising from an inconsistent treatment of core orbitals at different locations on a potential-energy
surface. For some systems and choices of core orbitals, errors as large as 19 kcal/mol can be introduced even when consistent
sets of orbitals are frozen, and the new method is shown to identify these cases of substantial core–valence mixing. Finally,
even when there is limited core–valence mixing, the frozen-core approximation can introduce errors of more than 5 kcal/mol,
which is much larger than the presumed accuracy of models such as G2 and CBS-QB3. The source of these errors includes interatomic
core–core and core–valence dispersion forces.
Received: 31 August 2001 / Accepted: 11 October 2001 / Published online: 9 January 2002 相似文献
3.
Kinetic isotope effects, KIEs, for hydrogen abstraction from C2H6 and C2D6 by chlorine atom have been studied by the dual-level direct dynamics approach. A low-level potential energy surface is obtained
with the MNDO-SRP method. High-level structural properties of the reactants, transition state, and products were obtained
at the MP2 level with the cc-pVDZ, aug-cc-pVDZ, and the cc-pVTZ basis sets. Using the variational transition state theory
with microcanonical optimized multidimensional tunneling, the values of deuterium KIE, at 300 K, range from 2.28 to 3.27,
in good agreement with the experimental values (2.69–5.88).
Received: 6 June 2001 / Accepted: 12 July 2001 / Published online: 19 November 2001 相似文献
4.
Rita Prosmiti Alexei A. Buchachenko Pablo Villarreal Gerardo Delgado-Barrio 《Theoretical chemistry accounts》2001,106(6):426-433
Ab initio molecular electronic structure calculations are performed for H5
+ at the QCISD(T) level of theory, using a correlation-consistent quadruple-zeta basis set. Structures, vibrational frequencies
and thermochemical properties are evaluated for ten stationary points of the H5
+ hypersurface and are compared with previous calculations. The features of the H3
+…H2 interaction at intermediate and large intermolecular distances are also investigated. Furthermore, an analytical functional
form for the potential-energy surface of H5
+ is derived using a first-order diatomics-in-molecule perturbation theory approach. Its topology is found to be qualitatively
correct for the short-range interaction region.
Received: 15 March 2001 / Accepted: 5 July 2001 / Published online: 11 October 2001 相似文献
5.
An efficient steepest descent algorithm for the integration of minimum energy paths, based on local quadratic approximations
of the potential energy surface, is presented. The algorithm incorporates a selection procedure for the points at which the
second derivatives of the energy are calculated fully or partially, thus minimizing the computational effort while maintaining
high accuracy. This makes the method especially well suited for application in variational transition state theory calculations
with tunnelling corrections, which have very high accuracy requirements. The performance of the algorithm is illustrated by
ab initio calculations for four chemical reactions of differing complexity. The overall computational cost is less than for,
or comparable to that of, first- or second-order algorithms published previously.
Received: 28 July 1998 / Accepted: 10 August 1998 / Published online: 28 October 1998 相似文献
6.
Ab initio calculations have been performed to study the molecular structures and vibrational levels of the four low-lying
ionic states (1, 22Π, and 1, 22Σ+) of carbonyl sulfide. The global regions of the potential-energy surfaces have been obtained by multireference single and
double excitation configuration interaction calculations. Vibrational calculations using explicit vibrational Hamiltonians
have been used for vibrational analysis. The equilibrium molecular structures and a vibrational analysis of the four ionic
states are presented. The theoretical ionization intensity curves including the vibrational structures of the ionic states
are also presented and are compared with the photoelectron spectrum.
Received: 20 January 2001 / Accepted: 22 August 2001 / Published online: 30 October 2001 相似文献
7.
The crystalline transformation in nylon 12 12 was monitored by variable-temperature wide-angle X-ray diffraction during heating,
isothermal and cooling processes. It could be found that the α phase of the sample transforms into a γ phase at about 130 °C
if the sample is heated from room temperature to a high temperature, which is the so-called Brill transition temperature.
In addition, nylon 12 12 was found to crystallize into the γ phase when isothermally crystallized at high temperature from
the melt. Upon further cooling from the crystallization temperature to room temperature, the reverse transition from the γ phase
to the α phase can be detected. Furthermore, in situ Fourier transform IR spectroscopy was also used to study the Brill transition
of nylon 12 12 samples on both heating and cooling. It is interesting that the hydrogen-bond strength in nylons decreases
dramatically and that some Brill bands of nylon 12 12 disappear abruptly during the Brill transition on heating.
Received: 3 April 2001 Accepted: 28 July 2001 相似文献
8.
Fernando R. Ornellas 《Theoretical chemistry accounts》2000,103(6):469-476
The [H,S,Cl] potential-energy surface has been investigated at the self-consistent field (SCF), complete active space self-consistent
field (CASSCF), second-order M?ller–Plesset, coupled-cluster single-double and perturbative triple excitation, [CCSD(T)]/6-31G(d,p),
6-31G(2df,2pd), and correlation-consistent polarized valence triple zeta (cc-pVTZ) levels of theory. CCSD(T)/ cc-pVTZ results
predict a very stable HSCl species, an isomer HClS, 51.84 kcal/mol higher in energy, and a transition state 57.68 kcal/mol
above HSCl. Independent of the level of theory, results with the smaller 6-31G(d,p) basis set turned out to be poor, especially
for HClS. Vibrational analysis indicates that both species can be easily differentiated if isolated. Bonding differences between
these molecules are illustrated by contour plots of valence orbitals. Viewed classically, bonding in HClS involves a dative
bond. Transition-state rate constants, and equilibrium constants for the HSCl ↔ HClS isomerization have been estimated for
various temperatures (200–1000 K). At 298.15 K, the forward rate is predicted to be 7.95 × 10−29 s−1, and the equilibrium constant to be 2.31 × 10−38. Tunneling corrections vary from 1.57 at 298.15 K to 1.05 at 1000 K. Activation energies have been obtained by a two-points
linear fit to the Arrhenius equation.
Received: 7 May 1999 / Accepted: 22 July 1999 / Published online: 4 October 1999 相似文献
9.
10.
The surface tensions of aqueous poly(N-vinylformamide) (PNVF), polyvinylamine (PVAM), and PNVF–PVAM copolymers were measured as functions of pH. The nonionic PNVF
gave a pH-independent surface tension of 68 mN/m. The surface tension of PVAM was pH dependent; at pH 10 it was 56 mN/m, whereas
it was 71.5 mN/m at pH 3.5. The transition from higher to lower surface tension values occurred most dramatically between
pH 8 and 9, reflecting the dissociation behavior of the amine groups. The copolymers showed intermediate behavior.
Received: 20 August 2001 Accepted: 26 September 2001 相似文献
11.
Jing-yao Liu Ze-sheng Li Zhen-wen Dai Xu-ri Huang Chia-chung Sun 《Theoretical chemistry accounts》2002,108(3):179-186
The hydrogen-abstraction reaction C2H5F+O → C2H4F+OH has been studied by a dual-level direct dynamics method. For the reaction, three reaction channels, one for α-abstraction
and two for β-abstraction, have been identified. The potential-energy surface information is obtained at the MP2(full)/6-311G(d,p)
and PMP2(full)/6-311G(3df,3pd) (single-point) levels. By canonical variational transition-state theory, rate constants for
each reaction channel are calculated with a small-curvature tunneling correction. The total rate constant is calculated from
the sum of the individual rate constants and the temperature dependence of the branching ratios is obtained over a wide range
of temperatures from 300 to 5,000 K. The agreement of the rate constants with experiment is good in the experimental temperature
range from 1,000 to 1,250 K. The calculated results indicate that at low temperatures α-abstraction is most likely to be the
major reaction channel, while β-abstraction channels will significantly contribute to the whole reaction rate as the temperature
increases.
Received: 23 January 2002 / Accepted: 23 June 2002 / Published online: 20 September 2002 相似文献
12.
Hybrid quantum mechanical/molecular mechanical electronic structure calculations reveal the transition state for C–H bond
cleavage in [(LCu)2 (μ-O)2]2+ (L=1,4,7-tribenzyl-1,4,7-triazacyclononane) to be consistent with a hydrogen-atom-transfer mechanism from carbon to oxygen.
At the MPW1K/double-zeta effective core potential(+)|univeral force field level, 0 K activation enthalpies for the parent,
p-CF3, and p-OH substituted benzyl systems are predicted to be 8.8, 9.5, and 7.8 kcal/mol. Using a one-dimensional Eckart potential to
estimate quantum effects on the reaction coordinate, reaction in the unsubstituted system is predicted to proceed with a primary
kinetic isotope effect of 22 at 233 K. Structural parameters associated with the hydrogen-atom transfer are consistent with
the Hammond postulate.
Received: 10 October 2000 / Accepted: 3 November 2000 / Published online: 3 April 2001 相似文献
13.
M. C. Bacchus-Montabonel 《Theoretical chemistry accounts》2000,104(3-4):296-301
Ab initio potential-energy curves and coupling matrix elements of the Σ and Π molecular states involved in the collision
of the Si2+, Si3+ and Si4+ multicharged ions on atomic hydrogen and helium have been determined by means of configuration interaction methods. The total
and partial electron capture cross sections have been determined using a semiclassical or a quantal approach in the 0.002–0.1 au
velocity range. A detailed comparison with very recent theoretical and experimental rate coefficient results is made.
Received: 14 September 1999 / Accepted: 3 February 2000 / Published online: 2 May 2000 相似文献
14.
Kazunari Naka Hirofumi Sato Akihiro Morita Fumio Hirata Shigeki Kato 《Theoretical chemistry accounts》1999,102(1-6):165-169
The free-energy profile for the Menshutkin-type reaction NH3 + CH3Cl → NH3CH3
+ + Cl− in aqueous solution is studied using the RISM-SCF method. The effect of electron correlation on the free-energy profile is
estimated by the RISM-MP2 method at the HF optimized geometries along the reaction coordinate. Solvation was found to have
a large influence on the vibrational frequencies at the reactant, transition state and product; these vibrational frequencies
are utilized to calculate the zero-point energy correction of the free-energy profile. The computed barrier height and reaction
exothermicity are in reasonable agreement with those of experiment and previous calculations. The change of solvation structure
along the reaction path is represented by radial distribution functions between solute-solvent atomic sites. The mechanisms
of the reaction are discussed from the view points of solute electronic and solvation structures.
Received: 26 June 1998/Accepted: 28 August 1998 / Published online: 2 November 1998 相似文献
15.
The rational function optimization algorithm is one of the widely used methods to search stationary points on surfaces. However,
one of the drawbacks of this method is the step reduction procedure to deal with the overstepping problem. We present and
comment on a method such that the step obtained from the solution of the rational function equations possesses the desired
correct length. The analysis and discussion of the method is mainly centered on the location and optimization of transition
states.
Received: 18 June 1998 / Accepted: 17 September 1998 / Published online: 23 November 1998 相似文献
16.
P. L. A. Popelier 《Theoretical chemistry accounts》2001,105(4-5):393-399
This work proposes a novel algorithm to compute atomic charges as defined by the theory of “atoms in molecules” (AIM). Using
the divergence theorem it is possible to express the 3D volume integral over an atomic basin purely in terms of 2D surface
integrals. Hence, it can be proven that an atomic charge is equal to the flux of the electric field of the whole molecule
through the atom's complete boundary. This boundary consists of the interatomic surfaces and the so-called outeratomic surface,
which is the open side of the atom. When fine-tuned the algorithm can generate atomic charges in the order of minutes without
introducing any approximations. Moreover, the problem of the geometrical cusp occurring in atomic basins and that of multiple
intersections is also eliminated. The computational overhead of computing the electric field (which is analytical) is compensated
by the gain in computing time by eliminating one dimension of quadrature. The proposed algorithm opens an avenue to invalidate
the oft-quoted drawback that AIM charges are computationally expensive. We explain the details of the implementation in MORPHY01
and illustrate the novel algorithm with a few examples.
Received: 1 June 2000 / Accepted: 4 October 2000 / Published online: 23 January 2001 相似文献
17.
The B3LYP/6-311G(d) and CCSD(T)/6-311G(2df) (single-point) methods have been used to investigate the singlet potential energy
surface of C2NP, in which seven stationary isomers and seventeen interconversion transition states are involved. At the final CCSD(T)/6-311G(2df)//B3LYP6-311G(d)
level with zero-point vibrational energy correction the lowest-lying isomer is a linear NCCP followed by two linear CNCP isomers
at 23.9 and CCNP at 65.8 kcal mol−1, respectively. The three isomers are kinetically very stable towards both isomerization and dissociation, and CCNP is even
more kinetically stable than CNCP – by 14.3 kcal mol−1 despite its high energy. Further comparative calculations were performed at the QCISD and QCISD(T) levels with the 6-311G(d)
and 6-311G(2d) basis sets to obtain more reliable structures and spectroscopy for the three isomers. The calculated bond lengths,
rotational constant, and dipole moment for NCCP were in close agreement with the experimentally determined values. Finally,
similarities and discrepancies between the potential energy surface of C2NP and those of the analogous species C2N2 and C2P2 were compared. The results presented in this paper might be helpful for future identification of the two still unknown yet
kinetically very stable isomers CNCP and CCNP, both in the laboratory and in interstellar space.
Received: 3 January 2001 / Accepted: 6 June 2001 / Published online: 30 October 2001 相似文献
18.
Alexandra Romina Albunia Raffaele Borrelli Andrea Peluso 《Theoretical chemistry accounts》2000,104(3-4):218-222
Electron transfer (ET) from toluene to the nitronium ion in the region of van der Waals intermolecular distances has been
investigated by a quantum dynamical analysis performed on potential-energy surfaces computed at the ab initio multireference
configuration interaction level. The results show that ET is very fast, occurring on a timescale of a few picoseconds. This
has important implications for the mechanism of aromatic nitration: the ET path can compete efficiently with the direct attack
of the nitronium ion to the aromatic substrate to yield the Wheland intermediate and that could explain some unsettled points
in the mechanism of aromatic nitration.
Received: 16 September 1999 / Accepted: 3 February 2000 / Published online: 12 May 2000 相似文献
19.
A method is described to extract a complete set of sequence-dependent energy parameters for a rigid base-pair model of DNA
from molecular dynamics (MD) simulations. The method is properly consistent with equilibrium statistical mechanics and leads
to effective inertia parameters for the base-pair units as well as stacking and stiffness parameters for the base-pair junctions.
We give explicit formulas that yield a complete set of base-pair model parameters in terms of equilibrium averages that can
be estimated from a time series generated in an MD simulation. The expressions to be averaged depend strongly both on the
choice of coordinates used to describe rigid-body orientations and on the choice of strain measures at each junction.
Received: 12 July 2000 / Accepted: 5 January 2001 / Published online: 3 May 2001 相似文献
20.
The electronic structure of the TiSi molecule was examined using two types of multireference single and double excitation
configuration interactions with highly extended basis sets, one including valence correlation and the other including valence
and core–valence correlation. A multireference coupled-pair approximation (MRCPA) was further applied to the latter. The calculations
suggest a 5Δ ground state, and the lowest excited state is 3Π and is only slightly (0.12 eV as estimated by MRCPA) above the ground state. The spectroscopic constants of the low-lying
1Δ, 3Δ, 1Π, 5Π, and 7Σ+ states as well as the 5Δ ground state and the 3Π excited states were evaluated, and we found that the molecule has only a weak σ bond and that six of the eight valence electrons
essentially do not contribute to the bonding. The bonding nature of TiSi in these states is discussed in comparison with the
TiC molecule.
Received: 7 October 2000 / Accepted: 8 January 2001 / Published online: 3 May 2001 相似文献