Researcher profile

Ivan Saika-Voivod

Ivan Saika-Voivod contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

cond-mat.mtrl-sci cond-mat.stat-mech cond-mat.soft cond-mat.dis-nn

Trust snapshot

Quick read

Trust 21 - Emerging

15works

0followers

4topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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.

preprint2022arXiv

Micromagnetic simulations of clusters of nanoparticles with internal structure: Application to magnetic hyperthermia

Micromagnetic simulation results on dynamic hysteresis loops of clusters of iron oxide nanoparticles (NPs) with internal structure composed of nanorods are compared with the widely used macrospin approximation. Such calculations allowing for nanorod-composed NPs is facilitated by a previously developed coarse-graining method based on the renormalization group approach. With a focus on applications to magnetic hyperthermia, we show that magnetostatic interactions improve the heating performance of NPs in chains and triangles, and reduce heating performance in fcc arrangements. Hysteresis loops of triangular and fcc systems of complex NPs are not recovered within the macrospin approximation, especially at smaller interparticle distances. For triangular arrangements, the macrospin approximation predicts that magnetostatic interactions reduce loop area, in contrast to the complex NP case. An investigation of the local hysteresis loops of individual NPs and macrospins in clusters reveals the impact of the geometry of their neighbours on individual versus collective magnetic response, inhomogenous heating within clusters, and further differences between simulating NPs with internal structure and the use of the macrospin approximation. Capturing the internal physical and magnetic structure of NPs is thus important for some applications.

preprint2015arXiv

Crystallization of Lennard-Jones nanodroplets: from near melting to deeply supercooled

We carry out molecular dynamics (MD) and Monte Carlo (MC) simulations to characterize nucleation in liquid clusters of 600 Lennard-Jones particles over a broad range of temperatures. We use the formalism of mean first-passage times to determine the rate and find that Classical Nucleation Theory (CNT) predicts the rate quite well, even when employing simple modelling of crystallite shape, chemical potential, surface tension and particle attachment rate, down to the temperature where the droplet loses metastability and crystallization proceeds through growth-limited nucleation in an unequilibrated liquid. Below this crossover temperature, the nucleation rate is still predicted when MC simulations are used to directly calculate quantities required by CNT. Discrepancy in critical embryo sizes obtained from MD and MC arises when twinned structures with five-fold symmetry provide a competing free energy pathway out of the critical region. We find that crystallization begins with hcp-fcc stacked precritical nuclei and differentiation to various end structures occurs when these embryos become critical. We confirm that using the largest embryo in the system as a reaction coordinate is useful in determining the onset of growth-limited nucleation and show that it gives the same free energy barriers as the full cluster size distribution once the proper reference state is identified. We find that the bulk melting temperature controls the rate, even though the solid-liquid coexistence temperature for the droplet is significantly lower. The value of surface tension that renders close agreement between CNT and direct rate determination is significantly lower than what is expected for the bulk system.

preprint2015arXiv

Free energy of formation of small ice nuclei near the Widom line in simulations of supercooled water

The ST2 interaction potential has been used in a large number of simulation studies to explore the possibility of a liquid-liquid phase transition (LLPT) in supercooled water. Using umbrella sampling Monte Carlo simulations of ST2 water, we evaluate the free energy of formation of small ice nuclei in the supercooled liquid in the vicinity of the Widom line, the region above the critical temperature of the LLPT where a number of thermodynamic anomalies occur. Our results show that in this region there is a substantial free energy cost for the formation of small ice nuclei, demonstrating that the thermodynamic anomalies associated with the Widom line in ST2 water occur in a well-defined metastable liquid phase. On passing through the Widom line, we identify changes in the free energy to form small ice nuclei that illustrate how the thermodynamic anomalies associated with the LLPT may influence the ice nucleation process.

preprint2014arXiv

Inverse melting in a two-dimensional off-lattice model

We carry out computer simulations of a simple, two-dimensional off-lattice model that exhibits inverse melting. The monodisperse system comprises core-softened disks interacting through a repulsive square shoulder located inside an attractive square well. By systematically varying the potential parameters, we increase the pressure range over which the liquid freezes to a crystal upon isobaric heating. The effect is largely controlled by the extent of the shoulder. Despite occurring in two dimensions, the melting transition is first order and to a liquid, rather than to a hexatic or quasicrystal phase. We also provide comment on a commonly employed correlation function used to determine the degree of translational ordering in a system.

preprint2013arXiv

Heterogeneous Nucleation in the Low Barrier Regime

