We study numerically the thermal buckling transition in graphene membranes under compressive strain and clamped boundaries employing an atomistic quasiharmonic model. The numerical simulations combine three different Monte Carlo methods, local moves, collective wave moves and parallel tempering. We determine the static and the dynamical critical exponents by finite-size scaling, and the nonlinear response to a transverse force near the transition. The correlation length exponent and the nonlinear response are in good agreement with recent renormalization-group calculations of elastic membranes. Despite the applied strain, we find a dynamical critical exponent and diffusion exponent of height fluctuations at the transition close to the value for freestanding graphene, z = 2(1 + ) and alpha = /( + 1), where is the static roughening exponent, as obtained in a recent study with a phase-field crystal model.

We present a mesoscale field theory unifying the modeling of growth, elasticity, and dislocations in quasicrystals. The theory is based on the amplitudes entering their density-wave representation. We introduce a free energy functional for complex amplitudes and assume nonconserved dissipative dynamics to describe their evolution. Elasticity, including phononic and phasonic deformations, along with defect nucleation and motion, emerges self-consistently by prescribing only the symmetry of quasicrystals. Predictions on the formation of semicoherent interfaces and dislocation kinematics are given.

We investigate the growth of two-dimensional (2D) crystals on fluctuating surfaces using a phase-field crystal model that is relevant on atomic length and diffusive time scales. Motivated by recent experiments, which achieved unprecedented fast growth of large-size high-quality 2D crystals on liquid substrates, we uncover novel effects of liquid surfaces on microstructural ordering. We find that substrate fluctuations generate shortranged noise that speeds up crystallization and grain growth of the overlayer, surpassing that of free-standing system. Coupling to the liquid substrate fluctuations can also modulate local randomness, leading to intriguing disordered structures with hidden spatial order, i.e., disordered hyperuniformity. These results reveal the physical mechanisms underlying the fast growth of 2D crystals on liquid surfaces, and demonstrate a novel strategy for synthesizing disordered hyperuniform thin film structures.

We study the effects of thermal fluctuations and pinned boundaries in graphene membranes by using a phase-field crystal model with out-of-plane deformations. For sufficiently long times, the linear diffusive behavior of height fluctuations in systems with free boundaries becomes a saturation regime, while at intermediate times the behavior is still subdiffusive as observed experimentally. Under compression, we find mirror buckling fluctuations where the average height changes from above to below the pinned boundaries, with the average time between fluctuations diverging below a critical temperature corresponding to a thermally induced buckling transition. Near the transition, we find a nonlinear height response in agreement with recent renormalization-group calculations and observed in experiments on graphene membranes under an external transverse force with clamped boundaries.

In this paper a systematic examination of graphene/hexagonal boron nitride (g/hBN) bilayers is presented, through a recently developed two-dimensional phase field crystal model that incorporates out-of-plane deformations. The system parameters are determined by closely matching the stacking energies and heights of graphene/hBN bilayers to those obtained from existing quantum-mechanical density functional theory calculations. Out-of-plane deformations are shown to reduce the energies of inversion domain boundaries in hBN, and the coupling between graphene and hBN layers leads to a bilayer defect configuration consisting of an inversion boundary in hBN and a domain wall in graphene. Simulations of twisted bilayers reveal the structure, energy, and elastic properties of the corresponding Moiré patterns, and show a crossover, as the misorientation angle between the layers increases, from a well-defined hexagonal network of domain boundaries and junctions to smeared-out patterns. The transition occurs when the thickness of domain walls approaches the size of the Moiré patterns, and coincides with the peaks in the average von Mises and volumetric stresses of the bilayer.

We study the dynamical behavior of the mean-square displacement of height fluctuations of freestanding graphene using a phase-field-crystal model introduced recently. We find that the dynamic scaling behavior obtained numerically at long times is well described by the scaling theory of polymerized membranes. The critical exponent characterizing the power-law increase with time depends only on the equilibrium roughening exponent xi as alpha = xi/(1 + xi). For sufficiently long times it crosses over to linear behavior for finite-size systems. The critical exponent a is in good agreement with the anomalous diffusion exponent observed experimentally in graphene, suggesting this is a property that could also be observable in other two-dimensional crystalline materials.

