BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Dmitry V. Matyushov

Dmitry V. Matyushov contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.softBiological Physicsphysics.chem-phcond-mat.stat-mechcond-mat.dis-nn

Trust snapshot

Quick read

Trust 21 - Emerging
18works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2026arXiv

Rotational Memory Function of SPC/E water

Memory effects are essential for dynamics of condensed materials and are responsible for non-exponential relaxation of correlation functions of dynamic variables through the memory function. Memory functions of dipole rotations for polar liquids have never been calculated. We present here calculations of memory functions for single-dipole rotations and for the overall dipole moment of the sample for SPC/E water. The memory functions for single-particle and collective dipole dynamics turn out to be nearly identical. This result validates theories of dielectric spectroscopy in terms of single-particle time correlation functions and the connection between the collective and single-particle relaxation times through the Kirkwood factor. The dielectric function in this formalism contains no new dynamic information that does not exist in the single-dipole correlation function. A short memory time, $\lesssim 1$ fs, justifies the use of rotational diffusion model to describe dynamics of a single molecular dipole moment in bulk water.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nonlinear dielectric relaxation of polar liquids

Molecular dynamics of two water models, SPC/E and TIP3P, at a number of temperatures is used to test the Kivelson-Madden equation connecting single-particle and collective dielectric relaxation times through the Kirkwood factor. The relation is confirmed by simulations and used to estimate the nonlinear effect of the electric field on the dielectric relaxation time. We show that the main effect of the field comes through slowing down of the single-particle rotational dynamics and the relative contribution of the field-induced alteration of the Kirkwood factor is insignificant for water. Theories of nonlinear dielectric relaxation need to mostly account for the effect of the field on rotations of a single dipole in a polar liquid.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dielectric constant of water in the interface

We define the dielectric constant (susceptibility) that should enter the Maxwell boundary value problem when applied to microscopic dielectric interfaces. The dielectric constant (susceptibility) of the interface is defined by exact linear-response equations involving correlations of statistically fluctuating interface polarization and the Coulomb interaction energy of external charges with the dielectric. The theory is applied to the interface between water and spherical solutes of altering size studied by molecular dynamics (MD) simulations. The effective dielectric constant of interfacial water is found to be significantly lower than its bulk value, and it also depends on the solute size. For TIP3P water used in MD simulations the interface dielectric constant changes from 9 to 4 when the effective solute radius is increased from ~5 to 18 A.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical and orientational structural crossovers in low-temperature glycerol

Mean square displacements of hydrogen atoms in glass-forming materials and proteins, as reported by incoherent elastic neutron scattering, show kinks in their temperature dependence. This crossover, known as the dynamical transition, connects two approximately linear regimes. It is often assigned to the dynamical freezing of subsets of molecular modes at the point of equality between their corresponding relaxation times and the instrumental observation window. The origin of the dynamical transition in glass-forming glycerol is studied here by extensive molecular dynamics simulations. We find the dynamical transition to occur for both the center of mass translations and the molecular rotations at the same temperature, insensitive to changes of the observation window. Both the translational and rotational dynamics of glycerol show a dynamic crossover from the structural to a secondary relaxation at the temperature of the dynamical transition. A significant and discontinuous increase in the orientational Kirkwood factor and in the dielectric constant is observed in the same range of temperatures. We, however, do not find any indications of a true thermodynamic transition to an ordered low-temperature phase. We therefore suggest that all observed crossovers are dynamic in character. The increase in the dielectric constant is related to the dynamic freezing of dipolar domains on the time-scale of simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mobility of nanometer-size solutes in water driven by electric field

We investigate the mobility of nanometer-size solutes in water in a uniform external electric field. General arguments are presented to show that a closed surface cutting a volume from a polar liquid will carry an effective non-zero surface charge density when preferential orientation of dipoles exists in the interface. This effective charge will experience a non-vanishing drag in an external electric field even in the absence of free charge carriers. Numerical simulations of model solutes are used to estimate the magnitude of the surface charge density. We find it to be comparable to the values typically reported from the mobility measurements. Hydrated ions can potentially carry a significant excess of the effective charge due to over-polarization of the interface. As a result, the electrokinetic charge can significantly deviate from the physical charge of free charge carriers. We propose to test the model by manipulating the polarizability of hydrated semiconductor nanoparticles with light. The inversion of the mobility direction can be achieved by photoexcitation, which increases the nanoparticle polarizability and leads to an inversion of the dipolar orientations of water molecules in the interface.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electrostatics of liquid interfaces