In simulations of the 2D Ising model, we examine heterogeneous nucleation induced by a small impurity consisting of a line of $l$ fixed spins. As $l$ increases, we identify a limit of stability beyond which the metastable phase is not defined. We evaluate the free energy barrier for nucleation of the stable phase, and show that contrary to expectation, the barrier does not vanish on approach to the limit of stability. We also demonstrate that our values for the height of the barrier yield predictions for the nucleation time (from transition state theory) and the size of the critical cluster (from the nucleation theorem) that are in excellent agreement with direct measurements, even near the limit of stability.

preprint2013arXiv

Understanding tetrahedral liquids through patchy colloids

We investigate the structural properties of a simple model for tetrahedral patchy colloids in which the patch width and the patch range can be tuned independently. For wide bond angles, a fully bonded network can be generated by standard Monte Carlo or molecular dynamics simulations of the model, providing a neat method for generating defect-free random tetrahedral networks. This offers the possibility of focusing on the role of the patch angular width on the structure of the fully bonded network. The analysis of the fully bonded configurations as a function of the bonding angle shows how the bonding angle controls the system compressibility, the strength of the pre-peak in the structure factor and ring size distribution. Comparison with models of liquid water and silica allows us to find the best mapping between these continuous potentials and the colloidal one. Building on previous studies focused on the connection between angular range and crystallization, the mapping makes it possible to shed new light on the glass-forming ability of network-forming tetrahedral liquids.

preprint2012arXiv

Free energy surface of ST2 water near the liquid-liquid phase transition

We carry out umbrella sampling Monte Carlo simulations to evaluate the free energy surface of the ST2 model of water as a function two order parameters, the density and a bond-orientational order parameter. We approximate the long-range electrostatic interactions of the ST2 model using the reaction-field method. We focus on state points in the vicinity of the liquid-liquid critical point proposed for this model in earlier work. At temperatures below the predicted critical temperature we find two basins in the free energy surface, both of which have liquid-like bond orientational order, but differing in density. The pressure and temperature dependence of the shape of the free energy surface is consistent with the assignment of these two basins to the distinct low density and high density liquid phases previously predicted to occur in ST2 water.

preprint2012arXiv

Phase diagram of a two-dimensional system with anomalous liquid properties

Using Monte Carlo simulation techniques, we calculate the phase diagram for a square shoulder-square well potential in two dimensions that has been previously shown to exhibit liquid anomalies consistent with a metastable liquid-liquid critical point. We consider the liquid, gas and five crystal phases, and find that all the melting lines are first order, despite a small range of metastability. One melting line exhibits a temperature maximum, as well as a pressure maximum that implies inverse melting over a small range in pressure.

preprint2011arXiv

Nucleation barriers in tetrahedral liquids spanning glassy and crystallizing regimes

Crystallization and vitrification of tetrahedral liquids are important both from a fundamental and a technological point of view. Here, we study via extensive umbrella sampling Monte Carlo computer simulations the nucleation barriers for a simple model for tetrahedral patchy particles in the regime where open tetrahedral crystal structures (namely cubic and hexagonal diamond and their stacking hybrids) are thermodynamically stable. We show that by changing the angular bond width, it is possible to move from a glass-forming model to a readily crystallizing model. From the shape of the barrier we infer the role of surface tension in the formation of the crystalline clusters. Studying the trends of the nucleation barriers with the temperature and the patch width, we are able to identify an optimal value of the patch size that leads to easy nucleation. Finally, we find that the nucleation barrier is the same, within our numerical precision, for both diamond crystals and for their stacking forms.

preprint2011arXiv

Simulation of a two-dimensional model for colloids in a uniaxial electric field

We perform Monte Carlo simulations of a simplified two-dimensional model for colloidal hard spheres in an external uniaxial AC electric field. Experimentally, the external field induces dipole moments in the colloidal particles, which in turn form chains. We therefore approximate the system as composed of well formed chains of dipolar hard spheres of a uniform length. The dipolar interaction between colloidal spheres gives rise to an effective interaction between the chains, which we treat as disks in a plane, that includes a short range attraction and long range repulsion. Hence, the system favors finite clustering over bulk phase separation and indeed we observe at low temperature and density that the system does form a cluster phase. As density increases, percolation is accompanied by a pressure anomaly. The percolated phase, despite being composed of connected, locally crystalline domains, does not bear the typical signatures of a hexatic phase. At very low densities, we find no indication of a "void phase" with a cellular structure seen recently in experiments.

preprint2010arXiv

The role of FCC tetrahedral subunits in the phase behavior of medium sized Lennard-Jones clusters

