Use of L matrix to obtain reliable ab initio force constants of polyatomic molecules: ethylene as a test system

Direct use of experimental eigenvalues of the vibrational secular equation on to the ab initio predicted eigenvector space is suggested as a means of obtaining a reliable set of intramolecular force constants. This method which we have termed RECOVES (recovery in the eigenvector space) is computationally simple and free from arbitrariness. The RECOVES force constants, by definition, reproduce the experimental vibrational frequencies of the parent molecule exactly. The ab initio calculations were carried out for ethylene as a test molecule and the force constants obtained by the present procedure also correctly predict the vibrational frequencies of the deuterated species. The RECOVES force constants for ethylene are compared with those obtained by using the SQM procedure.

Reactions of cyanogen radical with alkanes and an explanation for negative temperature dependence of rate constants

Reactions of cyanide radicals with alkanes have been investigated by ab initio methods. It is found that the potential energy surface for reaction of CN with a primary C-H bond in methane has a small positive barrier while reactions of CN with a secondary and a tertiary C-H bond in alkanes are barrierless at the correlated level. A simple explanation for the obtained negative temperature dependence of rate constants for reactions of CN with a secondary and a tertiary C-H bond in alkanes are given in terms of the collision theory of bimolecular reactions. It is shown that for barrierless reactions the negative temperature dependence of the rate constants is attributed to the variation of the pre-exponential factor with temperature.

SCF and electron correlation studies on structure, force constants and vibrational spectra of borane monoammoniate complex

Effects of basis set and electron correlation on the equilibrium geometry, force constants and vibrational spectra of BH3NH3 have been studied. A series of basis sets ranging from double zeta to triple zeta including polarization and diffuse functions have been utilized. All the SCF based calculations overestimate the dative B-N bond distance and considerable improvement occurs when the treatment for electron correlation is introduced. Detailed vibrational analysis for BH3NH3 has been carried out. The mean absolute percentage deviation of the ab initio predicted vibration frequencies of (BH3NH3)-B-11 from the experiment is about 10% for the SCF based calculations and the MP2 method shows better agreement, the overall deviation being 5-6%. The ground state effective force constants of BH3NH3 were obtained using RECOVES procedure. The RECOVES sets of force constants are found to be highly satisfactory for the prediction of the vibrational frequencies of different isotopomers of BH3NH3. The mean absolute percentage deviation of the calculated frequencies of different isotopomers from the experiment is much less than 1%. The RECOVES-MP2/augDZP set of force constants was found to be the best set among the different sets for this molecule. Theoretical infrared intensities are in fair agreement with the observed spectral features.

Dissociation constants of weak organic acids in protic solvents obtained from their first hyperpolarizabilities in solution

The first hyperpolarizabilities (beta) of some weak aromatic organic acids have been measured in protic solvents by the hyper-Rayleigh scattering (HRS) technique at low concentrations. The measured hyperpolarizability (beta(m)) varies between the two extreme limits: the hyperpolarizability of the acid form (beta(HA)) at the lower side and that of the basic form (beta(A-)) at the higher side. The degree of dissociation (alpha) of the acid in a solvent is related to the measured hyperpolarizability, beta(m), by the following relationship: beta(m)(2)=(1-alpha)beta(HA)(2)+alpha beta(A-)(2). The calculated beta's including solvent effects in terms of an Onsager field do not reproduce the experimentally measured hyperpolarizabilities. Other solvent-induced effects like hydrogen bonding and van der Waals interactions seem to influence the first hyperpolarizability and, thus, indirectly the extent of dissociation of these weak acids in these protic solvents.

Interplay of spin and lattice degrees of freedom in the frustrated antiferromagnet CdCr2O4: High-field and temperature-induced anomalies of the elastic constants

Temperature and magnetic field studies of the elastic constants of the chromium spinel CdCr2O4 show pronounced anomalies related to strong spin-phonon coupling in this frustrated antiferromagnet. A detailed comparison of the longitudinal acoustic mode propagating along the 111] direction with a theory based on an exchange-striction mechanism leads to an estimate of the strength of the magnetoelastic interaction. The derived spin-phonon coupling constant is in good agreement with previous determinations based on infrared absorption. Further insight is gained from intermediate and high magnetic field experiments in the field regime of the magnetization plateau. The role of the antisymmetric Dzyaloshinskii-Moriya interaction is discussed.

