RESUMEN

We extend Onsager's reciprocal relation to systems in a nonequilibrium steady state. While Onsager's reciprocal relation concerns the kinetic (Onsager) coefficient, the extended reciprocal relation concerns violation of the fluctuation response relation (FRR) for mechanical and thermal perturbations. This extended relation holds at each frequency when the extent of the FRR violation is expressed in a frequency domain. This nonintegral form distinguishes the extended relation from previous relations expressed by integration over a frequency. To obtain this relation, we consider one-particle one-dimensional systems described by an overdamped Langevin equation with a force driving the system away from equilibrium. We assume a special property of the potential in the system. From this Langevin equation, we obtain the Fokker-Planck (FP) equation describing the time evolution of the distribution function of the particle. Using the FP equation, we calculate the responses of the particle velocity and heat current by applying time-dependent perturbations of the driving force and temperature. We express the extent of the FRR violation in terms of these responses with time correlation functions and expand them in powers of the FP operator. This reciprocal relation is valid far from equilibrium. One can also confirm this reciprocal relation through experiments with systems such as colloidal suspensions because the FRR violation can be experimentally observed.

RESUMEN

A series of new Monte Carlo (MC) transition probabilities was investigated that could produce molecular trajectories statistically satisfying the diffusion equation with a position-dependent diffusion coefficient and potential energy. The MC trajectories were compared with the numerical solution of the diffusion equation by calculating the time evolution of the probability distribution and the mean first passage time, which exhibited excellent agreement. The method is powerful when investigating, for example, the long-distance and long-time global transportation of a molecule in heterogeneous systems by coarse-graining them into one-particle diffusive molecular motion with a position-dependent diffusion coefficient and free energy. The method can also be applied to many-particle dynamics.

Asunto(s)

RESUMEN

We investigate the dependence of the diffusion coefficient of a large solute particle on the solvation structure around a solute. The diffusion coefficient of a hard-sphere system is calculated by using a perturbation theory of large-particle diffusion with radial distribution functions around the solute. To obtain the radial distribution function, some integral equation theories are examined, such as the Percus-Yevick (PY), hypernetted-chain (HNC), and modified HNC theories using a bridge function proposed by Kinoshita (MHNC) closures. In one-component solvent systems, the diffusion coefficient depends on the first-minimum value of the radial distribution function. The results of the MHNC closure are in good agreement with those of calculation using the radial distribution functions of Monte Carlo simulations since the MHNC closure very closely reproduces the radial distribution function of Monte Carlo simulations. In binary-solvent mixtures, the diffusion coefficient is affected by the larger solvent density distribution in the short-range part, particularly the height and sharpness of the first peak and the depth of the first minimum. Since the HNC closure gives the first peak that is higher and sharper than that of the MHNC closure, the calculated diffusion coefficient is smaller than the MHNC closure result. In contrast, the results of the PY closure are qualitatively and quantitatively different from those of the MHNC and HNC closures.

RESUMEN

Solute-solvent reduced density profiles of hard-sphere fluids were calculated by using several integral equation theories for liquids. The traditional closures, Percus-Yevick (PY) and the hypernetted-chain (HNC) closures, as well as the theories with bridge functions, Verlet, Duh-Henderson, and Kinoshita (named MHNC), were used for the calculation. In this paper, a one-solute hard-sphere was immersed in a one-component hard-sphere solvent and various size ratios were examined. The profiles between the solute and solvent particles were compared with those calculated by Monte Carlo simulations. The profiles given by the integral equations with the bridge functions were much more accurate than those calculated by conventional integral equation theories, such as the Ornstein-Zernike (OZ) equation with the PY closure. The accuracy of the MHNC-OZ theory was maintained even when the particle size ratio of solute to solvent was 50. For example, the contact values were 5.7 (Monte Carlo), 5.6 (MHNC), 7.8 (HNC), and 4.5 (PY), and the first minimum values were 0.48 (Monte Carlo), 0.46 (MHNC), 0.54 (HNC), and 0.40 (PY) when the packing fraction of the hard-sphere solvent was 0.38 and the size ratio was 50. The asymptotic decay and the oscillation period for MHNC-OZ were also very accurate, although those given by the HNC-OZ theory were somewhat faster than those obtained by Monte Carlo simulations.

RESUMEN