The free energy of a 600-atom Lennard-Jones cluster is calculated as a function of surface and bulk crystallinity in order to study the structural transformations that occur in the core of medium sized clusters. Within the order parameter range studied, we find the existence of two free energy minima at temperatures near freezing. One minimum, at low values of both bulk and surface order, belongs to the liquid phase. The second minimum exhibits a highly ordered core with a disordered surface and is related to structures containing a single FCC-tetrahedral subunit, with an edge length of seven atoms (l=7), located in the particle core. At lower temperatures, a third minimum appears at intermediate values of the bulk order parameter which is shown to be related to the formation of multiple l=6 tetrahedra in the core of the cluster. We also use molecular dynamics simulations to follow a series of nucleation events and find that the clusters freeze to structures containing l=5,6,7 and 8 sized tetrahedra as well as those containing no tetrahedral units. The structural correlations between bulk and surface order with the size of the tetrahedral units in the cluster core are examined. Finally, the relationships between the formation of FCC tetrahedral subunits in the core, the phase behavior of medium sized clusters and the nucleation of noncrystalline global structures such as icosahedra and decahedra are discussed.

preprint2009arXiv

Crystal Nucleation in a Supercooled Liquid with Glassy Dynamics

In simulations of supercooled, high-density liquid silica we study a range of temperature T in which we find both crystal nucleation, as well as the characteristic dynamics of a glass forming liquid, including a breakdown of the Stokes-Einstein relation. We find that the liquid cannot be observed below a homogeneous nucleation limit (HNL) at which the liquid crystallizes faster than it can equilibrate. We show that the HNL would occur at lower T, and perhaps not at all, if the Stokes-Einstein relation were obeyed, and hence that glassy dynamics plays a central role in setting a crystallization limit on the liquid state in this case. We also explore the relation of the HNL to the Kauzmann temperature, and test for spinodal-like effects near the HNL.

preprint2005arXiv

Density minimum and liquid-liquid phase transition

We present a high-resolution computer simulation study of the equation of state of ST2 water, evaluating the liquid-state properties at 2718 state points, and precisely locating the liquid-liquid critical point (LLCP) occurring in this model. We are thereby able to reveal the interconnected set of density anomalies, spinodal instabilities and response function extrema that occur in the vicinity of a LLCP for the case of a realistic, off-lattice model of a liquid with local tetrahedral order. In particular, we unambiguously identify a density minimum in the liquid state, define its relationship to other anomalies, and show that it arises due to the approach of the liquid structure to a defect-free random tetrahedral network of hydrogen bonds.

preprint2004arXiv

Phase diagram of silica from computer simulation

We evaluate the phase diagram of the ``BKS'' potential [Van Beest, Kramer and van Santen, Phys. Rev. Lett. 64, 1955 (1990)], a model of silica widely used in molecular dynamics (MD) simulations. We conduct MD simulations of the liquid, and three crystals (beta-quartz, coesite and stishovite) over wide ranges of temperature and density, and evaluate the total Gibbs free energy of each phase. The phase boundaries are determined by the intersection of these free energy surfaces. Not unexpectedly for a classical pair potential, our results reveal quantitative discrepancies between the locations of the BKS and real silica phase boundaries. At the same time, we find that the topology of the real phase diagram is reproduced, confirming that the BKS model provides a satisfactory qualitative description of a silica-like material. We also compare the phase boundaries with the locations of liquid-state thermodynamic anomalies identified in previous studies of the BKS model.

preprint2003arXiv

Free energy and configurational entropy of liquid silica: fragile-to-strong crossover and polyamorphism

Recent molecular dynamics (MD) simulations of liquid silica, using the ``BKS'' model [Van Beest, Kramer and van Santen, Phys. Rev. Lett. {\bf 64}, 1955 (1990)], have demonstrated that the liquid undergoes a dynamical crossover from super-Arrhenius, or ``fragile'' behavior, to Arrhenius, or ``strong'' behavior, as temperature $T$ is decreased. From extensive MD simulations, we show that this fragile-to-strong crossover (FSC) can be connected to changes in the properties of the potential energy landscape, or surface (PES), of the liquid. To achieve this, we use thermodynamic integration to evaluate the absolute free energy of the liquid over a wide range of density and $T$. We use this free energy data, along with the concept of ``inherent structures'' of the PES, to evaluate the absolute configurational entropy $S_c$ of the liquid. We find that the temperature dependence of the diffusion coefficient and of $S_c$ are consistent with the prediction of Adam and Gibbs, including in the region where we observe the FSC to occur. We find that the FSC is related to a change in the properties of the PES explored by the liquid, specifically an inflection in the $T$ dependence of the average inherent structure energy. In addition, we find that the high $T$ behavior of $S_c$ suggests that the liquid entropy might approach zero at finite $T$, behavior associated with the so-called Kauzmann paradox. However, we find that the change in the PES that underlies the FSC is associated with a change in the $T$ dependence of $S_c$ that elucidates how the Kauzmann paradox is avoided in this system. Finally, we also explore the relation of the observed PES changes to the recently discussed possibility that BKS silica exhibits a liquid-liquid phase transition, a behavior that has been proposed to underlie the observed polyamorphism of amorphous solid silica.