Formulas for algorithmic computation of impedance matrix using general circuit constants

A set of formulas is derived from general circuit constants which facilitates formation of the impedance matrix of a power system by the bus-impedance method. The errors associated with the lumpedparameter representation of a transmission line are thereby eliminated. The formulas are valid for short lines also, if the relevant general circuit constants are employed. The mutual impedance between the added line and the existing system is not considered, but the approach suggested can well be extended to it.

Structure constants of beta deformed super Yang-Mills

We study the structure constants of the N = 1 beta deformed theory perturbatively and at strong coupling. We show that the planar one loop corrections to the structure constants of single trace gauge invariant operators in the scalar sector is determined by the anomalous dimension Hamiltonian. This result implies that 3 point functions of the chiral primaries of the theory do not receive corrections at one loop. We then study the structure constants at strong coupling using the Lunin-Maldacena geometry. We explicitly construct the supergravity mode dual to the chiral primary with three equal U(1) R-charges in the Lunin-Maldacena geometry. We show that the 3 point function of this supergravity mode with semi-classical states representing two other similar chiral primary states but with large U(1) charges to be independent of the beta deformation and identical to that found in the AdS(5) x S-5 geometry. This together with the one-loop result indicate that these structure constants are protected by a non-renormalization theorem. We also show that three point function of U(1) R-currents with classical massive strings is proportional to the R-charge carried by the string solution. This is in accordance with the prediction of the R-symmetry Ward identity.

Dissociation Constants of Weak Acids from ab Initio Molecular Dynamics Using Metadynamics: Influence of the Inductive Effect and Hydrogen Bonding on pK(a) Values

The theoretical estimation of the dissociation constant, or pK(a), of weak acids continues to be a challenging field. Here, we show that ab initio CarParrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free-energy profile of the dissociation reaction provide reasonable estimates of the pK(a) value. Water molecules, sufficient to complete the three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. The free-energy profiles exhibit two distinct minima corresponding to the dissociated and neutral states of the acid, and the difference in their values provides the estimate for pK(a). We show for a series of organic acids that CPMD simulations in conjunction with metadynamics can provide reasonable estimates of pK(a) values. The acids investigated were aliphatic carboxylic acids, chlorine-substituted carboxylic acids, cis- and trans-butenedioic acid, and the isomers of hydroxybenzoic acid. These systems were chosen to highlight that the procedure could correctly account for the influence of the inductive effect as well as hydrogen bonding on pK(a) values of weak organic acids. In both situations, the CPMD metadynamics procedure faithfully reproduces the experimentally observed trend and the magnitudes of the pK(a) values.

Isolated singularities of polyharmonic operator in even dimension

We consider the equation Delta(2)u = g(x, u) >= 0 in the sense of distribution in Omega' = Omega\textbackslash {0} where u and -Delta u >= 0. Then it is known that u solves Delta(2)u = g(x, u) + alpha delta(0) - beta Delta delta(0), for some nonnegative constants alpha and beta. In this paper, we study the existence of singular solutions to Delta(2)u = a(x) f (u) + alpha delta(0) - beta Delta delta(0) in a domain Omega subset of R-4, a is a nonnegative measurable function in some Lebesgue space. If Delta(2)u = a(x) f (u) in Omega', then we find the growth of the nonlinearity f that determines alpha and beta to be 0. In case when alpha = beta = 0, we will establish regularity results when f (t) <= Ce-gamma t, for some C, gamma > 0. This paper extends the work of Soranzo (1997) where the author finds the barrier function in higher dimensions (N >= 5) with a specific weight function a(x) = |x|(sigma). Later, we discuss its analogous generalization for the polyharmonic operator.

Optical Factors Determined by the T-matrix Method in Turbidity Measurement of Absolute Coagulation Rate Constants

