Liquid-liquid phase transition for soft-core attractive potential

Using event-driven molecular dynamics simulations, we study a three-dimensional one-component system of spherical particles interacting via a discontinuous potential combining a repulsive square soft core and an attractive square well. In the case of a narrow attractive well, it has been shown that this potential has two metastable gas-liquid critical points. Here we systematically investigate how the changes of the parameters of this potential affect the phase diagram of the system. We find a broad range of potential parameters for which the system has both a gas-liquid critical point C1 and a liquid-liquid critical point C2. For the liquid-gas critical point we find that the derivatives of the critical temperature and pressure, with respect to the parameters of the potential, have the same signs: they are positive for increasing width of the attractive well and negative for increasing width and repulsive energy of the soft core. This result resembles the behavior of the liquid-gas critical point for standard liquids. In contrast, for the liquid-liquid critical point the critical pressure decreases as the critical temperature increases. As a consequence, the liquid-liquid critical point exists at positive pressures only in a finite range of parameters. We present a modified van der Waals equation which qualitatively reproduces the behavior of both critical points within some range of parameters, and gives us insight on the mechanisms ruling the dependence of the two critical points on the potential¿s parameters. The soft-core potential studied here resembles model potentials used for colloids, proteins, and potentials that have been related to liquid metals, raising an interesting possibility that a liquid-liquid phase transition may be present in some systems where it has not yet been observed.

Dynamic properties of Lennard Jones fluids and liquid metals

The dependence of the dynamic properties of liquid metals and Lennard-Jones fluids on the characteristics of the interaction potentials is analyzed. Molecular-dynamics simulations of liquids in analogous conditions but assuming that their particles interact either through a Lennard-Jones or a liquid-metal potential were carried out. The Lennard-Jones potentials were chosen so that both the effective size of the particles and the depth of the potential well were very close to those of the liquid-metal potentials. In order to investigate the extent to which the dynamic properties of liquids depend on the short-range attractive interactions as well as on the softness of the potential cores, molecular-dynamics simulations of the same systems but assuming purely repulsive interactions with the same potential cores were also performed. The study includes both singleparticle dynamic properties, such as the velocity autocorrelation functions, and collective dynamic properties, such as the intermediate scattering funcfunctions, and collective dynamic properties, such as the intermediate scattering functions, the dynamic structure factors, the longitudinal and transverse current correlations, and the transport coefficients.

Finite-size scaling investigation of the liquid-liquid critical point in ST2 water and its stability with respect to crystallization.

The liquid-liquid critical point scenario of water hypothesizes the existence of two metastable liq- uid phases low-density liquid (LDL) and high-density liquid (HDL) deep within the supercooled region. The hypothesis originates from computer simulations of the ST2 water model, but the stabil- ity of the LDL phase with respect to the crystal is still being debated. We simulate supercooled ST2 water at constant pressure, constant temperature, and constant number of molecules N for N ≤ 729 and times up to 1 μs. We observe clear differences between the two liquids, both structural and dynamical. Using several methods, including finite-size scaling, we confirm the presence of a liquid-liquid phase transition ending in a critical point. We find that the LDL is stable with respect to the crystal in 98% of our runs (we perform 372 runs for LDL or LDL-like states), and in 100% of our runs for the two largest system sizes (N = 512 and 729, for which we perform 136 runs for LDL or LDL-like states). In all these runs, tiny crystallites grow and then melt within 1 μs. Only for N ≤ 343 we observe six events (over 236 runs for LDL or LDL-like states) of spontaneous crystal- lization after crystallites reach an estimated critical size of about 70 ± 10 molecules.

Towards and ab initio description of magnetism in ionic solids

The physical contributions to the KNiF3 magnetic exchange coupling integral have been obtained from specially designed ab initio cluster model calculations. Three important mechanisms have been identified. These are the delocalization of the magnetic orbitals into the anion p band, the variational contribution of the second-order interactions, and the many-body terms hidden in the two-body operator and the Heisenberg Hamiltonian.

Dynamics of transient pattern formation in nematic liquid crystals

We analyze the dynamics of a transient pattern formation in the Fréedericksz transition corresponding to a twist geometry. We present a calculation of the time-dependent structure factor based on a dynamical model which incorporates consistently the coupling of the director field with the velocity flow and also the effect of fluctuations. The appearance and development of a characteristic periodicity is described in terms of the time dependence of the maximum of the structure factor. We find a well-defined time for the appearance of the pattern and a subsequent stage of pattern development in which the characteristic periodicity tends to an asymptotic value.

Ab initio valence-bond cluster model for ionic solids: Alkaline-earth oxides

A linear M-O-M (M=metal, O=oxygen) cluster embedded in a Madelung field, and also including the quantum effects of the neighboring ions, is used to represent the alkaline-earth oxides. For this model an ab initio wave function is constructed as a linear combination of Slater determinants written in an atomic orbital basis set, i.e., a valence-bond wave function. Each valence-bond determinant (or group of determinants) corresponds to a resonating valence-bond structure. We have obtained ab initio valence-bond cluster-model wave functions for the electronic ground state and the excited states involved in the optical-gap transitions. Numerical results are reasonably close to the experimental values. Moreover, the model contains the ionic model as a limiting case and can be readily extended and improved.

Local character of magnetic coupling in ionic solids