Magnetic skyrmions are topologically distinct particles whose thermally activated motion could be used to implement probabilistic computing paradigms. While solid-liquid phase transitions in skyrmion lattices have been demonstrated, the behavior of a skyrmion liquid and the effects of pinning are largely unknown. Here we demonstrate the formation of a weakly pinned skyrmion liquid in a magnetic heterostructure. By inserting a Ru wedge layer at the ferromagnet/heavy metal interface we evaluate the dependence of skyrmion dynamics on the skyrmion size and density. Our experiments demonstrate that the diffusion of skyrmions is largest in dense liquids with small skyrmions. The thermal motion of skyrmions at room temperature easily overcomes the narrow distribution of pinning site energies in the granular film structure, satisfying a key requirement of probabilistic device architectures. Micromagnetic simulations support the findings and also reveal the existence of a thermally activated high-frequency collective oscillation.

Comprehensive investigations of crystalline systems often require methods bridging atomistic and continuum scales. In this context, coarse-grained mesoscale approaches are of particular interest as they allow the examination of large systems and time scales while retaining some microscopic details. The so-called phase-field crystal (PFC) model conveniently describes crystals at diffusive time scales through a continuous periodic field which varies on atomic scales and is related to the atomic number density. To go beyond the restrictive atomic length scales of the PFC model, a complex amplitude formulation was first developed by Goldenfeld et al (2005 Phys. Rev. E 72 020601). While focusing on length scales larger than the lattice parameter, this approach can describe crystalline defects, interfaces, and lattice deformations. It has been used to examine many phenomena including liquid/solid fronts, grain boundary energies, and strained films. This topical review focuses on this amplitude expansion of the PFC model and its developments. An overview of the derivation, connection to the continuum limit, representative applications, and extensions is presented. A few practical aspects, such as suitable numerical methods and examples, are illustrated as well. Finally, the capabilities and bounds of the model, current challenges, and future perspectives are addressed.

The velocity of dislocations is derived analytically to incorporate and predict the intriguing effects induced by the preferential solute segregation and Cottrell atmospheres in both two-dimensional and three-dimensional binary systems of various crystalline symmetries. The corresponding mesoscopic description of defect dynamics is constructed through the amplitude formulation of the phase-field crystal model, which has been shown to accurately capture elasticity and plasticity in a wide variety of systems. Modifications of the Peach-Koehler force as a result of solute concentration variations and compositional stresses are presented, leading to interesting new predictions of defect motion due to effects of Cottrell atmospheres. These include the deflection of dislocation glide paths, the variation of climb speed and direction, and the change or prevention of defect annihilation, all of which play an important role in determining the fundamental behaviors of complex defect network and dynamics. The analytic results are verified by numerical simulations.

In this paper, a two-dimensional phase field crystal model of graphene and hexagonal boron nitride (hBN) is extended to include out-of-plane deformations in stacked multilayer systems. As proof of principle, the model is shown analytically to reduce to standard models of flexible sheets in the small deformation limit. Applications to strained sheets, dislocation dipoles, and grain boundaries are used to validate the behavior of a single flexible graphene layer. For multilayer systems, parameters are obtained to match existing theoretical density functional theory calculations for graphene/graphene, hBN/hBN, and graphene/hBN bilayers. More precisely, it is shown that the parameters can be chosen to closely match the stacking energies and layer spacing calculated by Zhou et al. [Phys. Rev. B 92, 155438 (2015)]. Further validation of the model is presented in a study of rotated graphene bilayers and stacking boundaries. The flexibility of the model is illustrated by simulations that highlight the impact of complex microstructures in one layer on the other layer in a graphene/graphene bilayer.

