Effective force field for liquid hydrogen fluoride from ab initio molecular dynamics simulation using the force-matching method

A recently developed force-matching method for obtaining effective force fields for condensed matter systems from ab initio molecular dynamics (MD) simulations has been applied to fit a simple nonpolarizable two-site pairwise force field for liquid hydrogen fluoride. The ab initio MD in this case was a Car-Parrinello (CP) MD simulation of 64 HF molecules at nearly ambient conditions within the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure included a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The performance of the force-match potential was examined for the gas-phase dimer and for the liquid phase at various temperatures. The model was able to reproduce correctly the bent structure and energetics of the gas-phase dimer, while the results for the structural properties, self-diffusion, vibrational spectra, density, and thermodynamic properties of liquid HF were compared to both experiment and the CP MD simulation. The force-matching model performs well in reproducing nearly all of the liquid properties as well as the aggregation behavior at different temperatures. The model is computationally cheap and compares favorably to many more computationally expensive potential energy functions for liquid HF.     

Synthesis, structure and reactivity of binuclear metal-metal bonded molybdenum(V) and tungsten(V) thioselenohalides: Molecular structure of Mo2(μ-S2)2Cl6(SeCl2)2 and W2(μ-S2)2Cl6(SeCl2)2

Thioselenohalide complexes Mo₂(μ-S₂)₂Cl₆(SeCl₂)₂ (I), Mo₂(μ-S₂)₂Br₆(SeBr₂)₂ (II), and W₂(μ-S₂)₂Br₆(SeBr₂)₂ (III) were synthesized by the reactions of corresponding metal halides or carbonyls or molybdenum metal with excesses of S₂ X ₂+Se₂ X ₂ mixtures. The complex W₂(μ-S₂)₂Cl₆(SeCl₂)₂ (IV) was obtained by an exchange reaction between (III) and excess of Se₂Cl₂. Coordination of the neutral SeX ₂ ligands to thiohalidesM ₂(μ-S₂)₂ X ₆ results in higher thermal stability, and suggests the possibility to synthesize SeX ₂ complexes of the unstable parent tungsten thiohalides. An unusual oxidative addition reaction of (I) was detected: {fx27-1} Both (I) and (IV) were characterized by X-ray crystal structure analysis. They are isostructural and form discrete molecules. Bridging S ₂ ^2? ligands are coordinated perpendicularly to the metal-metal bond;d(M?M)=2.8066 Å and 2.793 Å for I and IV, respectively. Nonequivalence of chlorine atoms which are bound to the metal atom, relate to nonequivalence of halogen atoms in the complexesM ₂(μ?S₂)₂ X ₈ ^2? . Chlorine atomstrans to SeCl₂ ligands form short bonds with the metal; the corresponding³⁵Cl NQR frequency is increased. The selenium dichloride ligand is ambidentate. The selenium atom binds as a donor to the metal and as an acceptor to two chlorine atoms which are also bound covalently to the same metal atom.     

Forward and backward pericyclic photochemical reactions have intermediates in common, yet cyclobutenes break the rules.

Photochemical pericyclic reactions are believed to proceed via a so-called pericyclic minimum on the lowest excited potential surface (S(1)), which is common to both the forward and backward reactions. Such a common intermediate has never been directly detected. The photointerconversion of 1,3-butadiene and cyclobutene is the prevailing prototype for such reactions, yet only diene ring closure proceeds with the stereospecificity that the Woodward-Hoffmann rules predict. This contrast seems to exclude a common intermediate. Using ultrafast spectroscopy, we show that the excited states of two cyclobutene/diene isomeric pairs are linked by not one, but by two common minima, p* and ct*. Starting from the diene side (cyclohepta-1,3-diene and cycloocta-1,3-diene), electrocyclic ring closure passes via the pericyclic minimum p*, whereas ct* is mainly responsible for cis-trans isomerization. Starting from the corresponding cyclobutenes (bicyclo[3.2.0]heptene-6 and bicyclo[4.2.0]octene-7), the forbidden isomer is formed from ct*. The path branches at the first (S(2)/S(1)) conical intersection towards p* and ct*. The fact that the energetically unfavorable ct* path can compete is ascribed to a dynamic effect: the momentum in C=C twist direction, acquired--such as in other olefins--in the Franck-Condon region of the cyclobutenes.     

Development of a heterogeneous chemiluminescent flow immunoassay for DDT and related compounds

A heterogeneous chemiluminescent (CL) flow immunoassay for DDT was optimized comparing different types of immunoaffinity supports: beads, nylon coils and membranes (membranes HyBondN+). In order to characterize solid immunoaffinity supports two basic immunoassay formats were performed, using (1) enzyme-labeled secondary and (2) enzyme-labeled specific monoclonal antibodies (MAbs). In both formats, hapten DDT5 conjugated to ovalbumin immobilized on solid supports according to the appropriate immobilization procedure, enzyme label (horseradish peroxidase, HRP) and luminescent detection (luminol/H₂O₂/p-iodophenol) were used. The lowest limit of detection (LOD), 1 nM p,p-DDT, was obtained with a membrane-based flow immunoassay with HRP-labeled specific antibody. Beads and packed tubing were discarded as appropriate supports because of the difficulties encountered for reproducible packing and the occurrence of light scatterring (beads), which seriously compromised the performance and reproducibility of the flow immunoassay.     