Magnetic interactions in ionic solids are studied using parameter-free methods designed to provide accurate energy differences associated with quantum states defining the Heisenberg constant J. For a series of ionic solids including KNiF3, K2NiF4, KCuF3, K2CuF4, and high- Tc parent compound La2CuO4, the J experimental value is quantitatively reproduced. This result has fundamental implications because J values have been calculated from a finite cluster model whereas experiments refer to infinite solids. The present study permits us to firmly establish that in these wide-gap insulators, J is determined from strongly local electronic interactions involving two magnetic centers only thus providing an ab initio support to commonly used model Hamiltonians.

Neutral atoms in ionic lattices: Excited states of KCl:Ag(0)

The optical-absorption spectrum of a cationic Ag0 atom in a KCl crystal has been studied theoretically by means of a series of cluster models of increasing size. Excitation energies have been determined by means of a multiconfigurational self-consistent field procedure followed by a second-order perturbation correlation treatment. Moreover results obtained within the density-functional framework are also reported. The calculations confirm the assignment of bands I and IV to transitions of the Ag-5s electron into delocalized states with mainly K-4s,4p character. Bands II and III have been assigned to internal transitions on the Ag atom, which correspond to the atomic Ag-4d to Ag-5s transition. We also determine the lowest charge transfer (CT) excitation energy and confirm the assignment of band VI to such a transition. The study of the variation of the CT excitation energy with the Ag-Cl distance R gives additional support to a large displacement of the Cl ions due to the presence of the Ag0 impurity. Moreover, from the present results, it is predicted that on passing to NaCl:Ag0 the CT onset would be out of the optical range while the 5s-5p transition would undergo a redshift of 0.3 eV. These conclusions, which underline the different character of involved orbitals, are consistent with experimental findings. The existence of a CT transition in the optical range for an atom inside an ionic host is explained by a simple model, which also accounts for the differences with the more common 3d systems. The present study sheds also some light on the R dependence of the s2-sp transitions due to s2 ions like Tl+.

Rigorous characterization of oxygen vacancies in ionic oxides

Charged and neutral oxygen vacancies in the bulk and on perfect and defective surfaces of MgO are characterized as quantum-mechanical subsystems chemically bonded to the host lattice and containing most of the charge left by the removed oxygens. Attractors of the electron density appear inside the vacancy, a necessary condition for the existence of a subsystem according to the atoms in molecules theory. The analysis of the electron localization function also shows attractors at the vacancy sites, which are associated to a localization basin shared with the valence domain of the nearest oxygens. This polyatomic superanion exhibits chemical trends guided by the formal charge and the coordination of the vacancy. The topological approach is shown to be essential to understand and predict the nature and chemical reactivity of these objects. There is not a vacancy but a coreless pseudoanion that behaves as an activated host oxygen.

Dynamics of transient domain structures in nematic liquid crystals: The splay Fréedericksz transition

We discuss the dynamics of the transient pattern formation process corresponding to the splay Fréedericksz transition. The emergence and subsequent evolution of the spatial periodicity is here described in terms of the temporal dependence of the wave numbers corresponding to the maxima of the structure factor. Situations of perpendicular as well as oblique field-induced stripes relative to the initial orientation of the director are both examined with explicit indications of the time scales needed for their appearance and posterior development.

Transient patterns in nematic liquid crystals: Domain-wall dynamics

We study the dynamics of the late stages of the Fréedericksz transition in which a periodic transient pattern decays to a final homogeneous state. A stability analysis of an unstable stationary pattern is presented, and equations for the evolution of the domain walls are obtained. Using results of previous theories, we analyze the effect that the specific dynamics of the problem, incorporating hydrodynamic couplings, has on the expected logarithmic domain growth law.

Ionic-covalent transition in titanium oxides

The nature of the chemical bond in three titanium oxides of different crystal structure and different formal oxidation state has been studied by means of the ab initio cluster-model approach. The covalent and ionic contributions to the bond have been measured from different theoretical techniques. All the analysis is consistent with an increasing of covalence in the TiO, Ti2O3, and TiO2 series as expected from chemical intuition. Moreover, the use of the ab initio cluster-model approach combined with different theoretical techniques has permitted us to quantify the degree of ionic character, showing that while TiO can approximately be described as an ionic compound, TiO2 is better viewed as a rather covalent oxide.

Backflow-induced asymmetric collapse of disclination lines in liquid crystals

We present experiments where opposed pairs of planar parallel disclination lines of topological strength s=±1 move due to their mutual attraction. Our measurements show that their motion is clearly asymmetric, with +1 defects moving up to twice as fast as -1 ones. This is a clear indication of backflow, given the intrinsic isotropic elasticity of our system. A phenomenological model is able to account for the experimental observations by renormalizing the orientational diffusivity estimated from the velocity of each defect.

Chiroptical measurement of chiral aggregates at liquid-liquid interface in centrifugal liquid membrane cell by Mueller matrix and conventional circular dichroism methods

The centrifugal liquid membrane (CLM) cell has been utilized for chiroptical studies of liquid-liquid interfaces with a conventional circular dichroism (CD) spectropolarimeter. These studies required the characterization of optical properties of the rotating cylindrical CLM glass cell, which was used under the high speed rotation. In the present study, we have measured the circular and linear dichroism (CD and LD) spectra and the circular and linear birefringence (CB and LB) spectra of the CLM cell itself as well as those of porphyrine aggregates formed at the liquid-liquid interface in the CLM cell, applying Mueller matrix measurement method. From the results, it was confirmed that the CLM-CD spectra of the interfacial porphyrin aggregates observed by a conventional CD spectropolarimeter should be correct irrespective of LD and LB signals in the CLM cell.