A phase-field crystal model that spontaneously undergoes displacive phase transitions is introduced by explicitly studying a two-component two-dimensional square crystal reminiscent of a perovskite. When the intercomponent free energy is a simple polynomial, the crystal undergoes displacive transitions in the < 10 > and < 11 > directions. When the interaction is a correlation function, displacements in any direction can occur. This displacive phase-field crystal (DPFC) model maps to Landau-Ginzburg-Devonshire (LGD) theories for ferroelectrics, and the DPFC and LGD models are compared in terms of phase transitions and domain walls. Dynamical simulations of quadrijunctions were also performed and found stable spiraling quadrijunctions and unstable nonspiraling quadrijunctions. Last, domain coarsening across a small-angle grain boundary demonstrates multiple forms of non-mean-curvature-driven growth.

Using the phase field crystal model (PFC model), an analysis of slow and fast dynamics of solid-liquid interfaces in solidification and melting processes is presented. Dynamical regimes for cubic lattices invading metastable liquids (solidification) and liquids propagating into metastable crystals (melting) are described in terms of the evolving amplitudes of the density field. Dynamical equations are obtained for body-centered cubic (bcc) and face-centered cubic (fcc) crystal lattices in one- and two-mode approximations. A universal form of the amplitude equations is obtained for the three-dimensional dynamics for different crystal lattices and crystallographic directions. Dynamics of the amplitude's propagation for different lattices and PFC mode's approximations is qualitatively compared. The traveling-wave velocity is quantitatively compared with data of molecular dynamics simulation previously obtained by Mendelev et al. [Modell. Simul. Mater. Sci. Eng. 18, 074002 (2010)] for solidification and melting of the aluminum fcc lattice.

The phase-field crystal model in its amplitude equation approximation is shown to provide an accurate description of the deformation field in defected crystalline structures, as well as of dislocation motion. We analyze in detail the elastic distortion and stress regularization at a dislocation core and show how the Burgers vector density can be directly computed from the topological singularities of the phase-field amplitudes. Distortions arising from these amplitudes are then supplemented with non-singular displacements to enforce mechanical equilibrium. This allows for the consistent separation of plastic and elastic time scales in this framework. A finite element method is introduced to solve the combined amplitude and elasticity equations, which is applied to a few prototypical configurations in two spatial dimensions for a crystal of triangular lattice symmetry: i) the stress field induced by an edge dislocation with an analysis of how the amplitude equation regularizes stresses near the dislocation core, ii) the motion of a dislocation dipole as a result of its internal interaction, and iii) the shrinkage of a rotated grain. We also compare our results with those given by other extensions of classical elasticity theory, such as strain-gradient elasticity and methods based on the smoothing of Burgers vector densities near defect cores.

The ability to use external magnetic fields to influence the microstructure in polycrystalline materials has potential applications in microstructural engineering. To explore this potential and to understand the complex interactions between electromagnetic fields and solid-state matter transport we consider a phase-field-crystal model. Together with efficient and scalable numerical algorithms this allows the examination of the role that external magnetic fields play on the evolution of defect structures and grain boundaries, on diffusive timescales. Examples for planar and circular grain boundaries explain the essential atomistic processes and large scale simulations in 2D are used to obtain statistical data on grain growth under the influence of external fields.

We introduce a model which contains the essential elements to formulate and study antiferromagnetism, using the phase-field crystal framework. We focus on the question of how magneto-elastic coupling could lift the frustration in the two-dimensional hexagonal antiferromagnetic phase. Using simulations we observe a rich variety of different phases stable in this model. To characterize different phases we calculate the chiral order parameter and identify the scaling behavior of this order parameter. Furthermore, we observe that vortices appear and are stable close to the nonmagnetic defects. Finally, we studied the ferrimagnetic and spin-flop phase transition in the presence of an external magnetic field.

