Dynamical cluster approximation: Nonlocal dynamics of correlated electron systems

We recently introduced the dynamical cluster approximation (DCA), a technique that includes short-ranged dynamical correlations in addition to the local dynamics of the dynamical mean-field approximation while preserving causality. The technique is based on an iterative self-consistency scheme on a finite-size periodic cluster. The dynamical mean-field approximation (exact result) is obtained by taking the cluster to a single site (the thermodynamic limit). Here, we provide details of our method, explicitly show that it is causal, systematic, Phi derivable, and that it becomes conserving as the cluster size increases. We demonstrate the DCA by applying it to a quantum Monte Carlo and exact enumeration study of the two-dimensional Falicov-Kimball model. The resulting spectral functions preserve causality, and the spectra and the charge-density-wave transition temperature converge quickly and systematically to the thermodynamic limit as the cluster size increases.

Polarization relaxation, dielectric dispersion, and solvation dynamics in dense dipolar liquid

A unified treatment of polarization relaxation, dielectric dispersion and solvation dynamics in a dense, dipolar liquid is presented. It is shown that the information of solvent polarization relaxation that is obtained by macroscopic dielectric dispersion experiments is not sufficient to understand dynamics of solvation of a newly created ion or dipole. In solvation, a significant contribution comes from intermediate wave vector processes which depend critically on the short range (nearest‐neighbor) spatial and orientational order that are present in a dense, dipolar liquid. An analytic expression is obtained for the time dependent solvation energy that depends, in addition to the translational and rotational diffusion coefficients of the liquid, on the ratio of solute–solvent molecular sizes and on the microscopic structure of the polar liquid. Mean spherical approximation (MSA) theory is used to obtain numerical results for polarization relaxation, for wave vector and frequency dependent dielectric function and for time dependent solvation energy. We find that in the absence of translational contribution, the solvation of an ion is, in general, nonexponential. In this case, the short time decay is dominated by the longitudinal relaxation time but the long time decay is dominated by much slower large wave vector processes involving nearest‐neighbor molecules. The presence of a significant translational contribution drastically alters the decay behavior. Now, the long‐time behavior is given by the longitudinal relaxation time constant and the short time dynamics is controlled by the large wave vector processes. Thus, although the continuum model itself is conceptually wrong, a continuum model like result is recovered in the presence of a sizeable translational contribution. The continuum model result is also recovered in the limit of large solute to solvent size ratio. In the opposite limit of small solute size, the decay is markedly nonexponential (if the translational contribution is not very large) and a complete breakdown of the continuum model takes place. The significance of these results is discussed.

Dynamics and local structure of colossal magnetoresistance manganites

We report Extended X-ray Absorption Fine Structure and anelastic spectroscopy measurements on on hole doped manganese oxides La1-xCaxMnO3 which present the colossal magnetoresistance effect. EXAFS measurements were realized both in the absence and presence of an applied magnetic field of 1.1 Tesla, in a wide temperature range (between 330 and 77 K) and at various dopings (x = 0.25 and x = 0.33). The magnetic field orders the magnetic moments so favouring the electron mobility and the reduction of Mn-O octahedra distortions. We observe the presence of four short and two long Mn-O distances (1.93 and 2.05 Angstrom respectively) above and also below the metal-insulator phase transition. The overall distortion decreases but does not completely disappear in the metallic phase suggesting the possible coexistence of metallic and insulating regions at low temperatures. The magnetic field reduces the lattice distortions showing evidence of a microscopic counterpart of the macroscopic colossal magnetoresistance. We also present preliminary anelastic relaxation spectra in a wide temperature range from 900 K to 1 K on a sample with x = 0.40, in order to study the structural phase transitions and the lattice dynamics. A double peak has been observed at the metal-insulator transition in the imaginary part of Young's modulus. This double peak indicates that the metal-insulator transition could be a more complex phenomenon than a simple second order phase transition. In particular the peak at lower temperatures can be connected with the possible presence of inhomogeneous phase structures. Another intense dissipation peak has been observed corresponding to the structural orthorhombic-trigonal transition around 750 K.

On the generalized continuum model of dipolar solvation dynamics