The standard Maxwell formulation of the problem of polarized dielectrics suffers from a number of difficulties, both conceptual and practical. These difficulties are particularly significant in the case of liquid interfaces, where the ability of the interfacial multipoles to change their orientations to minimize their free energies leads to interfacial polarization localized within a thin microscopic layer. A formalism to capture this physical reality of localized interfacial polarization is proposed and is based on the surface charge as the source of microscopic electric fields in dielectrics. The surface charge density incorporates the local structure of the interface into electrostatic calculations. The corresponding surface susceptibility and interface dielectric constant provide local closures to the electrostatic boundary value problem. A robust approach to calculate the surface susceptibility from numerical simulations is proposed. The susceptibility can alternatively be extracted from a number of solution experiments, in particular those sensitive to the overall dipole moment of a closed dielectric surface. The theory is applied to the solvent-induced spectral shift and high-frequency dielectric response of solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Depolarized light scattering and dielectric response of a peptide dissolved in water

The density and orientational relaxation of bulk water can be separately studied by depolarized light scattering (DLS) and dielectric spectroscopy (DS), respectively. Here we ask the question of what are the leading collective modes responsible for polarization anisotropy relaxation (DLS) and dipole moment relaxation (DS) of solutions involving mostly hydrophobic solute-water interfaces. We study, by atomistic molecular dynamics simulations, the dynamics and structure of hydration water interfacing N-Acetyl-leucine-methylamide (NALMA) dipeptide. The DLS response of the solution is consistent with three relaxation processes: bulk water, rotations of single solutes, and collective dipole-induced-dipole polarizability of the solutes, with the time-scale of 130-200 ps. No separate DLS response of the hydration shell has been identified by our simulations. Density fluctuations of the hydration layer, which largely contribute to the response, do not produce a dynamical process distinct from bulk water. We find that the structural perturbation of the orientational distribution of hydration waters by the dipeptide solute is quite significant and propagates 3-5 hydration layers into the bulk. This perturbation is still below that produced by hydrated globular proteins. Despite this structural perturbation, there is little change in the orientational dynamics of the hydration layers, compared to the bulk, as probed by both single-particle orientational dynamics and collective dynamics of the dipole moment of the shells. There is a clear distinction between the perturbation of the interfacial structure by the solute-solvent interaction potential and the perturbation of the interfacial dynamics by the corresponding forces.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electrophoretic mobility without charge driven by spontaneous polarization of the nanoparticle/water interface

Polarization of the interface, spontaneously occurring when water is in contact with hydrophobic solutes or air, couples with the uniform external field to produce a non-zero force acting on a suspended particle. This force exists even in the absence of a net particle charge, and its direction is affected by the first-order, dipolar and the second-order, qudrupolar orientational order parameters of the interfacial water. The quadrupolar polarization gives rise to an effectively negative charge. The corresponding surface charge density is inversely proportional to the area of the shear surface. As a result, the overall contribution from the quadrupolar polarization to the particle mobility becomes negligible compared to experimentally reported values for particles exceeding a few nanometers in size. The dipolar order of the interface dominates the zero-charge mobility of sub-micron particles. The corresponding electrokinetic charge is determined by the preferential orientation of interfacial dipoles relative to the surface normal.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Theory of Dielectric Spectroscopy of Protein Solutions

We present a theory of the dielectric response of a solution containing large solutes, of a nanometer size, in a molecular solvent. It combines the molecular dipole moment of the solute with the polarization of a large subensemble of solvent molecules at the solute-solvent interface. The goal of the theory is two-fold: (i) to formulate the problem of the dielectric response avoiding the reliance on the cavity-field concepts of dielectric theories and (ii) to separate the non-additive polarization of the interface, jointly produced by the external field of the laboratory experiment and the solute, from specific solute-solvent interactions contributing to the dielectric signal. The theory is applied to experimentally reported frequency-dependent dielectric spectra of lysozyme in solution. The analysis of the data in the broad range of frequencies up to 700 GHz shows that the cavity field susceptibility, critical for the theory formulation, is consistent with the prediction of Maxwell's electrostatics in the frequency range of 10--200 GHz, but deviates from it outside this range. In particular, it becomes much smaller then the Maxwell result and shifts to negative values at small frequencies. The latter observation implies a dia-electric response, or negative dielectrophoresis, of hydrated lysozyme. It also implies that the effective protein dipole recorded by dielectric spectroscopy is much smaller than the value calculated from protein's charge distribution. We suggest an empirical equation that describes both the increment of the static dielectric constat and the decrement of the Debye water peak with increasing protein concentration. It gives fair agreement with broad-band dispersion and loss spectra of protein solutions, but misses the delta-dispersion region.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Solvated dissipative electro-elastic network model of hydrated proteins