Crystal lattice deformations can be described microscopically by explicitly accounting for the position of atoms or macroscopically by continuum elasticity. In this work, we report on the description of continuous elastic fields derived from an atomistic representation of crystalline structures that also include features typical of the microscopic scale. Analytic expressions for strain components are obtained from the complex amplitudes of the Fourier modes representing periodic lattice positions, which can be generally provided by atomistic modeling or experiments. The magnitude and phase of these amplitudes, together with the continuous description of strains, are able to characterize crystal rotations, lattice deformations, and dislocations. Moreover, combined with the so-called amplitude expansion of the phase-field crystal model, they provide a suitable tool for bridging microscopic to macroscopic scales. This study enables the in-depth analysis of elasticity effects for macroscale and mesoscale systems taking microscopic details into account.

Atomically thin two-dimensional heterostructures are a promising, novel class of materials with ground-breaking properties. The possibility of choosing many constituent components and their proportions allows optimization of these materials to specific requirements. The wide adaptability comes with a cost of large parameter space making it hard to experimentally test all the possibilities. Instead, efficient computational modeling is needed. However, large range of relevant time and length scales related to physics of polycrystalline materials poses a challenge for computational studies. To this end, we present an efficient and flexible phase-field crystal model to describe the atomic configurations of multiple atomic species and phases coexisting in the same physical domain. We extensively benchmark the model for two-dimensional binary systems in terms of their elastic properties and phase boundary configurations and their energetics. As a concrete example, we demonstrate modeling lateral heterostructures of graphene and hexagonal boron nitride. We consider both idealized bicrystals and large-scale systems with random phase distributions. We find consistent relative elastic moduli and lattice constants, as well as realistic continuous interfaces and faceted crystal shapes. Zigzag-oriented interfaces are observed to display the lowest formation energy.

The formation and dynamics of a wide variety of binary two-dimensional ordered structures and superlattices are investigated through a phase field crystal model with sublattice ordering. Various types of binary ordered phases, the phase diagrams, and the grain growth dynamics and structural transformation processes, including the emergence of topological defects, are examined. The results are compared to the ordering and assembly of two-component colloidal systems. Two factors governing the binary phase ordering are identified, the coupling and competition between the length scales of two sublattices and the selection of average particle densities of two components. The control and variation of these two factors lead to the prediction of various complex binary ordered patterns, with different types of sublattice ordering for integer vs. noninteger ratios of sublattice length scales. These findings will enable further systematic studies of complex ordering and assembly processes of binary systems particularly binary colloidal crystals.

Using a minimal model based on the phase-field-crystal formalism, we study the coupling between the density and magnetization in ferromagnetic solids. Analytical calculations for the square phase in two dimensions are presented and the small deformation properties of the system are examined. Furthermore, numerical simulations are conducted to study the influence of an external magnetic field on various phase transitions, the anisotropic properties of the free energy functional, and the scaling behaviour of the growth of the magnetic domains in a crystalline solid. It is shown that the energy of the system can depend on the direction of the magnetic moments, with respect to the crystalline direction. Furthermore, the growth of the magnetic domains in a crystalline solid is studied and is shown that the growth of domains is in agreement with expected behaviour.

We address a three-dimensional, coarse-grained description of dislocation networks at grain boundaries between rotated crystals. The so-called amplitude expansion of the phase-field crystal model is exploited with the aid of finite element method calculations. This approach allows for the description of microscopic features, such as dislocations, while simultaneously being able to describe length scales that are orders of magnitude larger than the lattice spacing. Moreover, it allows for the direct description of extended defects by means of a scalar order parameter. The versatility of this framework is shown by considering both fcc and bcc lattice symmetries and different rotation axes. First, the specific case of planar, twist grain boundaries is illustrated. The details of the method are reported and the consistency of the results with literature is discussed. Then, the dislocation networks forming at the interface between a spherical, rotated crystal embedded in an unrotated crystalline structure, are shown. Although explicitly accounting for dislocations which lead to an anisotropic shrinkage of the rotated grain, the extension of the spherical grain boundary is found to decrease linearly over time in agreement with the classical theory of grain growth and recent atomistic investigations. It is shown that the results obtained for a system with bcc symmetry agree very well with existing results, validating the methodology. Furthermore, fully original results are shown for fcc lattice symmetry, revealing the generality of the reported observations.