HOME

TheInfoList



OR:

Multiscale Green's function (MSGF) is a generalized and extended version of the classical
Green's function In mathematics, a Green's function is the impulse response of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions. This means that if \operatorname is the linear differentia ...
(GF) technique for solving mathematical equations. The main application of the MSGF technique is in modeling of
nanomaterials * Nanomaterials describe, in principle, materials of which a single unit is sized (in at least one dimension) between 1 and 100 nm (the usual definition of nanoscale). Nanomaterials research takes a materials science-based approach to nan ...
. These materials are very small – of the size of few
nanometer 330px, Different lengths as in respect to the molecular scale. The nanometre (international spelling as used by the International Bureau of Weights and Measures; SI symbol: nm) or nanometer (American and British English spelling differences#-re, ...
s. Mathematical modeling of nanomaterials requires special techniques and is now recognized to be an independent branch of science. A mathematical model is needed to calculate the displacements of atoms in a crystal in response to an applied static or time dependent force in order to study the mechanical and physical properties of nanomaterials. One specific requirement of a model for nanomaterials is that the model needs to be multiscale and provide seamless linking of different length scales.
Green's function In mathematics, a Green's function is the impulse response of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions. This means that if \operatorname is the linear differentia ...
(GF) was originally formulated by the British mathematical physicist George Green in the year 1828 as a general technique for solution of operator equations. It has been extensively used in mathematical
Physics Physics is the natural science that studies matter, its fundamental constituents, its motion and behavior through space and time, and the related entities of energy and force. "Physical science is that department of knowledge which relat ...
over the last almost two hundred years and applied to a variety of fields. Reviews of some applications of GFs such as for many body theory and
Laplace equation In mathematics and physics, Laplace's equation is a second-order partial differential equation named after Pierre-Simon Laplace, who first studied its properties. This is often written as \nabla^2\! f = 0 or \Delta f = 0, where \Delta = \nab ...
are available in the Wikipedia. The GF based techniques are used for modeling of various physical processes in materials such as
phonon In physics, a phonon is a collective excitation in a periodic, elastic arrangement of atoms or molecules in condensed matter, specifically in solids and some liquids. A type of quasiparticle, a phonon is an excited state in the quantum mechanical ...
s,
Electronic band structure In solid-state physics, the electronic band structure (or simply band structure) of a solid describes the range of energy levels that electrons may have within it, as well as the ranges of energy that they may not have (called ''band gaps'' or '' ...
and elastostatics.


Application of the MSGF method for modeling nanomaterials

The MSGF method is a relatively new GF technique for mathematical modeling of nanomaterials. Mathematical models are used for calculating the response of materials to an applied force in order to simulate their mechanical properties. The MSGF technique links different length scales in modeling of nanomaterials. Nanomaterials are of atomistic dimensions and need to be modeled at the length scales of nanometers. For example, a silicon
nanowire A nanowire is a nanostructure in the form of a wire with the diameter of the order of a nanometre (10−9 metres). More generally, nanowires can be defined as structures that have a thickness or diameter constrained to tens of nanometers or less ...
, whose width is about five nanometers, contains just 10 – 12 atoms across its width. Another example is
graphene Graphene () is an allotrope of carbon consisting of a single layer of atoms arranged in a hexagonal lattice nanostructure.
and many new two-dimensional (2D) solids. These new materials are ultimate in thinness because they are just one or two atoms thick. Multiscale modeling is needed for such materials because their properties are determined by the discreteness of their atomistic arrangements as well as their overall dimensions. The MSGF method is multiscale in the sense that it links the response of materials to an applied force at atomistic scales to their response at the macroscopic scales. The response of materials at the macroscopic scales is calculated by using the continuum model of solids. In the continuum model, the discrete atomistic structure of solids is averaged out into a continuum. Properties of nanomaterials are sensitive to their atomistic structure as well as their overall dimensions. They are also sensitive to the macroscopic structure of the host material in which they are embedded. The MSGF method is used to model such composite systems. The MSGF method is also used for analyzing behavior of crystals containing lattice defects such as vacancies, interstitials, or foreign atoms. Study of these lattice defects is of interest as they play a role in materials technology. Presence of a defect in a lattice displaces the host atoms from their original position or the lattice gets distorted. This is shown in Fig 1 for a 1D lattice as an example. Atomistic scale modeling is needed to calculate this distortion near the defect, whereas the continuum model is used to calculate the distortion far away from the defect. The MSGF links these two scales seamlessly.


MSGF for nanomaterials

The MSGF model of nanomaterials accounts for multiparticles as well as multiscales in materials. It is an extension of the lattice statics Green’s function (LSGF) method that was originally formulated at the Atomic Energy Research Establishment Harwell in U.K. in 1973. It is also referred to as Tewary method in the literature The LSGF method complements
molecular dynamics Molecular dynamics (MD) is a computer simulation method for analyzing the physical movements of atoms and molecules. The atoms and molecules are allowed to interact for a fixed period of time, giving a view of the dynamic "evolution" of the s ...
(MD) method for modeling multiparticle systems. The LSGF method is based upon the use of the Born von Karman (BvK) model and can be applied to different lattice structures and defects. The MSGF method is an extended version of the LSGF method and has been applied to many nanomaterials and 2D materials At the atomistic scales, a crystal or a crystalline solid is represented by a collection of interacting atoms located at discrete sites on a geometric lattice. A perfect crystal consists of a regular and periodic geometric lattice. The perfect lattice has translation symmetry, which means that all the unit cells are identical. In a perfect periodic lattice, which is assumed to be infinite, all atoms are identical. At equilibrium each atom is assumed to be located at its lattice site. The force at any atom due to other atoms just cancels out so the net force at each atom is zero. These conditions break down in a distorted lattice in which atoms get displaced from their positions of equilibrium. The lattice distortion may be caused by an externally applied force. The lattice can also be distorted by introducing a defect in the lattice or displacing an atom that disturbs the equilibrium configuration and induces a force on the lattice sites. This is shown in Fig. 1. The objective of the mathematical model is to calculate the resulting values of the atomic displacements. The GF in the MSGF method is calculated by minimizing the total energy of the lattice. The potential energy of the lattice in the form of an infinite Taylor series in atomic displacements in the harmonic approximation as follows : W = -\sum_ f_a(L)u_a(L) + \frac 1 2 \sum_ \sum_ K_(L,L') u_a(L) u_b(L') \qquad\qquad (1) where ''L'' and ''L''′ label the atoms, ''a'' and ''b'' denote the Cartesian coordinates, ''u'' denotes the atomic displacement, and −''f'' and ''K'' are the first and second coefficients in the Taylor series. They are defined by : f_a(L) = -\frac,\qquad\qquad(2) and : K_(L,L')= \frac, \qquad\qquad(3) where the derivatives are evaluated at zero displacements. The negative sign is introduced in the definition of ''f'' for convenience. Thus ''f''(''L'') is a 3D vector that denotes the force at the atom L. Its three Cartesian components are denoted by fa(L) where ''a'' = ''x'', ''y'', or ''z''. Similarly K(L,L’) is a 3x3 matrix, which is called the force- constant matrix between the atoms at L and L’. Its 9 elements are denoted by ''K''''ab''(''L'',''L''′) for ''a'', ''b'' = ''x'', ''y'', or ''z''. At equilibrium, the energy W is minimum. Accordingly, the first derivative of W with respect to each u must be zero. This gives the following relation from Eq. (1) : \sum_ K_(L,L')u_b(L') = f_a(L).\qquad\qquad(4) It can be shown by direct substitution that the solution of Eq. (4) can be written as : u_a(L) = \sum_ G_(L,L') f_b(L') \qquad\qquad(5) where G is defined by the following inversion relation : \sum_ K_(L,L'') G_(L'',L') = \delta(a,b) \delta(L,L'). \qquad\qquad(6) In Eq. (6), ''δ''(''m'',''n'') is the discrete delta function of two discrete variable m and ''n''. Similar to the case of
Dirac delta function In mathematics, the Dirac delta distribution ( distribution), also known as the unit impulse, is a generalized function or distribution over the real numbers, whose value is zero everywhere except at zero, and whose integral over the entire ...
for continuous variables, it is defined to be 1 if ''m'' = ''n'' and 0 otherwise. Equations (4)–(6) can be written in the matrix notation as follows: : \begin K u & = f & & (7) \\ u & = G f & & (8) \\ \text G & = K^ & & (9) \end The matrices ''K'' and ''G'' in the above equations are 3''N'' × 3''N'' square matrices and ''u'' and ''f'' are 3''N''-dimensional column vectors, where N is the total number of atoms in the lattice. The matrix ''G'' is the multiparticle GF and is referred to as the lattice statics Green’s function (LSGF),. If ''G'' is known, the atomic displacements for all atoms can be calculated by using Eq. (8). One of the main objectives of modeling is the calculation of the atomistic displacements u caused by an applied force f. The displacements, in principle, are given by Eq. (8). However, it involves inversion of the matrix K which is 3N x 3N. For any calculation of practical interest N ~ 10,000 but preferably a million for more realistic simulations. Inversion of such a large matrix is computationally extensive and special techniques are needed for the calculation of u’s. For regular periodic lattices, LSGF is one such technique. It consists of calculating G in terms of its Fourier transform, and is similar to the calculation of the phonon GF. The LSGF method has now been generalized to include the multiscale effects in the MSGF method. The MSGF method is capable of linking length scales seamlessly. This property has been used in developing a hybrid MSGF method that combines the GF and the MD methods and has been used for simulating less symmetric nanoinclusions such as quantum dots in semiconductors. For a perfect lattice without defects, the MSGF links directly the atomistic scales in LSGF to the macroscopic scales through the continuum model. A perfect lattice has full translation symmetry so all the atoms are equivalent. In this case any atom can be chosen as the origin and G(L,L') can be expressed by a single index (L'-L) defined as : G(L, L') = G(0, L-L') = G(L'-L) \qquad\qquad (10) The asymptotic limit of ''G''(''L''), that satisfies Eq. (10), for large values of ''R''(''L'') is given by : \lim_ (0,L)\rightarrow G_c(x) + O(1/x^4) \qquad\qquad (11) where ''x'' = ''R''(''L'') is the position vector of the atom ''L'', and ''G''''c''(''x'') is the continuum Green's function (CGF), which is defined in terms of the elastic constants and used in modeling of conventional bulk materials at macroscales. In Eq. (11), O(1/''x''''n'') is the standard mathematical notation for a term of order 1/''x''''n'' and higher. The magnitude of ''G''''c''(''x'') is O(1/''x''2). The LSGF ''G''(0,''L'') in this equation reduces smoothly and automatically to the CGF for large enough ''x'' as terms O(1/''x''4) become gradually small and negligible. This ensures the seamless linkage of the atomistic length scale to the macroscopic continuum scale. Equations (8) and (9) along with the limiting relation given by Eq. (11), form the basic equations for the MSGF. Equation (9) gives the LSGF, which is valid at the atomistic scales and Eq. (11) relates it to the CGF, which is valid at the macro continuum scales. This equation also shows that the LSGF reduces seamlessly to the CGF.


MSGF method for calculating the effect of defects and discontinuities in nanomaterials

If a lattice contains defects, its translation symmetry is broken. Consequently, it is not possible to express ''G'' in terms of a single distance variable ''R''(''L''). Hence Eq. (10) is not valid anymore and the correspondence between the LSGF and the CGF, needed for their seamless linking breaks down. In such cases the MSGF links the lattice and the continuum scales by using the following procedure: If p denotes the change in the matrix K, caused by the defect(s), the force constant matrix K* for the defective lattice is written as : K^* = K - p \qquad\qquad (12) As in the case for the perfect lattice in Eq. (9), the corresponding defect GF is defined as the inverse of the full ''K*'' matrix. Use of Eq. (12), then leads to the following Dyson’s equation for the defect LSGF: : G^* = G + GpG^* \qquad\qquad (13) The MSGF method consists of solving Eq. (13) for ''G*'' by using the matrix partitioning technique or double Fourier transform. Once ''G*'' is known, the displacement vector is given by the following GF equation similar to Eq. (8): u= G* f (14) Equation (14) gives the desired solution, that is, the atomic displacements or the lattice distortion caused by the force f. However, it does not show the linkage of the lattice and the continuum multiple scales, because Eqs. (10) and (11) are not valid for the defect LSGF G*. The linkage between the lattice and the continuum model in case of lattice with defects is achieved by using an exact transformation described below. Using Eq.(13), Eq. (14) can be written in the following exactly equivalent form: u = Gf + G p G* f . (15) Use of Eq. (14) again on the right hand side of Eq. (15) gives, u = G f* (16) where f* = f + p u. (17) Note that Eq. (17) defines an effective force f* such that Eqs. (14) and (16) are exactly equivalent. Equation (16) expresses the atomic displacements u in terms of G, the perfect LSGF even for lattices with defects. The effect of the defects is included exactly in f*. The LSGF G is independent of f or f* and reduces to the CGF asymptotically and smoothly as given in Eq. (11). The effective force f* can be determined in a separate calculation by using an independent method if needed, and the lattice statics or the continuum model can be used for G. This is the basis of a hybrid model that combines MSGF and MD for simulating a germanium quantum dot in a silicon lattice. Equation (16) is the master equation of the MSGF method. It is truly multiscale. All the discrete atomistic contributions are included in f*. The Green's function G can be calculated independently, which can be fully atomistic for nanomaterials or partly or fully continuum for macroscales to account for surfaces and interfaces in material systems as needed


GF as Physical Characteristic of Solids

Tewary, Quardokus and DelRio have suggested that Green's function is not just a mathematical artefact, but a physical characteristic of the solid, which can be measured by using
scanning probe microscopy Scan may refer to: Acronyms * Schedules for Clinical Assessment in Neuropsychiatry (SCAN), a psychiatric diagnostic tool developed by WHO * Shared Check Authorization Network (SCAN), a database of bad check writers and collection agency for ba ...
. Their suggestion is based upon the fact that any process of measurement on a system, requires quantifying the response of the system to an external probe and the Green’s function gives the total response of the system to an applied probe. For example, if we want to measure the characteristics of a spring, we fix it at one end and apply a force f or f* at the other end. A measurement of the stretching u of the spring will then give the value of the Green’s function of the spring by using Eq. 14 or 16. In this example, the applied force is the probe and the stretching of the spring is its response to the probe. If the measured values of the Green’s functions of a solid are available, it will accurately characterize the response of a solid for engineering applications. For this reason, the Green’s function is also called the response function.


Causal Green's Function Molecular Dynamics

An important application of the MSGF method is in modeling the temporal (time-dependent) processes in solids, especially in nanomaterials. This is needed in diverse applications like testing and characterization of materials, propagation of waves and heat in nanomaterials, and modeling of radiation damage in semiconductors. These processes need to be simulated over a wide range of times from
femtoseconds A femtosecond is a unit of time in the International System of Units (SI) equal to 10 or of a second; that is, one quadrillionth, or one millionth of one billionth, of a second. For context, a femtosecond is to a second as a second is to about 31. ...
to
nanoseconds A nanosecond (ns) is a unit of time in the International System of Units (SI) equal to one billionth of a second, that is, of a second, or 10 seconds. The term combines the SI prefix ''nano-'' indicating a 1 billionth submultiple of an SI unit ( ...
or even
microseconds A microsecond is a unit of time in the International System of Units (SI) equal to one millionth (0.000001 or 10−6 or ) of a second. Its symbol is μs, sometimes simplified to us when Unicode is not available. A microsecond is equal to 1000 n ...
which is a challenging multiscale problem for nanomaterials. It was shown by Tewary that the use of the causal Green’s functions in molecular dynamics can significantly accelerate the temporal convergence of molecular dynamics. The new method, called CGFMD (Causal Green's Function Molecular Dynamics) is the temporal equivalent of the MSGF and is based upon the use of
causal Causality (also referred to as causation, or cause and effect) is influence by which one event, process, state, or object (''a'' ''cause'') contributes to the production of another event, process, state, or object (an ''effect'') where the ca ...
or retarded Green's functions. It has been applied to simulate the propagation of ripples in graphene, where it has been shown that the CGFMD can model time scales over 6 to 9 orders of magnitude at the atomistic level. At least in some idealized cases such as propagation of ripples in graphene, the CGFMD can bridge the time scales from femto to microseconds. The CGFMD has been refined and further developed in papers by Coluci, Dantas and Tewary.}


See also

*
Microscale and macroscale models Microscale models form a broad class of computational models that simulate fine-scale details, in contrast with macroscale models, which amalgamate details into select categories. Microscale and macroscale models can be used together to understan ...


References

{{Reflist Generalized functions