Perturbative growth of cosmological clustering. I: Formalism

Here we rederive the hierarchy of equations for the evolution of distribution functions of various orders using a convenient parameterization. We use this to obtain equations for two- and three-point correlation functions in powers of a small parameter, viz., the initial density contrast. The correspondence of the lowest order solutions of these equations to the results from the linear theory of density perturbations is shown for an OMEGA = 1 universe. These equations are then used to calculate, to the lowest order, the induced three-point correlation function that arises from Gaussian initial conditions in an OMEGA = 1 universe. We obtain an expression which explicitly exhibits the spatial structure of the induced three-point correlation function. It is seen that the spatial structure of this quantity is independent of the value of OMEGA. We also calculate the triplet momentum. We find that the induced three-point correlation function does not have the ''hierarchical'' form often assumed. We discuss possibilities of using the induced three-point correlation to interpret observational data. The formalism developed here can also be used to test a validity of different schemes to close the

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.

Perturbative growth of cosmological clustering .2. The two-point correlation

We use the BBGKY hierarchy equations to calculate, perturbatively, the lowest order nonlinear correction to the two-point correlation and the pair velocity for Gaussian initial conditions in a critical density matter-dominated cosmological model. We compare our results with the results obtained using the hydrodynamic equations that neglect pressure and find that the two match, indicating that there are no effects of multistreaming at this order of perturbation. We analytically study the effect of small scales on the large scales by calculating the nonlinear correction for a Dirac delta function initial two-point correlation. We find that the induced two-point correlation has a x(-6) behavior at large separations. We have considered a class of initial conditions where the initial power spectrum at small k has the form k(n) with 0 < n less than or equal to 3 and have numerically calculated the nonlinear correction to the two-point correlation, its average over a sphere and the pair velocity over a large dynamical range. We find that at small separations the effect of the nonlinear term is to enhance the clustering, whereas at intermediate scales it can act to either increase or decrease the clustering. At large scales we find a simple formula that gives a very good fit for the nonlinear correction in terms of the initial function. This formula explicitly exhibits the influence of small scales on large scales and because of this coupling the perturbative treatment breaks down at large scales much before one would expect it to if the nonlinearity were local in real space. We physically interpret this formula in terms of a simple diffusion process. We have also investigated the case n = 0, and we find that it differs from the other cases in certain respects. We investigate a recently proposed scaling property of gravitational clustering, and we find that the lowest order nonlinear terms cause deviations from the scaling relations that are strictly valid in the linear regime. The approximate validity of these relations in the nonlinear regime in l(T)-body simulations cannot be understood at this order of evolution.

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.

Low tension strings on a cosmological singularity

It has recently been argued that the singularity of the Milne orbifold can be resolved in higher spin theories. In string theory scattering amplitudes, however, the Milne singularity gives rise to ultraviolet divergences that signal uncontrolled backreaction. Since string theory in the low tension limit is expected to be a higher spin theory (although precise proposals only exist in special cases), we investigate what happens to these scattering amplitudes in the low tension limit. We point out that the known problematic ultraviolet divergences disappear in this limit. In addition we systematically identify all divergences of the simplest 2-to-2 tree-level string scattering amplitude on the Milne orbifold, and argue that the divergences that survive in the low tension limit have sensible infrared interpretations.

Exact N-wave solutions of generalized Burgers equations

Exact N-wave solutions for the generalized Burgers equation u(t) + u(n)u(x) + (j/2t + alpha) u + (beta + gamma/x) u(n+1) = delta/2u(xx),where j, alpha, beta, and gamma are nonnegative constants and n is a positive integer, are obtained. These solutions are asymptotic to the (linear) old-age solution for large time and extend the validity of the latter so as to cover the entire time regime starting where the originally sharp shock has become sufficiently thick and the viscous effects are felt in the entire N wave.

Nonlinear Enzyme Kinetics Can Lead to High Metabolic Flux Control Coefficients: Implications for the Evolution of Dominance

In a classic study, Kacser & Burns (1981, Genetics 97, 639-666) demonstrated that given certain plausible assumptions, the flux in a metabolic pathway was more or less indifferent to the activity of any of the enzymes in the pathway taken singly. It was inferred from this that the observed dominance of most wild-type alleles with respect to loss-of-function mutations did not require an adaptive, meaning selectionist, explanation. Cornish-Bowden (1987, J. theor. Biol. 125, 333-338) showed that the Kacser-Burns inference was not valid when substrate concentrations were large relative to the relevant Michaelis constants. We find that in a randomly constructed functional pathway, even when substrate levels are small, one can expect high values of control coefficients for metabolic flux in the presence of significant nonlinearities as exemplified by enzymes with Hill coefficients ranging from two to six, or by the existence of oscillatory loops. Under these conditions the flux can be quite sensitive to changes in enzyme activity as might be caused by inactivating one of the two alleles in a diploid. Therefore, the phenomenon of dominance cannot be a trivial ''default'' consequence of physiology but must be intimately linked to the manner in which metabolic networks have been moulded by natural selection.