Vili Heinonen | Massachusetts Institute of Technology (MIT) (original) (raw)
Papers by Vili Heinonen
The Journal of Chemical Physics, 2013
We determine the orientation-resolved interfacial free energy between a body-centered-cubic (bcc)... more We determine the orientation-resolved interfacial free energy between a body-centered-cubic (bcc) crystal and the coexisting fluid for a many-particle system interacting via a Yukawa pair potential. For two different screening strengths, we compare results from molecular dynamics computer simulations, density functional theory, and a phase-field-crystal approach. Simulations predict an almost orientationally isotropic interfacial free energy of 0.12k B T/a 2 (with k B T denoting the thermal energy and a the mean interparticle spacing), which is independent of the screening strength. This value is in reasonable agreement with our Ramakrishnan-Yussouff density functional calculations, while a high-order fitted phase-field-crystal approach gives about 2−3 times higher interfacial free energies for the Yukawa system. Both field theory approaches also give a considerable anisotropy of the interfacial free energy. Our result implies that, in the Yukawa system, bcc crystal-fluid free energies are a factor of about 3 smaller than face-centered-cubic crystal-fluid free energies.
Physical Review Letters, 2010
Even below its normal melting temperature, ice melts when subjected to high pressure and refreeze... more Even below its normal melting temperature, ice melts when subjected to high pressure and refreezes once the pressure is lifted. A classic demonstration of this regelation phenomenon is the passing of a thin wire through a block of ice when sufficient force is exerted. Here we present a molecular-dynamics study of a nanowire cutting through ice to unravel the molecular level mechanisms responsible for regelation. In particular, we show that the transition from a stationary to a moving wire due to increased driving force changes from symmetric and continuous to asymmetric and discontinuous as a hydrophilic wire is replaced by a hydrophobic one. This is explained at the molecular level in terms of the wetting properties of the wire.
Phase Field Crystal (PFC) models and their amplitude expansions are a novel attempt to bridge the... more Phase Field Crystal (PFC) models and their amplitude expansions are a novel attempt to bridge the gap between atomistic and continuum models in materials modeling. The studied quantity is the atomic density field that varies in time and space. Not only do these new models allow longer length scales but also longer time scales: many interesting phenomena happen over diffusive time scales that are beyond the reach of classical molecular dynamics or Monte Carlo methods. As with Dynamic Density Field theory a local equilibrium in the system is assumed leading into diffusive dynamics. This implies that the time scale under study is lot slower than time scale of elastic excitations in the system. In other words, we are assuming that phonon modes die out instantaneously. However, it turns out that the system exhibits elastic excitations that have to be relaxed separately. We propose a physical constraint to the time evolution of the density field, which ensures elastic equilibrium at all t...
Physical Review Letters, 2010
Even below its normal melting temperature, ice melts when subjected to high pressure and refreeze... more Even below its normal melting temperature, ice melts when subjected to high pressure and refreezes once the pressure is lifted. A classic demonstration of this regelation phenomenon is the passing of a thin wire through a block of ice when sufficient force is exerted. Here we present a molecular-dynamics study of a nanowire cutting through ice to unravel the molecular level mechanisms responsible for regelation. In particular, we show that the transition from a stationary to a moving wire due to increased driving force changes from symmetric and continuous to asymmetric and discontinuous as a hydrophilic wire is replaced by a hydrophobic one. This is explained at the molecular level in terms of the wetting properties of the wire.
Physical Review E, 2014
Phase field crystal (PFC) models constitute a field theoretical approach to solidification, melti... more Phase field crystal (PFC) models constitute a field theoretical approach to solidification, melting and related phenomena at atomic length and diffusive time scales. One of the advantages of these models is that they naturally contain elastic excitations associated with strain in crystalline bodies. However, instabilities that are diffusively driven towards equilibrium are often orders of magnitude slower than the dynamics of the elastic excitations, and are thus not included in the standard PFC model dynamics. We derive a method to isolate the time evolution of the elastic excitations from the diffusive dynamics in the PFC approach and set up a two-stage process, in which elastic excitations are equilibrated separately. This ensures mechanical equilibrium at all times. We show concrete examples demonstrating the necessity of the separation of the elastic and diffusive time scales. In the small deformation limit this approach is shown to agree with the theory of linear elasticity. arXiv:1311.7336v1 [cond-mat.mtrl-sci]
The Journal of Chemical Physics, 2013
We determine the orientation-resolved interfacial free energy between a body-centered-cubic (bcc)... more We determine the orientation-resolved interfacial free energy between a body-centered-cubic (bcc) crystal and the coexisting fluid for a many-particle system interacting via a Yukawa pair potential. For two different screening strengths, we compare results from molecular dynamics computer simulations, density functional theory, and a phase-field-crystal approach. Simulations predict an almost orientationally isotropic interfacial free energy of 0.12k B T/a 2 (with k B T denoting the thermal energy and a the mean interparticle spacing), which is independent of the screening strength. This value is in reasonable agreement with our Ramakrishnan-Yussouff density functional calculations, while a high-order fitted phase-field-crystal approach gives about 2−3 times higher interfacial free energies for the Yukawa system. Both field theory approaches also give a considerable anisotropy of the interfacial free energy. Our result implies that, in the Yukawa system, bcc crystal-fluid free energies are a factor of about 3 smaller than face-centered-cubic crystal-fluid free energies.
Computer Physics Communications, 2012
The three-dimensional Mercedes-Benz model was recently introduced to account for the structural a... more The three-dimensional Mercedes-Benz model was recently introduced to account for the structural and thermodynamic properties of water. It treats water molecules as point-like particles with four dangling bonds in tetrahedral coordination, representing H-bonds of water. Its conceptual simplicity renders the model attractive in studies where complex behaviors emerge from H-bond interactions in water, e.g., the hydrophobic effect. A molecular dynamics (MD) implementation of the model is non-trivial and we outline here the mathematical framework of its force-field. Useful routines written in modern Fortran are also provided. This open source code is free and can easily be modified to account for different physical context. The provided code allows both serial and MPI-parallelized execution.
Journal of Chemical Physics 138, 044705 (2013)
We determine the orientation-resolved interfacial free energy between a body-centered-cubic (bcc)... more We determine the orientation-resolved interfacial free energy between a body-centered-cubic (bcc) crystal and the coexisting fluid for a many-particle system interacting via a Yukawa pair potential. For two different screening strengths, we compare results from molecular dynamics computer simulations, density functional theory, and a phase-field-crystal approach. Simulations predict an almost orientationally isotropic interfacial free energy of 0.12k B T/a 2 (with k B T denoting the thermal energy and a the mean interparticle spacing), which is independent of the screening strength. This value is in reasonable agreement with our Ramakrishnan-Yussouff density functional calculations, while a high-order fitted phase-field-crystal approach gives about 2−3 times higher interfacial free energies for the Yukawa system. Both field theory approaches also give a considerable anisotropy of the interfacial free energy. Our result implies that, in the Yukawa system, bcc crystal-fluid free energies are a factor of about 3 smaller than face-centered-cubic crystal-fluid free energies."
The Journal of Chemical Physics, 2013
We determine the orientation-resolved interfacial free energy between a body-centered-cubic (bcc)... more We determine the orientation-resolved interfacial free energy between a body-centered-cubic (bcc) crystal and the coexisting fluid for a many-particle system interacting via a Yukawa pair potential. For two different screening strengths, we compare results from molecular dynamics computer simulations, density functional theory, and a phase-field-crystal approach. Simulations predict an almost orientationally isotropic interfacial free energy of 0.12k B T/a 2 (with k B T denoting the thermal energy and a the mean interparticle spacing), which is independent of the screening strength. This value is in reasonable agreement with our Ramakrishnan-Yussouff density functional calculations, while a high-order fitted phase-field-crystal approach gives about 2−3 times higher interfacial free energies for the Yukawa system. Both field theory approaches also give a considerable anisotropy of the interfacial free energy. Our result implies that, in the Yukawa system, bcc crystal-fluid free energies are a factor of about 3 smaller than face-centered-cubic crystal-fluid free energies.
Physical Review Letters, 2010
Even below its normal melting temperature, ice melts when subjected to high pressure and refreeze... more Even below its normal melting temperature, ice melts when subjected to high pressure and refreezes once the pressure is lifted. A classic demonstration of this regelation phenomenon is the passing of a thin wire through a block of ice when sufficient force is exerted. Here we present a molecular-dynamics study of a nanowire cutting through ice to unravel the molecular level mechanisms responsible for regelation. In particular, we show that the transition from a stationary to a moving wire due to increased driving force changes from symmetric and continuous to asymmetric and discontinuous as a hydrophilic wire is replaced by a hydrophobic one. This is explained at the molecular level in terms of the wetting properties of the wire.
Phase Field Crystal (PFC) models and their amplitude expansions are a novel attempt to bridge the... more Phase Field Crystal (PFC) models and their amplitude expansions are a novel attempt to bridge the gap between atomistic and continuum models in materials modeling. The studied quantity is the atomic density field that varies in time and space. Not only do these new models allow longer length scales but also longer time scales: many interesting phenomena happen over diffusive time scales that are beyond the reach of classical molecular dynamics or Monte Carlo methods. As with Dynamic Density Field theory a local equilibrium in the system is assumed leading into diffusive dynamics. This implies that the time scale under study is lot slower than time scale of elastic excitations in the system. In other words, we are assuming that phonon modes die out instantaneously. However, it turns out that the system exhibits elastic excitations that have to be relaxed separately. We propose a physical constraint to the time evolution of the density field, which ensures elastic equilibrium at all t...
Physical Review Letters, 2010
Even below its normal melting temperature, ice melts when subjected to high pressure and refreeze... more Even below its normal melting temperature, ice melts when subjected to high pressure and refreezes once the pressure is lifted. A classic demonstration of this regelation phenomenon is the passing of a thin wire through a block of ice when sufficient force is exerted. Here we present a molecular-dynamics study of a nanowire cutting through ice to unravel the molecular level mechanisms responsible for regelation. In particular, we show that the transition from a stationary to a moving wire due to increased driving force changes from symmetric and continuous to asymmetric and discontinuous as a hydrophilic wire is replaced by a hydrophobic one. This is explained at the molecular level in terms of the wetting properties of the wire.
Physical Review E, 2014
Phase field crystal (PFC) models constitute a field theoretical approach to solidification, melti... more Phase field crystal (PFC) models constitute a field theoretical approach to solidification, melting and related phenomena at atomic length and diffusive time scales. One of the advantages of these models is that they naturally contain elastic excitations associated with strain in crystalline bodies. However, instabilities that are diffusively driven towards equilibrium are often orders of magnitude slower than the dynamics of the elastic excitations, and are thus not included in the standard PFC model dynamics. We derive a method to isolate the time evolution of the elastic excitations from the diffusive dynamics in the PFC approach and set up a two-stage process, in which elastic excitations are equilibrated separately. This ensures mechanical equilibrium at all times. We show concrete examples demonstrating the necessity of the separation of the elastic and diffusive time scales. In the small deformation limit this approach is shown to agree with the theory of linear elasticity. arXiv:1311.7336v1 [cond-mat.mtrl-sci]
The Journal of Chemical Physics, 2013
We determine the orientation-resolved interfacial free energy between a body-centered-cubic (bcc)... more We determine the orientation-resolved interfacial free energy between a body-centered-cubic (bcc) crystal and the coexisting fluid for a many-particle system interacting via a Yukawa pair potential. For two different screening strengths, we compare results from molecular dynamics computer simulations, density functional theory, and a phase-field-crystal approach. Simulations predict an almost orientationally isotropic interfacial free energy of 0.12k B T/a 2 (with k B T denoting the thermal energy and a the mean interparticle spacing), which is independent of the screening strength. This value is in reasonable agreement with our Ramakrishnan-Yussouff density functional calculations, while a high-order fitted phase-field-crystal approach gives about 2−3 times higher interfacial free energies for the Yukawa system. Both field theory approaches also give a considerable anisotropy of the interfacial free energy. Our result implies that, in the Yukawa system, bcc crystal-fluid free energies are a factor of about 3 smaller than face-centered-cubic crystal-fluid free energies.
Computer Physics Communications, 2012
The three-dimensional Mercedes-Benz model was recently introduced to account for the structural a... more The three-dimensional Mercedes-Benz model was recently introduced to account for the structural and thermodynamic properties of water. It treats water molecules as point-like particles with four dangling bonds in tetrahedral coordination, representing H-bonds of water. Its conceptual simplicity renders the model attractive in studies where complex behaviors emerge from H-bond interactions in water, e.g., the hydrophobic effect. A molecular dynamics (MD) implementation of the model is non-trivial and we outline here the mathematical framework of its force-field. Useful routines written in modern Fortran are also provided. This open source code is free and can easily be modified to account for different physical context. The provided code allows both serial and MPI-parallelized execution.
Journal of Chemical Physics 138, 044705 (2013)
We determine the orientation-resolved interfacial free energy between a body-centered-cubic (bcc)... more We determine the orientation-resolved interfacial free energy between a body-centered-cubic (bcc) crystal and the coexisting fluid for a many-particle system interacting via a Yukawa pair potential. For two different screening strengths, we compare results from molecular dynamics computer simulations, density functional theory, and a phase-field-crystal approach. Simulations predict an almost orientationally isotropic interfacial free energy of 0.12k B T/a 2 (with k B T denoting the thermal energy and a the mean interparticle spacing), which is independent of the screening strength. This value is in reasonable agreement with our Ramakrishnan-Yussouff density functional calculations, while a high-order fitted phase-field-crystal approach gives about 2−3 times higher interfacial free energies for the Yukawa system. Both field theory approaches also give a considerable anisotropy of the interfacial free energy. Our result implies that, in the Yukawa system, bcc crystal-fluid free energies are a factor of about 3 smaller than face-centered-cubic crystal-fluid free energies."