We have studied the diffusion of a large hard-sphere solute immersed in binary hard-sphere mixtures. We reveal how the boundary condition at the solute surface is affected by the solvent density around the solute. Solving equations for a binary compressible mixture by perturbation expansions, we obtain the boundary condition depending on the size ratio of binary solvent spheres. When the size ratio is 1:2, the boundary condition lies close to the slip boundary condition. By contrast, when the size ratio becomes large, the boundary condition approaches the stick boundary condition with the addition of larger solvent spheres. We find that the transition to the stick boundary condition is caused by the increase in the solvent density around the solute due to an entropic effect.

RESUMEN

Insertion of a solute into a vessel comprising biopolymers is a fundamental function in a biological system. The entropy originating from the translational displacement of solvent particles plays an essential role in the insertion. Here we study the dynamics of entropic insertion of a large spherical solute into a cylindrical vessel. The solute and the vessel are immersed in small spheres forming the solvent. We develop a theoretical method formulated using the Fokker-Planck equation. The spatial distribution of solute-vessel entropic potential, which is calculated by the three-dimensional integral equation theory combined with rigid-body models, serves as input data. The key quantity analyzed is the density of the probability of finding the solute at any position at any time. It is found that the solute is inserted along the central axis of the vessel cavity and trapped at a position where the entropic potential takes a local minimum value. The solute keeps being trapped without touching the vessel inner surface. In a significantly long time τ, the solute transfers to the position in contact with the vessel bottom possessing the global potential minimum along the central axis. As the solute size increases, τ becomes remarkably longer. We also discuss the relevance of our result to the functional expression of a chaperonin/cochaperonin in the assistance of protein folding.

Asunto(s)

RESUMEN

We study phase stability of a system with double-minimum interaction potential in a wide range of parameters by a thermodynamic perturbation theory. The present double-minimum potential is the Lennard-Jones-Gauss potential, which has a Gaussian pocket as well as a standard Lennard-Jones minimum. As a function of the depth and position of the Gaussian pocket in the potential, we determine the coexistence pressure of crystals (fcc and bcc). We show that the fcc crystallizes even at zero pressure when the position of the Gaussian pocket is coincident with the first or third nearest neighbor site of the fcc crystal. The bcc crystal is more stable than the fcc crystal when the position of the Gaussian pocket is coincident with the second nearest neighbor sites of the bcc crystal. The stable crystal structure is determined by the position of the Gaussian pocket. These results show that we can control the stability of the solid phase by tuning the potential function.

RESUMEN

Exact equations of motion for the microscopically defined collective density ρ(x,t) and the momentum density g(x,t) of a fluid have been obtained in the past starting from the corresponding Langevin equations representing the dynamics of the fluid particles. In the present work we average these exact equations of microscopic dynamics over the local equilibrium distribution to obtain stochastic partial differential equations for the coarse-grained densities with smooth spatial and temporal dependence. In particular, we consider Dean's exact balance equation for the microscopic density of a system of interacting Brownian particles to obtain the basic equation of the dynamic density functional theory with noise. Our analysis demonstrates that on thermal averaging the dependence of the exact equations on the bare interaction potential is converted to dependence on the corresponding thermodynamic direct correlation functions in the coarse-grained equations.

RESUMEN

We have calculated the dielectric relaxation of water around an ion using molecular dynamics simulations. The collective motion of water near the ion showed fast relaxation, whereas the reorientational motion of individual water molecules does not have the fast component. The ratio of the relaxation time for the fast component and the bulk water was consistent with the experimental results, known as hyper-mobile water, for alkali halide aqueous solution.

RESUMEN

Exploiting the thermodynamic potential functional provided by density functional theory, we determine analytically the free-energy landscape (FEL) in a hard-sphere fluid. The FEL is represented in the three-dimensional coordinate space of the tagged particle. We also analyze the distribution of the free-energy barrier between adjacent basins and show that the most provable value and the average of the free-energy barrier are increasing functions of the density. Since the size of the cooperatively rearranging region (CRR) is also increased as the density is raised [Yoshidome, Phys. Rev. E 76, 021506 (2007)], the present result is consistent with the Adam-Gibbs theory in which the increase of the activation energy is due to the increase of the size of the CRR.

RESUMEN