Elastic netwok models coarse grain proteins into a network of residue beads connected by springs. We add dissipative dynamics to this mechanical system by applying overdamped Langevin equations of motion to normal-mode vibrations of the network. In addition, the network is made heterogeneous and softened at the protein surface by accounting for hydration of the ionized residues. Solvation changes the network Hessian in two ways. Diagonal solvation terms soften the spring constants and off-diagonal dipole-dipole terms correlate displacements of the ionized residues. The model is used to formulate the response functions of the electrostatic potential and electric field appearing in theories of redox reactions and spectroscopy. We also formulate the dielectric response of the protein and find that solvation of the surface ionized residues leads to a slow relaxation peak in the dielectric loss spectrum, about two orders of magnitude slower than the main peak of protein relaxation. Finally, the solvated network is used to formulate the allosteric response of the protein to ion binding. The global thermodynamics of ion binding is not strongly affected by the network solvation, but it dramatically enhances conformational changes in response to placing a charge at the active site of the protein.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dipolar response of hydrated proteins

The paper presents an analytical theory and numerical simulations of the dipolar response of hydrated proteins. The effective dielectric constant of the solvated protein, representing the average dipole moment induced at the protein by a uniform external field, shows a remarkable variation among the proteins studied by numerical simulations. It changes from 0.5 for ubiquitin to 640 for cytochrome c. The former value implies a negative dipolar susceptibility of ubiquitin, that is a dia-electric dipolar response and negative dielectrophoresis. It means that a protein carrying an average dipole of ~240 D is expected to repel from the region of a stronger electric field. This outcome is the result of a negative cross-correlation between the protein and water dipoles, compensating for the positive variance of the protein dipole in the overall dipolar susceptibility. This phenomenon can be characterized as overscreening of protein's dipole by the hydration shell. In contrast to the neutral ubiquitin, charged proteins studied here show para-electric dipolar response and positive dielectrophoresis. The protein-water dipolar cross-correlations are long-ranged, extending approximately 2 nm from the protein surface into the bulk. The analysis of numerical simulations suggests that the polarization of the protein-water interface is strongly affected by the distribution of the protein surface charge. This component of the protein dipolar response gains in importance for high frequencies, above the protein Debye peak, when the response of the protein dipole becomes dynamically arrested. The interface response found in simulations suggests a possibility of a positive increment of the high-frequency dielectric constant of the solution compared to the dielectric constant of the solvent, in support of the observed THz absorbance of protein solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dissipative electro-elastic network model of protein electrostatics

We propose a dissipative electro-elastic network model (DENM) to describe the dynamics and statistics of electrostatic fluctuations at active sites of proteins. The model combines the harmonic network of residue beads with overdamped dynamics of the normal modes of the network characterized by two friction coefficients. The electrostatic component is introduced to the model through atomic charges of the protein force field. The overall effect of the electrostatic fluctuations of the network is recorded through the frequency-dependent response functions of the electrostatic potential and electric field at the active site. We also consider the dynamics of displacements of individual residues in the network and the dynamics of distances between pairs of residues. The model is tested against loss spectra of residue displacements and the electrostatic potential and electric field at the heme's iron from all-atom molecular dynamics simulations of three hydrated globular proteins.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Electric field inside a "Rossky cavity" in uniformly polarized water

Electric field produced inside a solute by a uniformly polarized liquid is strongly affected by dipolar polarization of the liquid at the interface. We show, by numerical simulations, that the electric "cavity" field inside a hydrated non-polar solute does not follow the predictions of standard Maxwell's electrostatics of dielectrics. Instead, the field inside the solute tends, with increasing solute size, to the limit predicted by the Lorentz virtual cavity. The standard paradigm fails because of its reliance on the surface charge density at the dielectric interface determined by the boundary conditions of the Maxwell dielectric. The interface of a polar liquid instead carries a preferential in-plane orientation of the surface dipoles thus producing virtually no surface charge. The resulting boundary conditions for electrostatic problems differ from the traditional recipes, affecting the microscopic and macroscopic fields based on them. We show that relatively small differences in cavity fields propagate into significant differences in the dielectric constant of an ideal mixture. The slope of the dielectric increment of the mixture versus the solute concentration depends strongly on which polarization scenario at the interface is realized. A much steeper slope found in the case of Lorentz polarization also implies a higher free energy penalty for polarizing such mixtures.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-Gaussian statistics of electrostatic fluctuations of hydration shells

We report the statistics of electric field fluctuations produced by SPC/E water inside a Kihara solute given as a hard-sphere core with a Lennard-Jones layer at its surface. The statistics of electric field fluctuations, obtained from numerical simulations, are studied as a function of the magnitude of a point dipole placed close to the solute-water interface. The free energy surface as a function of the electric field projected on the dipole direction shows a cross-over with the increasing dipole magnitude. While it is a single-well harmonic function at low dipole values, it becomes a double-well surface at intermediate dipole moment magnitudes, transforming to a single-well surface, with a non-zero minimum position, at still higher dipoles. A broad intermediate region where the interfacial waters fluctuate between the two minima is characterized by intense field fluctuations, with non-Gaussian statistics and the variance far exceeding the linear-response expectations. The excited state of the surface water is found to be lifted above the ground state by the energy required to break approximately two hydrogen bonds. This state is pulled down in energy by the external electric field of the solute dipole, making it readily accessible to thermal excitations. The excited state is a localized surface defect in the hydrogen-bond network creating a stress in the nearby network, but otherwise relatively localized in the region closest to the solute dipole.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Potential-tuning molecular dynamics studies of fusion, and the question of ideal glassformers: (I) The Gay-Berne model