The continuum model of dipolar solvation dynamics is reviewed. The effects of non-spherical molecular shapes, of non-Debye dielectric relaxation of the polar solvent and of dielectric inhomogeneity of the solvent around the solute dipole are investigated. Several new theoretical results are presented. It is found that our generalized continuum model, which takes into account the dielectric inhomogeneity of the surrounding solvent, provides a description of solvation dynamics consistent with recent experimental results.

Solvation of an ion and of a dipole in a dipolar liquid: How different are the dynamics

Theoretical expressions for the time-dependent solvation energy of an ion and of a dipole in a dense dipolar liquid are derived from microscopic considerations. We show that in contradiction to the prediction of the continuum models, the dynamics of these two species are significantly different from each other. Especially, the zero wavevector contribution, which is significant for ions, is totally absent for dipoles. Dipolar solvation may be profoundly influenced by the translational modes of the host solvent.

Dynamics of polar solvation: Route to single exponential relaxation via translational diffusion

A microscopic theoretical calculation of time-dependent solvation energy shows that the solvation of an ion or a dipole is dominated by a single relaxation time if the translational contribution to relaxation is significant.

Phase space methods and the Hamilton-Jacobi form of dynamics

A general analysis of the Hamilton-Jacobi form of dynamics motivated by phase space methods and classical transformation theory is presented. The connection between constants of motion, symmetries, and the Hamilton-Jacobi equation is described.

Monte carlo simulation for molecular gas dynamics

The dynamics of low-density flows is governed by the Boltzmann equation of the kinetic theory of gases. This is a nonlinear integro-differential equation and, in general, numerical methods must be used to obtain its solution. The present paper, after a brief review of Direct Simulation Monte Carlo (DSMC) methods due to Bird, and Belotserkovskii and Yanitskii, studies the details of theDSMC method of Deshpande for mono as well as multicomponent gases. The present method is a statistical particle-in-cell method and is based upon the Kac-Prigogine master equation which reduces to the Boltzmann equation under the hypothesis of molecular chaos. The proposed Markoff model simulating the collisions uses a Poisson distribution for the number of collisions allowed in cells into which the physical space is divided. The model is then extended to a binary mixture of gases and it is shown that it is necessary to perform the collisions in a certain sequence to obtain unbiased simulation.

Effect of √-irradiation on the structure and dynamics of domains in triglycine selenate

A comparative study of the switching properties of pure and √-irradiated TGSe crystals has been carried out to see the effect of irradiation on the structure and dynamics of domains. The switching behaviour of √-irradiated TGSe has been found to be qualitatively similar to that of unirradiated crystal and this has been interpreted in terms of structural inhibition caused by the formation of radiolysis products as well as the difference between the domain structures of the unirradiated and irradiated samples. Confirmation of this has been obtained by studying the domain patterns using the etch method.

Graphene analogue BCN: Femtosecond nonlinear optical susceptibility and hot carrier dynamics

Third-order nonlinear absorption and refraction coefficients of a few-layer boron carbon nitride (BCN) and reduced graphene oxide (RGO) suspensions have been measured at 3.2 eV in the femtosecond regime. Optical limiting behavior is exhibited by BCN as compared to saturable absorption in RGO. Nondegenerate time-resolved differential transmissions from BCN and RGO show different relaxation times. These differences in the optical nonlinearity and carrier dynamics are discussed in the light of semiconducting electronic band structure of BCN vis-a-vis the Dirac linear band structure of graphene. (C) 2010 Elsevier B.V. All rights reserved.

Investigations on the melting and bending modulus of polymer grafted bilayers using dissipative particle dynamics

Understanding the influence of polymer grafted bilayers on the physicomechanical properties of lipid membranes is important while developing liposomal based drug delivery systems. The melting characteristics and bending moduli of polymer grafted bilayers are investigated using dissipative particle dynamics simulations as a function of the amount of grafted polymer and lipid tail length. Simulations are carried out using a modified Andersen barostat, whereby the membrane is maintained in a tensionless state. For lipids made up of four to six tail beads, the transition from the low temperature L-beta phase to the L-alpha phase is lowered only above a grafting fraction of G(f)=0.12 for polymers made up of 20 beads. Below G(f)=0.12 small changes are observed only for the HT4 bilayer. The bending modulus of the bilayers is obtained as a function of G(f) from a Fourier analysis of the height fluctuations. Using the theory developed by Marsh Biochim. Biophys. Acta 1615, 33 (2003)] for polymer grafted membranes, the contributions to the bending modulus due to changes arising from the grafted polymer and bilayer thinning are partitioned. The contributions to the changes in kappa from bilayer thinning were found to lie within 11% for the lipids with four to six tail beads, increasing to 15% for the lipids containing nine tail beads. The changes in the area stretch modulus were also assessed and were found to have a small influence on the overall contribution from membrane thinning. The increase in the area per head group of the lipids was found to be consistent with the scalings predicted by self-consistent mean field results. (C) 2010 American Institute of Physics.