Validity of the centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD) in quantum liquids is studied on an assumption that momenta of liquid particles relax fast. The projection operator method allows one to derive the generalized Langevin equation including a memory effect for the full-quantum canonical (Kubo-transformed) correlation function. Similar equations for the CMD and RPMD correlation functions can be derived too. The comparison of these equations leads to conditions under which the RPMD and CMD correlation functions agree approximately with the full-quantum canonical correlation function. The condition for the RPMD is that the memory effects of the full-quantum and RPMD equations vanish quickly with the same time constants. The CMD correlation function requires additional conditions concerning static correlation.

RESUMEN

Exploiting the density functional theory, we calculate the free energy landscape (FEL) of the hard sphere glass in three dimensions. From the FEL, we estimate the number of the particles in the cooperatively rearranging region (CRR). We find that the density dependence of the number of the particles in the CRR is expressed as a power law function of the density. Analyzing the relaxation process in the CRR, we also find that the string motion is the elementary process for the structural relaxation, which leads to the natural definition of the simultaneously rearranging region as the particles displaced in the string motion.

RESUMEN

The authors applied the time dependent density functional method (TDDFM) and a linear model to solvation dynamics in simple binary solvents. Changing the solute-solvent interactions at t=0, the authors calculated the time evolution of density fields for solvent particles after the change (t>0) by the TDDFM and linear model. First, the authors changed the interaction of only one component of solvents. In this case, the TDDFM showed that the solvation time decreased monotonically with a mole fraction of the solvent strongly interacting with the solute. The monotonical decreases agreed with experimental results, while the linear model did not reproduce these results. The authors also calculated the solvation time by changing the interaction of both components. The calculation showed that the mole fraction dependence had the peak. The TDDFM presented a much higher peak than the linear model. The difference between the TDDFM and the linear model was caused by a nonlinear effect on an exchange process of solvent particles.

RESUMEN

Time-dependent density functional methods (TDDFM) are studied from the microscopic viewpoint using projection operator methods in classical liquids. A density field is defined without averaging, so that a time evolution equation of the density field is derived with a random force. The derived equation includes a free energy functional, which is different from that defined in the TDDFM. The projection operator method provides the exact expression of the free energy functional. Another definition of the density field by an average leads to the equation of the TDDFM. In addition, an equation describing fluctuations is also derived.

RESUMEN

In order to understand the behavior of thermodynamic quantities near the glass transition temperature, we put the energy landscape picture and the particle's jump motion together and calculate the specific heat of a nonequilibrium system. Taking the finite observation time into account, we study the observation time dependence of the specific heat. We assume the Einstein oscillators for the dynamics of each basin in the landscape structure of phase space and calculate the specific heat of a system with 20 basins. For a given observation time, a transition from annealed to quenched system occurs at the temperature when the time scale of jumps exceeds the observation time. The transition occurs at lower temperature and becomes sharper for longer observation time.

RESUMEN

Uncovering why spatial mosaics of mimetic morphs are maintained in a Müllerian mimicry system has been a challenging issue in evolutionary biology. In this article, we analyze the reaction diffusion system that describes two-species Müllerian mimicry in one- and two-dimensional habitats. Due to positive frequency-dependent selection, a local population first approaches the state where one of the comimicking patterns predominates, which is followed by slow movement of boundaries where different patterns meet. We then analyze the interfacial dynamics of the boundaries to find whether a stable cline is maintained and to obtain the wave speed if the cline is unstable. The results are: (1) In a spatially uniform habitat the morph with greater base fitness spreads both in one and two species system. (2) The strength of cross-species interaction determines whether the mimetic morph clines of model and mimic species coalesce into the same geographical region or pass through each other. The joint wave speed of clines decreases by increasing the number of comimicking species in the mimicry ring. (3) In spatial heterogeneous habitats, stable clines can be maintained due to the balance between the base fitness gradient and the biased gene flow by negative curvature of boundary. This allows the persistence of a spatial mosaic even if one of the morphs is in every place advantageous over the other. A balanced cline is also maintained if there is a gradient in the population density. (4) A new advantageous morph occurring at a local region is doomed to go to extinction in a finite time if the "radius" of initial distribution is below a threshold. Possible applications to the heliconiine butterfly mimicry ring, heterozygous disadvantage systems of chromosomal rearrangement and hybrid zone, the third phase of Wright's Shifting Balance theory, and cytoplasmic incompatibility are discussed.

Asunto(s)