Turbidity measurement for the absolute coagulation rate constant of suspensions has been extensively adopted because of its simplicity and easy implementation. A key factor to derive the rate constant from experimental data is how to theoretically evaluate the so-called optical factor involved in calculating the extinction cross section of doublets formed in the aggregation. In a previous paper, we have shown that compared with other theoretical approaches, the T-matrix method provides a robust solution to this problem and is effective in extending the applicability range of the turbidity methodology as well as increasing measurement accuracy. This paper will provide a more comprehensive discussion about the physical insight of using the T-matrix method in turbidity measurement and associated technical details. In particular, the importance of ensuring the correct value for the refractive indices for colloidal particles and the surrounding medium used in the calculation is addressed because the indices generally vary with the wavelength of the incident light. The comparison of calculated results with experiments shows that the T-matrix method can correctly calculate optical factors even for large particles, whereas other existing theories cannot. In addition, the calculated data of the optical factor by the T-matrix method for a range of particle radii and incident light wavelengths are listed.

Toward an understanding of the turbidity measurement of fleterocoagulation rate constants of dispersions containing particles of different sizes

Our previous studies have shown that the determination of coagulation rate constants by turbidity measurement becomes impossible for a certain operating wavelength (that is, its blind point) because at this wavelength the change in the turbidity of a dispersion completely loses its response to the coagulation process. Therefore, performing the turbidity measurement in the wavelength range near the blind point should be avoided. In this article, we demonstrate that the turbidity measurement of the rate constant for coagulation of a binary dispersion containing particles of two different sizes (heterocoagulation) presents special difficulties because the blind point shifts with not only particle size but also with the component fraction. Some important aspects of the turbidity measurement for the heterocoagulation rate constant are discussed and experimentally tested. It is emphasized that the T-matrix method can be used to correctly evaluate extinction cross sections of doublets formed during the heterocoagulation process, which is the key data determining the rate constant from the turbidity measurement, and choosing the appropriate operating wavelength and component fraction are important to achieving a more accurate rate constant. Finally, a simple scheme in experimentally determining the sensitivity of the turbidity changes with coagulation over a wavelength range is proposed.

A new expression of hardening coefficients for fcc-crystal and calibration of the material constants

In order to describe the effect of latent hardening on the macro-plastic behavior of foc-crystal, a new expression for hardening coefficient is proposed in which there are 12 material constants, each having clear physical meaning. And a method of material constant calibration is suggested and used to determine the material constants of copper and aluminum crystal. The simulated load-elongation curves along various crystallographic orientations are comparable with the experimental ones.

Mechanisms underlying two kinds of surface effects on elastic constants

ABSTRACT Recently, people are confused with two opposite variations of elastic modulus with decreasing size of nano scale sample: elastic modulus either decreases or increases with decreas- ing sample size. In this paper, based on intermolecular potentials and a one dimensional model, we provide a unified understanding of the two opposite size effects. Firstly, we analyzed the mi- crostructural variation near the surface of an fcc nanofilm based on the Lennard-Jones potential. It is found that the atomic lattice near the surface becomes looser in comparison with the bulk, indicating that atoms in the bulk are located at the balance of repulsive forces, resulting in the decrease of the elastic moduli with the decreasing thickness of the film accordingly. In addition, the decrease in moduli should be attributed to both the looser surface layer and smaller coor- dination number of surface atoms. Furthermore, it is found that both looser and tighter lattice near the surface can appear for a general pair potential and the governing mechanism should be attributed to the surplus of the nearest force to all other long range interactions in the pair po- tential. Surprisingly, the surplus can be simply expressed by a sum of the long range interactions and the sum being positive or negative determines the looser or tighter lattice near surface re- spectively. To justify this concept, we examined ZnO in terms of Buckingham potential with long range Coulomb interactions. It is found that compared to its bulk lattice, the ZnO lattice near the surface becomes tighter, indicating the atoms in the bulk located at the balance of attractive forces, owing to the long range Coulomb interaction. Correspondingly, the elastic modulus of one- dimensional ZnO chain increases with decreasing size. Finally, a kind of many-body potential for Cu was examined. In this case, the surface layer becomes tighter than the bulk and the modulus increases with deceasing size, owing to the long range repulsive pair interaction, as well as the cohesive many-body interaction caused by the electron redistribution.