Researcher profile

Amir Haji-Akbari

Amir Haji-Akbari contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

9 published item(s)

preprint2021arXiv

How to Quantify and Avoid Finite Size Effects in Computational Studies of Crystal Nucleation: The Case of Heterogeneous Ice Nucleation

Computational studies of crystal nucleation can be impacted by finite size effects, primarily due to unphysical interactions between crystalline nuclei and their periodic images. It is, however, not always feasible to systematically investigate the sensitivity of nucleation kinetics and mechanism to system size due to large computational costs of nucleation studies. Here, we use jumpy forward flux sampling to accurately compute the rates of heterogeneous ice nucleation in the vicinity of square-shaped model structureless ice nucleating particles (INPs) of different sizes, and identify three distinct regimes for the dependence of rate on the INP dimension, $L$. For small INPs, the rate is a strong function of $L$ due to artificial spanning of critical nuclei across the periodic boundary. Intermediate-sized INPs, however, give rise to the emergence of non-spanning 'proximal` nuclei that are close enough to their periodic images to fully structure the intermediary liquid. While such proximity can facilitate nucleation, its effect is offset by the higher density of the intermediary liquid, leading to artificially small nucleation rates overall. The critical nuclei formed at large INPs are neither spanning nor proximal. Yet, the rate is a weak function of $L$, with its logarithm scaling linearly with $1/L$. The key heuristic emerging from these observations is that finite size effects will be minimal if critical nuclei are neither spanning nor proximal, and if the intermediary liquid has a region that is structurally indistinguishable from the supercooled liquid under the same conditions.

preprint2021arXiv

How to quantify and avoid finite size effects in computational studies of crystal nucleation: The case of homogeneous crystal nucleation

Finite size artifacts arise in molecular simulations of nucleation when critical nuclei are too close to their periodic images. A rigorous determination of what constitutes too close is, however, a major challenge. Recently, we devised rigorous heuristics for detecting such artifacts based on our investigation of how system size impacts the rate of heterogeneous ice nucleation (Hussain, Haji-Akbari, \emph{J. Chem. Phys.} \textbf{154}, 014108, \textbf{2021}). We identified the prevalence of critical nuclei spanning across the periodic boundary, and the thermodynamic and structural properties of the liquid occupying the inter-image region as indicators of finite size artifacts. Here, we further probe the performance of such heuristics by examining the dependence of homogeneous crystal nucleation rates in the Lennard-Jones liquid on system size. The rates depend non-monotonically on system size and vary by almost six orders of magnitude for the range of system sizes considered here. We confirm that the prevalence of spanning critical nuclei is the primary indicator of finite size artifacts and almost fully explains the observed variations in rate. Proximity, or structuring of the inter-image liquid, however, is not as strong of an indicator due to the fragmented nature of crystalline nuclei. As a result, the dependence of rate on system size is subtle for the systems with a minuscule fraction of spanning critical nuclei. These observations indicate that our heuristics are universally applicable to different modes of nucleation (homogeneous and heterogeneous) in different systems even if they might be overly stringent for homogeneous nucleation,~e.g.,~in the LJ system.

preprint2021arXiv

Role of Nanoscale Interfacial Proximity in Contact Freezing in Water

Contact freezing is a mode of atmospheric ice nucleation in which a collision between a dry ice nucleating particle (INP) and a water droplet results in considerably faster heterogeneous nucleation. The molecular mechanism of such enhancement is, however, still a mystery. While earlier studies had attributed it to collision-induced transient perturbations, recent experiments point to the pivotal role of nanoscale proximity of the INP and the free interface. By simulating heterogeneous nucleation of ice within INP-supported nanofilms of two model water-like tetrahedral liquids, we demonstrate that such nanoscale proximity is sufficient for inducing rate increases commensurate with those observed in contact freezing experiments, but only if the free interface has a tendency to enhance homogeneous nucleation. Water is suspected of possessing this latter property, known as surface freezing propensity. Our findings therefore establish a connection between surface freezing propensity and kinetic enhancement during contact nucleation. We also observe that faster nucleation proceeds through a mechanism markedly distinct from classical heterogeneous nucleation, involving the formation of hourglass-shaped crystalline nuclei that conceive at either interface, and that have a lower free energy of formation due to the nanoscale proximity of the interfaces and the modulation of the free interfacial structure by the INP. In addition to providing valuable insights into the physics of contact nucleation, our findings can assist in improving the accuracy of heterogeneous nucleation rate measurements in experiments, and in advancing our understanding of ice nucleation on nonuniform surfaces such as organic, polymeric and biological materials.

preprint2020arXiv

Induced Charge Anisotropy: a Hidden Variable Affecting Ion Transport through Membranes

The ability of semipermeable membranes to selectively impede the transport of undesirable solutes is key to many applications. Yet, obtaining a systematic understanding of how membrane structure affects selectivity remains elusive due to the insufficient spatiotemporal resolution of existing experimental techniques, and the inaccessibility of relevant solute transport timescales to conventional molecular simulations. Here, we utilize jumpy forward-flux sampling to probe the transport of sodium and chloride ions through a graphitic membrane with sub-nm pores. We find chlorides to traverse the pore at rates over two orders of magnitude faster than sodiums. We also identify two major impediments to the transport of both ion types. In addition to the partial dehydration of the leading ion, its traversal induces charge anisotropy at its rear, which exerts a net restraining force on the ion. Charge anisotropy is therefore a crucial hidden variable controlling the kinetics of ion transport through nanopores.

preprint2019arXiv

Studying Rare Events using Forward-Flux Sampling: Recent Breakthroughs and Future Outlook