Reconstruction of the water-oxidizing complex in manganese-depleted Photosystem II using synthetic manganese complexes

The efficiency of a trinuclear and two binuclear manganese complexes in reconstituting electron transport and O(2) evolution activity in Mn-depleted Photosystem II preparations is analyzed. The trinuclear Mn-complex is more efficient than two binuclear Mn-complexes in restoring oxygen evolution, but it is less effective as an electron donor than binuclear Mn-complexes. It is inferred from our results that recovery of electron transport and O(2) evolution with polynuclear Mn-complexes is affected with different factors. Moreover, the trinuclear Mn-complex is extremely sensitive to the addition of CaCl(2). It is suggested that there is an interaction between Ca(2+) and carboxyl within the trinuclear Mn-complex during photoactivation and this interaction benefits the ligation of Mn atom to the apo-WOC and form an active WOC. Binuclear Mn(III)Mn(III) complex shows slightly higher efficiency than binuclear Mn(III)Mn(IV) complex in restoration of O(2) evolution activity. The efficiency of three Mn-complexes in the reconstitution of WOC is in an order: trinuclear Mn(3)(III)>binuclear Mn(III)Mn(III)>binuclear Mn(III)Mn(IV).     

The studies of the effect of photostimulated changes in thin carbon films

Photostimulated changes in thin carbon films deposited onto fused silica substrates by reactive magnetron sputtering and irradiated by the light of the first (1060nm) and the second (530nm) harmonics of a QW-switched Nd³⁺ AVG laser were investigated. The changes are supposed to be connected with the reordering of chemical bounds.     

Integrable structure of conformal field theory,quantum KdV theory and Thermodynamic Bethe Ansatz

We construct the quantum versions of the monodromy matrices of KdV theory. The traces of these quantum monodromy matrices, which will be called as “T-operators,” act in highest weight Virasoro modules. TheT-operators depend on the spectral parameter λ and their expansion around λ=∞ generates an infinite set of commuting Hamiltonians of the quantum KdV system. TheT-operators can be viewed as the continuous field theory versions of the commuting transfermatrices of integrable lattice theory. In particular, we show that for the values $c = 1 - 3\frac{{3(2n + 1)^2 }}{{2n + 3}}$ ,n=1,2,3 .... of the Virasoro central charge the eigenvalues of theT-operators satisfy a closed system of functional equations sufficient for determining the spectrum. For the ground-state eigenvalue these functional equations are equivalent to those of the massless Thermodynamic Bethe Ansatz for the minimal conformal field theoryM _2,2n+3; in general they provide a way to generalize the technique of the Thermodynamic Bethe Ansatz to the excited states. We discuss a generalization of our approach to the cases of massive field theories obtained by perturbing these Conformal Field Theories with the operator Φ_1,3. The relation of theseT-operators to the boundary states is also briefly described.     

The regularity of solutions of reaction-diffusion equations via Lagrangian coordinates

In this paper we study the regularity of nonnegative solutions and their interfaces for the nonlinear reaction-diffusion equation $$u_t = \left( {u^m } \right)_{xx} + f\left( u \right),\left( E \right)$$ wherem>1 andf(u) is aC ¹ function withf(0)=0 and is subject to some other technical conditions. This equation has the property of finite propagation which gives rise to interfaces separating regions whereu=0 andu>0. The analysis is carried out by means of Lagrangian coordinates, formally viewing the reaction-diffusion equation as the equation governing the evolution of the density of a certain continuum. Lagrangian coordinates have been successfully applied to study nonlinear diffusion equations posed in one space dimension. The usual formulation applies to equations which can be written in the form of a conservation law, which is not the case here because of the reaction term. In problems exhibiting interfaces such technique has the merit of rendering the interfaces straight lines, much simplifying the analysis. In this paper we present anon-standard Lagrangian formulation that works innon-conservation cases. Equation (E) is then translated into this framework and we find in a natural way the necessary estimates to prove theC ¹ regularity of moving interfaces and the regularity of the weak solution near such an interface, that allows us to establish the dynamic properties of the interface for the solutions. We end the paper by describing how the method can be applied to similar problems inseveral space dimensions with radial symmetry.     

New formulation of restricted growth processes

We present a new formulation of a class of growth models-those which evolve according to an exclusion process. This formulation is based upon a transformation of the probability distribution function which involves Grassmann variables. This method is very general and enables one to derive an exact stochastic differential equation for the model of interest. We describe this method using the traffic model as an example.     

LiNa2AlF6: a powder structure solution

Lithium sodium aluminium fluoride was obtained as a white powder by melting a stoichiometric mixture of AlF₃, NaF and LiF at 1223 K, and then cooling to 923 K and sintering at this temperature for 4 h. The ab initio crystal structure determination was carried out using X‐ray powder diffraction techniques. The monoclinic structure of LiNa₂AlF₆ can be related to cubic elpasolite. The Li and Al atoms lie on inversion centres. The main octahedral AlF₆ structural elements are not deformed, but are rotated slightly with respect to the unit‐cell axes. The Li atoms have octahedral coordinations, whereas the Na atoms have cubo‐octahedral coordinations. The Na coordination polyhedron is distorted in comparison with that of elpasolite.