The ability of some liquids to vitrify during supercooling is usually seen as a consequence of the rates of crystal nucleation (and/or crystal growth) becoming small- thus a matter of kinetics. However there is evidence, dating back to the empirics of coal briquetting for maximum trucking efficiency, that ellipsoids pack efficiently when disordered. Noting that key studies of non-spherical object packing have never been followed from hard ellipsoids or spherocylinders (diatomics excepted) into the world of molecules with attractive forces, we have made a molecular dynamics MD study of crystal melting and glass formation on the Gay- Berne (G-B) model of ellipsoidal objects across the aspect ratio range of the hard ellipsoid studies. Here we report that, in the aspect ratio range of maximum ellipsoid packing efficiency, various G-B crystalline states, that cannot be obtained directly from the liquid, disorder spontaneously near 0 K and transform to liquids without any detectable enthalpy of fusion. Without claiming to have proved the existence of single component examples, we use the present observations, together with our knowledge of non-ideal mixing effects, to discuss the probable existence of "ideal glassformers" - single or multicomponent liquids that vitrify before ever becoming metastable with respect to crystals. The existence of crystal-free routes to the glassy state removes any precrystalline fluctuation perspective from the "glass problem". Unexpectedly we find that liquids with aspect ratios in the "crystallophobic" range also behave in an unusual (non-hysteritic) way during temperature cycling through the glass transition. We link this to the highly volume fraction-sensitive ("fragile") behavior observed in recent hard dumbbell studies at similar length/diameter ratios.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cavity-water interface is polar

We present the results of numerical simulations of the electrostatics and dynamics of water hydration shells surrounding Kihara cavities given by a Lennard-Jones (LJ) layer at the surface of a hard-sphere cavity. The local dielectric response of the hydration layer substantially exceeds that of bulk water, with the magnitude of the dielectric constant peak in the shell increasing with the growing cavity size. The polar shell propagates into bulk water to approximately the cavity radius. The statistics of the electrostatic field produced by water inside the cavity follow linear response and approach the prediction of continuum electrostatics with increasing cavity size.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Electrostatic fluctuations promote the dynamical transition in proteins

Atomic displacements of hydrated proteins are dominated by phonon vibrations at low temperatures and by dissipative large-amplitude motions at high temperatures. A crossover between the two regimes is known as a dynamical transition. Recent experiments indicate a connection between the dynamical transition and the dielectric response of the hydrated protein. We analyze two mechanisms of the coupling between the protein atomic motions and the protein-water interface. The first mechanism considers viscoelastic changes in the global shape of the protein plasticized by its coupling to the hydration shell. The second mechanism involves modulations of the motions of partial charges inside the protein by electrostatic fluctuations. The model is used to analyze mean square displacements of iron of metmyoglobin reported by Moessbauer spectroscopy. We show that high flexibility of heme iron at physiological temperatures is dominated by electrostatic fluctuations. Two onsets, one arising from the viscoelastic response and the second from electrostatic fluctuations, are seen in the temperature dependence of the mean square displacements when the corresponding relaxation times enter the instrumental resolution window.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Gaussian excitations model for glass-former dynamics and thermodynamics

We describe a model for the thermodynamics and dynamics of glass-forming liquids in terms of excitations from an ideal glass state to a Gaussian manifold of configurationally excited states. The quantitative fit of this three parameter model to the experimental data on excess entropy and heat capacity shows that ``fragile'' behavior, indicated by a sharply rising excess heat capacity as the glass transition is approached from above, occurs in anticipation of a first-order transition -- usually hidden below the glass transition -- to a ``strong'' liquid state of low excess entropy. The dynamic model relates relaxation to a hierarchical sequence of excitation events each involving the probability of accumulating sufficient kinetic energy on a separate excitable unit. Super-Arrhenius behavior of the relaxation rates, and the known correlation of kinetic with thermodynamic fragility, both follow from the way the rugged landscape induces fluctuations in the partitioning of energy between vibrational and configurational manifolds. A relation is derived in which the configurational heat capacity, rather than the configurational entropy of the Adam Gibbs equation, controls the temperature dependence of the relaxation times, and this gives a comparable account of the experimental observations.

0 citations0 reviewsNo code linkedOpen access