Rare events are processes that occur upon the emergence of unlikely fluctuations. Unlike what their name suggests, rare events are fairly ubiquitous in nature, as the occurrence of many structural transformations in biology and material sciences is predicated upon crossing large free energy barriers. Probing the kinetics and uncovering the molecular mechanisms of possible barrier crossings in a system is critical to predicting and controlling its structural and functional properties. Due to their activated nature, however, rare events are exceptionally difficult to study using conventional experimental and computational techniques. In recent decades, a wide variety of specialized computational techniques-- known as advanced sampling techniques-- have been developed to systematically capture improbable fluctuations relevant to rare events. In this perspective, we focus on a technique called forward flux sampling (Allen~\emph{et al.,}~\emph{J. Chem. Phys.}, {\bf 124}: 024102, 2006), and overview its recent methodological variants and extensions. We also provide a detailed overview of its application to study a wide variety of rare events, and map out potential avenues for further explorations.

preprint2013arXiv

Packing and Self-assembly of Truncated Triangular Bipyramids

Motivated by breakthroughs in the synthesis of faceted nano- and colloidal particles, as well as theoretical and computational studies of their packings, we investigate a family of truncated triangular bipyramids. We report dense periodic packings with small unit cells that were obtained via numerical and analytical optimization. The maximal packing fraction $ϕ_{\max}$ changes continuously with the truncation parameter $t$. Eight distinct packings are identified based on discontinuities in the first and second derivatives of $ϕ_{\max}(t)$. These packings differ in the number of particles in the fundamental domain (unit cell) and the type of contacts between the particles. In particular, we report two packings with four particles in the unit cell for which both $ϕ_{\max}(t)$ and $ϕ'_{\max}(t)$ are continuous and the discontinuity occurs in the second derivative only. In the self-assembly simulations that we perform for larger boxes with 2048 particles, only one out of eight packings is found to assemble. In addition, the degenerate quasicrystal reported previously for triangular bipyramids without truncation [{http://prl.aps.org/abstract/PRL/v107/i21/e215702}{Haji-Akbari \emph{et al.}, \emph{Phys. Rev. Lett.} \textbf{107}: 215702 (2011)}] assembles for truncations as high as 0.45. The self-assembly propensities for the structures formed in the thermodynamic limit are explained using the isoperimetric quotient of the particles and the coordination number in the disordered fluid and in the assembled structure.

preprint2011arXiv

Degenerate Quasicrystal of Hard Triangular Bipyramids

We report a degenerate quasicrystal in Monte Carlo simulations of hard triangular bipyramids each composed of two regular tetrahedra sharing a single face. The dodecagonal quasicrystal is similar to that recently reported for hard tetrahedra [Haji-Akbari et al., Nature (London) 462, 773 (2009)] but degenerate in the pairing of tetrahedra, and self-assembles at packing fractions above 54%. This notion of degeneracy differs from the degeneracy of a quasiperiodic random tiling arising through phason flips. Free energy calculations show that a triclinic crystal is preferred at high packing fractions.

preprint2011arXiv

Phase Diagram of Hard Tetrahedra

Advancements in the synthesis of faceted nanoparticles and colloids have spurred interest in the phase behavior of polyhedral shapes. Regular tetrahedra have attracted particular attention because they prefer local symmetries that are incompatible with periodicity. Two dense phases of regular tetrahedra have been reported recently. The densest known tetrahedron packing is achieved in a crystal of triangular bipyramids (dimers) with packing density 4000/4671=85.63%. In simulation a dodecagonal quasicrystal is observed; its approximant, with periodic tiling (3.4.3^2.4), can be compressed to a packing fraction of 85.03%. Here, we show that the quasicrystal approximant is more stable than the dimer crystal for packing densities below 84% using Monte Carlo computer simulations and free energy calculations. To carry out the free energy calculations, we use a variation of the Frenkel-Ladd method for anisotropic shapes and thermodynamic integration. The enhanced stability of the approximant can be attributed to a network substructure, which maximizes the free volume (and hence the 'wiggle room') available to the particles and facilitates correlated motion of particles, which further contributes to entropy and leads to diffusion for packing densities below 65%. The existence of a solid-solid transition between structurally distinct phases not related by symmetry breaking -- the approximant and the dimer crystal-- is unusual for hard particle systems.

preprint2010arXiv

Disordered, Quasicrystalline and Crystalline Phases of Densely Packed Tetrahedra

All hard, convex shapes are conjectured by Ulam to pack more densely than spheres, which have a maximum packing fraction of ϕ = π/\sqrt18 ~ 0.7405. For many shapes, simple lattice packings easily surpass this packing fraction. For regular tetrahedra, this conjecture was shown to be true only very recently; an ordered arrangement was obtained via geometric construction with ϕ = 0.7786, which was subsequently compressed numerically to ϕ = 0.7820. Here we show that tetrahedra pack much better than this, and in a completely unexpected way. Following a conceptually different approach, using thermodynamic computer simulations that allow the system to evolve naturally towards high-density states, we observe that a fluid of hard tetrahedra undergoes a first-order phase transition to a dodecagonal quasicrystal, which can be compressed to a packing fraction of ϕ = 0.8324. By compressing a crystalline approximant of the quasicrystal, the highest packing fraction we obtain is ϕ = 0.8503. If quasicrystal formation is suppressed, the system remains disordered, jams, and compresses to ϕ = 0.7858. Jamming and crystallization are both preceded by an entropy-driven transition from a simple fluid of independent tetrahedra to a complex fluid characterized by tetrahedra arranged in densely packed local motifs that form a percolating network at the transition. The quasicrystal that we report represents the first example of a quasicrystal formed from hard or non-spherical particles. Our results demonstrate that particle shape and entropy can produce highly complex, ordered structures.