A molecular dynamics study of cage-to-cage migration in sodium Y zeolite: role of surface-mediated diffusion

The details of cage-to-cage migration have been obtained from an analysis of the molecular dynamics trajectory of a probe adsorbate. It is observed that particles utilize the region within a radius of 2 angstrom from the window center but with diffusion taking place predominantly at 1.6 angstrom from the window center and a potential energy of nearly -12 kJ/mol. A barrier of about 0.5 kJ/mol is observed for surface-mediated diffusion. Surprisingly, for diffusion without surface mediation for a particle going from one cage center to another, there is an attractive well near the window instead of a barrier. At low adsorbate concentrations and room temperature, the predominant mode for cage-to-cage migration is surface-mediated diffusion. The analysis suggests that particles slide along the surface of the inner walls of the alpha-cages during migration from one cage to another.

X-ray and molecular-dynamics studies on Mycobacterium leprae single-stranded DNA-binding protein and comparison with other eubacterial SSB structures

The crystal structures of two forms of Mycobacterium leprae single-stranded DNA-binding protein (SSB) have been determined at 2.05 and 2.8 A resolution. Comparison of these structures with the structures of other eubacterial SSBs indicates considerable variation in their quaternary association, although the DNA-binding domains in all of them exhibit the same OB-fold. This variation has no linear correlation with sequence variation, but could be related to variation in protein stability. Molecular-dynamics simulations have been carried out on tetrameric molecules derived from the two forms and the prototype Escherichia coli SSB and the individual subunits of both proteins. Together, the X-ray studies and molecular-dynamics simulations yield information on the relatively rigid and flexible regions of the molecule and on the effect of oligomerization on flexibility. The simulations provide insight into the changes in subunit structure on oligomerization. They also provide insight into the stability and time evolution of the hydrogen bonds/water bridges that connect the two pairs of monomers in the tetramer.

Dynamics of sheared inelastic dumbbells

We study the dynamical properties of the homogeneous shear flow of inelastic dumbbells in two dimensions as a first step towards examining the effect of shape on the properties of flowing granular materials. The dumbbells are modelled as smooth fused disks characterized by the ratio of the distance between centres (L) and the disk diameter (D), with an aspect ratio (L/D) varying between 0 and 1 in our simulations. Area fractions studied are in the range 0.1-0.7, while coefficients of normal restitution (e(n)) from 0.99 to 0.7 are considered. The simulations use a modified form of the event-driven methodology for circular disks. The average orientation is characterized by an order parameter S, which varies between 0 (for a perfectly disordered fluid) and 1 (for a fluid with the axes of all dumbbells in the same direction). We investigate power-law fits of S as a function of (L D) and (1 - e(n)(2)) There is a gradual increase in ordering as the area fraction is increased, as the aspect ratio is increased or as the coefficient of restitution is decreased. The order parameter has a maximum value of about 0.5 for the highest area fraction and lowest coefficient of restitution considered here. The mean energy of the velocity fluctuations in the flow direction is higher than that in the gradient direction and the rotational energy, though the difference decreases as the area fraction increases, due to the efficient collisional transfer of energy between the three directions. The distributions of the translational and rotational velocities are Gaussian to a very good approximation. The pressure is found to be remarkably independent of the coefficient of restitution. The pressure and dissipation rate show relatively little variation when scaled by the collision frequency for all the area fractions studied here, indicating that the collision frequency determines the momentum transport and energy dissipation, even at the lowest area fractions studied here. The mean angular velocity of the particles is equal to half the vorticity at low area fractions, but the magnitude systematically decreases to less than half the vorticity as the area fraction is increased, even though the stress tensor is symmetric.