Disordered Energy Minimizing Configurations
|
|
- Clifton Osborne
- 5 years ago
- Views:
Transcription
1 . p. /4 Disordered Energy Minimizing Configurations Salvatore Torquato Department of Chemistry, Department of Physics, Princeton Institute for the Science and Technology of Materials, and Program in Applied & Computational Mathematics Princeton University Joint work with Ge Zhang and Frank Stillinger, Princeton University
2 Slow and Rapid Cooling of a Liquid Classical ground states are those classical particle configurations with minimal potential energy per particle.. p. 2/4
3 . p. 2/4 Slow and Rapid Cooling of a Liquid Classical ground states are those classical particle configurations with minimal potential energy per particle. super cooled liquid liquid Volume glass crystal rapid quench very slow cooling glass transition (Tg) Temperature freezing point (Tf) Typically, ground states are periodic with high crystallographic symmetries.
4 . p. 2/4 Slow and Rapid Cooling of a Liquid Classical ground states are those classical particle configurations with minimal potential energy per particle. super cooled liquid liquid Volume glass crystal rapid quench very slow cooling glass transition (Tg) Temperature freezing point (Tf) Typically, ground states are periodic with high crystallographic symmetries. Can classical ground states ever be disordered?
5 Definitions A point process in d-dimensional Euclidean space R d is a distribution of an infinite number of points in R d with configuration r,r 2,... with a well-defined number density ρ (number of points per unit volume). This is statistically described by the n-particle correlation function g n (r,,...,r n ). A lattice Λ in d-dimensional Euclidean space R d is the set of points that are integer linear combinations of d basis (linearly independent) vectors a i, i.e., {n a + n 2 a n d a d n,...,n d Z} The space R d can be geometrically divided into identical regions F called fundamental cells, each of which contains just one point. For example, in R 2 : Every lattice has a dual (or reciprocal) lattice Λ. A periodic point distribution in R d is a fixed but arbitrary configuration of N points (N ) in each fundamental cell of a lattice.. p. 3/4
6 . p. 4/4 Definitions For statistically homogeneous and isotropic point processes in R d at number density ρ, g 2 (r) is a nonnegative radial function that is proportional to the probability density of pair distances r. We call h(r) g 2 (r) the total correlation function. When there is no long-range order in the system, h(r) [or g 2 (r) ] in the large-r limit. We call a point process disordered if h(r) tends to zero sufficiently rapidly such that it is integrable over all space. The nonnegative structure factor S(k) is defined in terms of the Fourier transform of h(r), which we denote by h(k): where k denotes wavenumber. S(k) + ρ h(k), When there is no long-range order in the system, S(k) in the large-k limit, the dual-space analog of the aforementioned direct space condition. In some generalized sense, S(k) can be viewed as a probability density of pair distances of points in reciprocal space.
7 . p. 5/4 Pair Statistics for Spatially Uncorrelated and Ordered Point Processes Poisson Distribution (Ideal Gas) g 2 (r) S(k) r Lattice 2 3 k g 2 (r) S(k) 2 3 r 4π/ 3 8π/ 3 2π/ 3 k
8 . p. 6/4 Hyperuniform Point Processes Torquato and Stillinger, PRE (23) Points can represent molecules of a material, stars in a galaxy, or trees in a forest. Let Ω represent a spherical window of radius R in d-dimensional Euclidean space R d. R Ω R Ω Denote by σ 2 (R) N 2 (R) N(R) 2 the number variance.
9 . p. 6/4 Hyperuniform Point Processes Torquato and Stillinger, PRE (23) Points can represent molecules of a material, stars in a galaxy, or trees in a forest. Let Ω represent a spherical window of radius R in d-dimensional Euclidean space R d. R Ω R Ω Denote by σ 2 (R) N 2 (R) N(R) 2 the number variance. For a Poisson point pattern and many disordered point patterns, σ 2 (R) R d.
10 Hyperuniform Point Processes Torquato and Stillinger, PRE (23) Points can represent molecules of a material, stars in a galaxy, or trees in a forest. Let Ω represent a spherical window of radius R in d-dimensional Euclidean space R d. R Ω R Ω Denote by σ 2 (R) N 2 (R) N(R) 2 the number variance. For a Poisson point pattern and many disordered point patterns, σ 2 (R) R d. We call point patterns whose variance grows more slowly than R d (window volume) hyperuniform (infinite-wavelength density fluctuations vanish). This implies that structure factor S(k) for k. All perfect crystals and quasicrystals are hyperuniform such that σ 2 (R) R d : number variance grows like window surface area.. p. 6/4
11 . p. 7/4 DISORDERED HYPERUNIFORM STATES OF MATTER Disordered hyperuniform many-particle systems can be regarded to be new distinguishable states of disordered matter in that they (i) behave more like perfect crystals or quasicrystals in the manner in which they suppress large-scale density fluctuations, and yet are also like liquids and glasses because they are statistically isotropic structures with no Bragg peaks; (ii) can exist as both as equilibrium and quenched nonequilibrium phases; (iii) and, appear to be endowed with unique bulk physical properties. Understanding such states of matter, which turn out to have technological importance, require new theoretical tools.
12 DISORDERED HYPERUNIFORM STATES OF MATTER Disordered hyperuniform many-particle systems can be regarded to be new distinguishable states of disordered matter in that they (i) behave more like perfect crystals or quasicrystals in the manner in which they suppress large-scale density fluctuations, and yet are also like liquids and glasses because they are statistically isotropic structures with no Bragg peaks; (ii) can exist as both as equilibrium and quenched nonequilibrium phases; (iii) and, appear to be endowed with unique bulk physical properties. Understanding such states of matter, which turn out to have technological importance, require new theoretical tools. The hyperuniformity concept provides a unified means of categorizing and characterizing crystals, quasicrystals and such special disordered systems according to the degree to which they suppress large-scale density fluctuations.. p. 7/4
13 . p. 8/4 Outline Brief Review of Collective-Coordinate Simulations Ensemble Theory of Disordered Ground States
14 . p. 9/4 Creation of Disordered Hyperuniform Ground States Uche, Stillinger & Torquato, Phys. Rev. E 24 Batten, Stillinger & Torquato, Phys. Rev. E 28 Collective-Coordinate Simulations Consider a system of N particles with configuration r N r,r 2,... in a fundamental region Ω (under periodic boundary conditions). The complex collective density variable is defined by ρ(k) = NX exp(ik r j ) j= The real nonnegative structure factor is S(k) = ρ(k) 2 N
15 . p. 9/4 Creation of Disordered Hyperuniform Ground States Uche, Stillinger & Torquato, Phys. Rev. E 24 Batten, Stillinger & Torquato, Phys. Rev. E 28 Collective-Coordinate Simulations Consider a system of N particles with configuration r N r,r 2,... in a fundamental region Ω (under periodic boundary conditions). The complex collective density variable is defined by ρ(k) = NX exp(ik r j ) j= The real nonnegative structure factor is S(k) = ρ(k) 2 N Consider stable bounded pair potentials v(r) in R d with Fourier tranform ṽ(k). The total energy is Φ N (r N ) = X i<j v(r ij ) = N 2 Ω X ṽ(k)s(k) + constant k When ṽ(k) positive k K and zero otherwise, finding configurations in which S(k) is constrained to be zero where ṽ(k) has support is equivalent to globally minimizing Φ(r N ). These hyperuniform ground states are called stealthy and generally highly degenerate.
16 Creation of Disordered Hyperuniform Ground States Uche, Stillinger & Torquato, Phys. Rev. E 24 Batten, Stillinger & Torquato, Phys. Rev. E 28 Collective-Coordinate Simulations Consider a system of N particles with configuration r N r,r 2,... in a fundamental region Ω (under periodic boundary conditions). The complex collective density variable is defined by ρ(k) = NX exp(ik r j ) j= The real nonnegative structure factor is S(k) = ρ(k) 2 N Consider stable bounded pair potentials v(r) in R d with Fourier tranform ṽ(k). The total energy is Φ N (r N ) = X i<j v(r ij ) = N 2 Ω X ṽ(k)s(k) + constant k When ṽ(k) positive k K and zero otherwise, finding configurations in which S(k) is constrained to be zero where ṽ(k) has support is equivalent to globally minimizing Φ(r N ). These hyperuniform ground states are called stealthy and generally highly degenerate. Stealthy patterns can be tuned by varying the parameter χ: number of independently (real) constrained degrees of freedom to the total number of degrees of freedom d(n ).. p. 9/4
17 . p. /4 Creation of Disordered Stealthy Ground States Unconstrained Region K Exclusion Zone S=
18 . p. /4 Creation of Disordered Stealthy Ground States Unconstrained Region K Exclusion Zone S= One class of stealthy potentials involves the following power-law form: ṽ(k) = v ( k/k) n Θ(K k), where n is any whole number. The special case n = is just the simple step function.. v(r) d=3 m= m=2 m= r k In the large-system (thermodynamic) limit with n = and n = 4, we have the following large-r v(k) asymptotic behavior, respectively: v(r) cos(r) r 2 (n = ) ~.5.5 v(r) r 4 (n = 4) m= m=2 m=4
19 Creation of Disordered Stealthy Ground States Unconstrained Region K Exclusion Zone S= One class of stealthy potentials involves the following power-law form: ṽ(k) = v ( k/k) n Θ(K k), where n is any whole number. The special case n = is just the simple step function.. v(r) d=3 m= m=2 m=4 ~ v(k).5.5 m= m=2 m= r k In the large-system (thermodynamic) limit with n = and n = 4, we have the following large-r asymptotic behavior, respectively: v(r) cos(r) r 2 (n = ) v(r) r 4 (n = 4) While the specific forms of these stealthy potentials lead to the same convergent ground-state energies, this will not be the case for the pressure and other thermodynamic quantities.. p. /4
20 . p. /4 Creation of Disordered Stealthy Ground States In our previous simulations, we began with an initial random distribution of N points and then found the energy minimizing configurations (with extremely high precision) using optimization techniques.
21 . p. /4 Creation of Disordered Stealthy Ground States In our previous simulations, we began with an initial random distribution of N points and then found the energy minimizing configurations (with extremely high precision) using optimization techniques. For χ <.5, the stealthy ground states are degenerate, disordered and isotropic..5 S(k) k The success rate to achieve any disordered ground state is %, independent of the dimension. It turns out that χ = /2 is the limit at which there are no longer degrees of freedom that can be constrained independently.
22 . p. /4 Creation of Disordered Stealthy Ground States In our previous simulations, we began with an initial random distribution of N points and then found the energy minimizing configurations (with extremely high precision) using optimization techniques. For χ <.5, the stealthy ground states are degenerate, disordered and isotropic..5 S(k) k The success rate to achieve any disordered ground state is %, independent of the dimension. It turns out that χ = /2 is the limit at which there are no longer degrees of freedom that can be constrained independently. These exotic disordered ground states, which possess a hidden order on large length scales, can be viewed as many-particle analogs of spin liquids.
23 . p. 2/4 Creation of Disordered Stealthy Ground States Two Dimensions Three Dimensions (a) χ=.467 (b) χ =.47 2 g 2 (r) r / L x
24 Creation of Disordered Stealthy Ground States Two Dimensions Three Dimensions (a) χ=.467 (b) χ =.47 2 g 2 (r) r / L x Increasing χ induces strong short-range ordering, and eventually the system runs out of degrees of freedom that can be constrained and crystallizes. Maximum χ S(k) 2 3. p. 2/4
25 . p. 3/4 Ensemble Theory of Disordered Ground States Can one formulate an ensemble theory for stealthy degenerate ground states? Which ensemble? A challenging issue is the determination of the entropically favored states at each χ or ρ, i.e., weights associated with ground states. For some ensemble at fixed density ρ, the average energy per particle u for radial potentials in the thermodynamic limit is given by u Φ(rN ) N = ρ 2 where v ṽ(k = ). Z R d v(r)g 2 (r)dr = ρ 2 v 2 v(r = ) + 2(2π) d Z R d ṽ(k)s(k)dk.
26 . p. 3/4 Ensemble Theory of Disordered Ground States Can one formulate an ensemble theory for stealthy degenerate ground states? Which ensemble? A challenging issue is the determination of the entropically favored states at each χ or ρ, i.e., weights associated with ground states. For some ensemble at fixed density ρ, the average energy per particle u for radial potentials in the thermodynamic limit is given by u Φ(rN ) N = ρ 2 where v ṽ(k = ). Z R d v(r)g 2 (r)dr = ρ 2 v 2 v(r = ) + 2(2π) d Z R d ṽ(k)s(k)dk. Consider the same class of stealthy radial potential functions ṽ(k) in R d. Whenever particle configurations in R d exist such that S(k) is constrained to be its minimum value of zero where ṽ(k) has support, the system must be at its ground state or global energy minimum: u = ρ 2 v v(r = ) 2 Remark: Ground-state manifold, which is generally degenerate, is invariant to the specific form of the stealthy potential because the energy is a configuration-independent constant at fixed ρ.
27 Ensemble Theory of Disordered Ground States Can one formulate an ensemble theory for stealthy degenerate ground states? Which ensemble? A challenging issue is the determination of the entropically favored states at each χ or ρ, i.e., weights associated with ground states. For some ensemble at fixed density ρ, the average energy per particle u for radial potentials in the thermodynamic limit is given by u Φ(rN ) N = ρ 2 where v ṽ(k = ). Z R d v(r)g 2 (r)dr = ρ 2 v 2 v(r = ) + 2(2π) d Z R d ṽ(k)s(k)dk. Consider the same class of stealthy radial potential functions ṽ(k) in R d. Whenever particle configurations in R d exist such that S(k) is constrained to be its minimum value of zero where ṽ(k) has support, the system must be at its ground state or global energy minimum: u = ρ 2 v v(r = ) 2 Remark: Ground-state manifold, which is generally degenerate, is invariant to the specific form of the stealthy potential because the energy is a configuration-independent constant at fixed ρ. In the thermodynamic limit, parameter χ is related to the number density ρ in any dimension d via ρ χ = V (K) 2d(2π) d, where V (K) is the volume of a d-dimensional sphere of radius K. Remarks: We see that χ and ρ are inversely proportional to one another. Thus, for fixed K and d, as χ tends to zero, ρ tends to infinity, which corresponds counterintutively to the uncorrelated ideal-gas limit (Poisson distribution). As χ increases from zero, the density ρ decreases.. p. 3/4
28 . p. 4/4 Ensemble Theory of Disordered Ground States Any periodic crystal with a finite basis is a stealthy ground state for all positive χ up to a maximum χ max (ρ min ) determined by its first positive Bragg peak (minimal vector in Fourier space). Lemma: At fixed K, a configuration comprised of the union of m different stealthy ground-state configurations in R d with χ, χ 2,..., χ m, respectively, is itself stealthy with a value χ value given by " m # X χ =, χ i i= which is the harmonic mean of the χ i divided by m. These last two facts can be used to see how complex aperiodic patterns are possible ground states.
29 Ensemble Theory of Disordered Ground States Any periodic crystal with a finite basis is a stealthy ground state for all positive χ up to a maximum χ max (ρ min ) determined by its first positive Bragg peak (minimal vector in Fourier space). Lemma: At fixed K, a configuration comprised of the union of m different stealthy ground-state configurations in R d with χ, χ 2,..., χ m, respectively, is itself stealthy with a value χ value given by " m # X χ =, χ i i= which is the harmonic mean of the χ i divided by m. These last two facts can be used to see how complex aperiodic patterns are possible ground states. Table 2: Periodic stealthy ground states in R with K =. Structure χ max ρ min Integer lattice Periodic with n-particle basis n n 2π 2π = = ( )n Table 2: Periodic stealthy ground states in R 2 with K =. Structure χ max ρ min Kagomé crystal Honeycomb crystal Square lattice Triangular lattice π 3 = π 8π 2 = = π 2 = π 4 = π 2 = π = π 2 = p. 4/4
30 . p. 5/4 Ensemble Theory of Disordered Ground States Table 3: Periodic stealthy ground states in R 3 with K =. Structure χ max ρ min Pyrochlore crystal Diamond crystal Simple hexagonal lattice SC lattice π 4 = π 3 = π 3 2 = = π 3 3 π 9 = π 4 3 π 3 = = π 3 = π 3 3 HCP crystal = = π 3 π FCC lattice = = π 3 MCC lattice BCC lattice 2 2π 9 = π 3 = Table 4: Periodic stealthy ground states in R 4 with K =. Structure χ max ρ min Kag 4 crystal Dia 4 crystal Z 4 lattice D 4 lattice π 2 4 = π 3 = π 2 6 = π 3 = π 2 6 = π 3 = π 2 8 = π 3 =
31 . p. 6/4 Canonical Ensemble Theory of Disordered Ground States We consider the Gibbs canonical ensemble in which the partition function Z is a function of ρ and absolute temperature T. Our main interest is in a theory in the limit T, i.e., the entropically favored ground states in the canonical ensemble.
32 . p. 6/4 Canonical Ensemble Theory of Disordered Ground States We consider the Gibbs canonical ensemble in which the partition function Z is a function of ρ and absolute temperature T. Our main interest is in a theory in the limit T, i.e., the entropically favored ground states in the canonical ensemble. Ground-State Pressure and Isothermal Compressibility Energy Route: The pressure in the thermodynamic limit at T = can be obtained from the energy per particle (taking v = ) via the relation Therefore, for stealthy potentials, p = ρ 2 u ρ The isothermal compressibility κ T ρ ρ p p = ρ2 2. T «κ T = ρ 2. T. of such a system is
33 Canonical Ensemble Theory of Disordered Ground States We consider the Gibbs canonical ensemble in which the partition function Z is a function of ρ and absolute temperature T. Our main interest is in a theory in the limit T, i.e., the entropically favored ground states in the canonical ensemble. Ground-State Pressure and Isothermal Compressibility Energy Route: The pressure in the thermodynamic limit at T = can be obtained from the energy per particle (taking v = ) via the relation Therefore, for stealthy potentials, p = ρ 2 u ρ The isothermal compressibility κ T ρ ρ p p = ρ2 2. T «κ T = ρ 2. T. of such a system is Virial Route: An alternative route to the pressure is through the virial equation, which at T = is given by p = ρ2 2d s ()» = ρ2 2d Z r d dv dr g 2(r)dr Z f(k) h(k)dk R d f(k = ) + (2π) d where f(k) is the Fourier transform of f(r) rdv/dr, when it exists.. p. 6/4
34 . p. 7/4 Canonical Ensemble Theory of Disordered Ground States Consider a radial function u(r) in R d such that its Fourier transform ũ(k) exists with compact support in the interval [, ). Let w(r) = rdu/dr such that its Fourier transform w(k) exists. Proposition : The Fourier transform of the radial function w(r) in R d is given by w(k) = d ũ(k) dũ dk. Corollary: It immediately follows from Proposition that w(k) has the same support as ũ(k) and w(k = ) = d ũ(k = ). Substitution of the general form h(k) = Θ(K k) + P(k) (where P(k) has support for k > K) into virial expression and use of Proposition yields p = 2d» ρ 2 ṽ(k = ) + ρ (2π) d Z R d f(k)dk ρ (2π) d Z f(k) P(k)dk R d = ρ2 2. Remarks: This agrees with pressure obtained via energy route whenever f(k) exists. Since the latter is a stronger condition than the existence of ṽ(k), it is possible to devise ṽ(k) for which f(k) does not exist and hence a virial pressure that either diverges or is nonconvergent.
35 . p. 8/4 Pair Statistics of Disordered Ground States In the ground state of a stealthy potential, the structure factor attains its minimum value S(k) = for < k K, and hence has the general form S(k) = Θ(k K)[ + Q(k)], where Θ(x) = (, x <,, x. Therefore, ρ h(k) = f(k) + P(k), where f(k) = Θ(K k) and P(k) = Θ(k K) Q(k). Remark: The function f(k) is identical to the Mayer-f function for a Gibbs-ensemble hard-sphere system in direct space. It easily shown that the functions v(r) and P(r) are orthogonal to one another: Z R d v(r)p(r)dr =.
36 Pair Statistics of Disordered Ground States Behavior of the Pair Correlation Function Near the Origin Using the general form for S(k) for stealthy ground states and elementary properties of Fourier transforms, it is easy to show pair correlation function at the origin is given by Z g 2 (r = ) = 2dχ + 2d 2 χ K k d Q(k)dk Since g 2 (r) must be nonnegative for all r, we have the following integral condition on Q(k): 2d 2 χ Z Hence, this integral must be positive for χ 2d. K k d Q(k)dk 2dχ Moreover, if g 2 (r) has positive curvature at the origin, Q(k) must obey the additional integral condition: (d + 2) Z K k d+ Q(k)dk. Finally, when g 2 (r = ) =, the results above yield the equality Z K k d Q(k)dk = 2dχ 2d 2 χ and, because the curvature must be positive in this instance, the previous inequality must generally be obeyed.. p. 9/4
37 . p. 2/4 Pseudo-Hard Spheres in Fourier Space From previous considerations, we see that that an important contribution to S(k) is a simple hard-core step function Θ(k K), which can be viewed as a disordered hard-sphere system in Fourier space in the limit that χ tends to zero, i.e., as the number density ρ tends to infinity..5.5 S(k).5 g 2 (r) k 2 3 r That the structure factor must have the behavior S(k) Θ(k K), χ is perfectly reasonable; it is a perturbation about the ideal-gas limit in which S(k) = for all k.
38 Pseudo-Hard Spheres in Fourier Space From previous considerations, we see that that an important contribution to S(k) is a simple hard-core step function Θ(k K), which can be viewed as a disordered hard-sphere system in Fourier space in the limit that χ tends to zero, i.e., as the number density ρ tends to infinity..5.5 S(k).5 g 2 (r) k 2 3 r That the structure factor must have the behavior S(k) Θ(k K), χ is perfectly reasonable; it is a perturbation about the ideal-gas limit in which S(k) = for all k. Imagine carrying out a series expansion of S(k) about χ =. We make the ansatz that for sufficiently small χ, S(k) in the canonical ensemble for a stealthy potential can be mapped to g 2 (r) for an effective disordered hard-sphere system for sufficiently small density.. p. 2/4
39 Let us define Pseudo-Hard Spheres in Fourier Space H(k) ρ h(k) = h HS (r = k) There is an Ornstein-Zernike integral equation that defines the appropriate direct correlation function C(k): H(k) = C(k) + η H(k) C(k), where η = bχ is an effective packing fraction, where b is a constant to be determined. Therefore, H(r) = C(r) (2π) d η C(r). This mapping enables us to exploit the well-developed accurate theories of standard Gibbsian disordered hard spheres in direct space. To get an idea of the large-r asymptotics, consider case χ for any d: ρh(r) r d/2 J d/2(r) (χ ) ρh(r) r cos(r (d + )π/4) (r + ) (d+)/2. p. 2/4
40 . p. 22/4 Pair Statistics: Theory vs. Simulations.5 d=3, χ=.5.5 d=3, χ=. S(k).5 Theoretical Simulation S(k).5 Theoretical Simulation k k.5 χ=.5.5 χ=. g 2 (r).5 d=, Simulation d=2, Simulation d=3, Simulation d=, Theoretical d=2, Theoretical d=3, Theoretical 5 r g 2 (r).5 d=, Simulation d=2, Simulation d=3, Simulation d=, Theoretical d=2, Theoretical d=3, Theoretical 5 r
41 . p. 23/4 Ensemble Theory for Number Variance For an ensemble in d-dimensional Euclidean space R d : h Z i σ 2 (R) = N(R) + ρ h(r)α(r; R)dr R d α(r; R)- scaled intersection volume of 2 windows of radius R separated by r R r α(r;r) Spherical window of radius R.8.6 d=.4 d= r/(2r)
42 . p. 23/4 Ensemble Theory for Number Variance For an ensemble in d-dimensional Euclidean space R d : h Z i σ 2 (R) = N(R) + ρ h(r)α(r; R)dr R d α(r; R)- scaled intersection volume of 2 windows of radius R separated by r R r α(r;r) Spherical window of radius R.8.6 d=.4 d= r/(2r) For large R and under certain conditions, we can show h «R d «R d «R d i σ 2 (R) = 2 d φ A + B + o, D D D where A and B are the volume and Z surface-area coefficients: Z A = S(k = ) = + ρ h(r)dr, B(R) = c(d) h(r)rdr, R d R d D: microscopic length scale, φ: dimensionless density Hyperuniform: A =, B >
43 . p. 23/4 Ensemble Theory for Number Variance For an ensemble in d-dimensional Euclidean space R d : h Z i σ 2 (R) = N(R) + ρ h(r)α(r; R)dr R d α(r; R)- scaled intersection volume of 2 windows of radius R separated by r R r α(r;r) Spherical window of radius R.8.6 d=.4 d= r/(2r) For large R and under certain conditions, we can show h «R d «R d «R d i σ 2 (R) = 2 d φ A + B + o, D D D where A and B are the volume and Z surface-area coefficients: Z A = S(k = ) = + ρ h(r)dr, B(R) = c(d) h(r)rdr, R d R d Definition: D: microscopic length scale, φ: dimensionless density Hyperuniform: A =, B > B = lim L L R L B(R)dR Finding the point process that minimizes B is equivalent to finding the ground state of overlap potential, which is the dual of minimizer of the Epstein zeta function with s = (d + )/2.
44 . p. 24/4 σ 2 (R)/R d Number Variance σ 2 (R): Theory vs. Simulations χ=.5 d=, Simulation d=2, Simulation d=3, Simulation d=, Theoretical d=2, Theoretical d=3, Theoretical σ 2 (R)/R d χ=. d=, Simulation d=2, Simulation d=3, Simulation d=, Theoretical d=2, Theoretical d=3, Theoretical R R σ 2 (R)/R d χ=.2 d=, Simulation d=2, Simulation d=3, Simulation d=, Theoretical d=2, Theoretical d=3, Theoretical R 6 8
45 . p. 25/4 Number Variance σ 2 (R): Theory vs. Simulations Table 5: The coefficent B for various values of χ across dimensions. χ d = d = 2 d = Integer lattice.67 Triangular lattice.58 BCC lattice.245
46 . p. 26/4 Definitions: Nearest-Neighbor Functions H P (r)dr = Probability that a point of the point process lies at a distance between r and r + from another point of the point process. E P (r) = Probability of finding a spherical cavity of radius r centered at an arbitrary point of the point process empty of any other points. Exact series representation: E P (r) = + X ( ) k ρk k! k= H P (r) = E P r. Z R d g k+ (r k+ ) k+ Y j=2 Θ(r r j )dr j. Upper and lower bounds: E P (r) Z(r) E P (r) exp[ Z(r)] where Z(r) = ρs () Z r x d g 2 (x)dx, is the cumulative coordination number: expected number of points within a sphere of radius r.
47 . p. 27/4 Nearest-Neighbor Function: Theory vs. Simulations.5 d=2, χ=.5.5 d=2, χ=. E p (r).5 Simulation Theoretical Lower Bound Theoretical Upper Bound E p (r).5 Simulation Theoretical Lower Bound Theoretical Upper Bound 2 3 r r.5 d=3, χ=.5.5 d=3, χ=. E p (r).5 Simulation Theoretical Lower Bound Theoretical Upper Bound E p (r).5 Simulation Theoretical Lower Bound Theoretical Upper Bound r r
48 T. p. 28/4 Remark About Excited States For a system in equilibrium, the compressibility relation relates the isothermal compressibility κ T to the structure factor at k = : S() = ρk B Tκ T, where k B is Boltzmann s constant. Because κ T is is bounded according to κ T = ρ 2 for finite ρ, S() = because T =, which clearly must be the case for stealthy ground states. Now consider excited states infinitesimally close to the stealthy ground states, i.e., when temperature T is positive and infinitesimally small. Under the assumption that the structure of such excited states will be infintesimally near the ground-state configurations, we can estimate to an excellent approximation how S() varies with T for such excited states: S() c(d) χ T in units k B = v =. Morover, it is expected that this positive value of S() will be the uniform value of S(k) for k K. This behavior of S(k) has indeed been verified by molecular dynamics simulations in the canonical ensemble: S() d=2, χ=. Simulation Theoretical
49 . p. 29/4 CONCLUSIONS Hyperuniformity has connections to physics (e.g., ground states, quantum systems, random matrices, novel materials, etc.), mathematics (e.g., geometry and number theory), and biology. Disordered hyperuniform systems are new distinguishable states of matter somewhere between a liquid and a perfect crystal, and appear to be endowed with novel physical properties Disordered ground states are necessarily highly degenerate and their existence for radial pair potentials in low d is surprising. Devising a corresponding ensemble theory is highly nontrivial (many ensembles from which to choose) and until recently had been lacking. We have derived exact general results that apply to any ensemble of stealthy ground states, disordered or not. For the specific case of the canonical ensemble, our theory enables us to obtain predictive analytical results for the structure and thermodynamics of the entropically favored disordered ground states. Proving the existence of disordered stealthy ground states is an outstanding problem. Extensions of these results to compact spaces would be very interesting.
50 . p. 3/4 Hyperuniformity and Number Theory Averaging fluctuating quantity Λ(R) gives coefficient of interest: Λ = lim L L L Λ(R)dR
51 . p. 3/4 Hyperuniformity and Number Theory Averaging fluctuating quantity Λ(R) gives coefficient of interest: Λ = lim L We showed that for a lattice σ 2 (R) = q L L Λ(R)dR ( ) 2πR d [J d/2(qr)] 2, Λ = 2 d π d q Epstein zeta function for a lattice is defined by Z(s) = q q 2s, Re s > d/2. q Summand can be viewed as an inverse power-law potential. For q d+. lattices, minimizer of Z(d + ) is the lattice dual to the minimizer of Λ. Surface-area coefficient Λ provides useful way to rank order crystals, quasicrystals and special correlated disordered point patterns.
52 . p. 3/4 Quantifying Suppression of Density Fluctuations at Large Scales: D The surface-area coefficient Λ for some crystal, quasicrystal and disordered one-dimensional hyperuniform point patterns. Pattern Λ Integer Lattice / Step+Delta-Function g 2 3/6 =.875 Fibonacci Chain.2 Step-Function g 2 /4 =.25 Randomized Lattice / Zachary & Torquato (29)
53 . p. 32/4 Quantifying Suppression of Density Fluctuations at Large Scales: 2D The surface-area coefficient Λ for some crystal, quasicrystal and disordered two-dimensional hyperuniform point patterns. 2D Pattern Λ/φ /2 Triangular Lattice Square Lattice.564 Honeycomb Lattice Kagomé Lattice Penrose Tiling Step+Delta-Function g 2.62 Step-Function g Zachary & Torquato (29)
54 . p. 33/4 Quantifying Suppression of Density Fluctuations at Large Scales: 3D Contrary to conjecture that lattices associated with the densest sphere packings have smallest variance regardless of d, we have shown that for d = 3, BCC has a smaller variance than FCC.
55 . p. 33/4 Quantifying Suppression of Density Fluctuations at Large Scales: 3D Contrary to conjecture that lattices associated with the densest sphere packings have smallest variance regardless of d, we have shown that for d = 3, BCC has a smaller variance than FCC. Pattern Λ/φ 2/3 BCC Lattice FCC Lattice HCP Lattice SC Lattice.2892 Diamond Lattice.4892 Wurtzite Lattice.4284 Damped-Oscillating g Step+Delta-Function g Step-Function g Carried out analogous calculations in high d (Zachary & Torquato, 29), of importance in communications. Disordered point patterns may win in high d (Torquato & Stillinger, 26).
56 Hyperuniform Diffraction Patterns. p. 34/4
57 . p. 35/4 General Scaling Behaviors Hyperuniform particle distributions possess structure factors with a small-wavenumber scaling S(k) k α, α >, including the special case α = + for periodic crystals. Hence, number variance σ 2 (R) increases for large R asymptotically as (Zachary and Torquato, 2) σ 2 (R) R d ln R, α = R d α, α < R d, α > (R + ). Until recently, all known hyperuniform configurations pertained to α.
58 . p. 36/4 Targeted Spectra S k α k /2 S(k) k k 2 k 4 k K Configurations are ground states of an interacting many-particle system with up to four-body interactions.
59 . p. 37/4 Targeted Spectra S k α with α Uche, Stillinger & Torquato (26) Figure : One of them is for S(k) k 6 and other for S(k) k.
60 . p. 38/4 Targeted Spectra S k α with α < Zachary & Torquato (2) (a) (b) Figure 2: Both configurations exhibit strong local clustering of points and possess a highly irregular local structure; however, only one of them is hyperuniform (with S k /2 ).
61 . p. 39/4 Stealthy Hyperuniform Ground and Excited States These soft, long-ranged potentials in two and three dimensions exhibit unusual low-t behavior. Batten, Stillinger and Torquato, Phys. Rev. Lett. (29) Have used disordered stealthy materials to design photonic materials with large complete band gaps. Florescu, Torquato and Steinhardt, PNAS (29)
62 . p. 39/4 Stealthy Hyperuniform Ground and Excited States These soft, long-ranged potentials in two and three dimensions exhibit unusual low-t behavior. Batten, Stillinger and Torquato, Phys. Rev. Lett. (29) Have used disordered stealthy materials to design photonic materials with large complete band gaps. Florescu, Torquato and Steinhardt, PNAS (29) Disordered Hyperuniform Biological Systems
63 . p. 39/4 Stealthy Hyperuniform Ground and Excited States These soft, long-ranged potentials in two and three dimensions exhibit unusual low-t behavior. Batten, Stillinger and Torquato, Phys. Rev. Lett. (29) Have used disordered stealthy materials to design photonic materials with large complete band gaps. Florescu, Torquato and Steinhardt, PNAS (29) Disordered Hyperuniform Biological Systems Which is the Cucumber Image?
64 . p. 4/4 Pair Correlation (Radial Distribution) Function g 2 (r) g 2 (r) Probability density function associated with finding a point at a radial distance r from a given point. Expected cumulative coordination number at number density ρ Z(r) = ρs () r xd g 2 (x)dx: Expected number of points contained in a spherical region of radius r.
65 Pair Correlation (Radial Distribution) Function g 2 (r) g 2 (r) Probability density function associated with finding a point at a radial distance r from a given point. Expected cumulative coordination number at number density ρ Z(r) = ρs () r xd g 2 (x)dx: Expected number of points contained in a spherical region of radius r. Poisson Point Distribution in R 3 Maximally Random Jammed State in R 3 g 2 (r) g 2 (r) r 2 3 r. p. 4/4
66 Diffraction Pattern for Poisson Point Distribution (Ideal Gas) Ground states are highly degenerate. 2.8 Thermodynamic limit.6.4 S(k) k. p. 4/4
Disordered Hyperuniformity: Liquid-like Behaviour in Structural Solids, A New Phase of Matter?
Disordered Hyperuniformity: Liquid-like Behaviour in Structural Solids, A New Phase of Matter? Kabir Ramola Martin Fisher School of Physics, Brandeis University August 19, 2016 Kabir Ramola Disordered
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 23 Nov 2003
Local Density Fluctuations, Hyperuniformity, and Order Metrics arxiv:cond-mat/0311532v1 [cond-mat.stat-mech] 23 Nov 2003 Salvatore Torquato Department of Chemistry and Princeton Materials Institute, Princeton
More informationEquilibrium Phase Behavior and Maximally. Random Jammed State of Truncated Tetrahedra
Equilibrium Phase Behavior and Maximally Random Jammed State of Truncated Tetrahedra Duyu Chen,, Yang Jiao, and Salvatore Torquato,,,,, Department of Chemistry, Princeton University, Princeton, New Jersey
More informationarxiv: v1 [cond-mat.stat-mech] 19 Aug 2015
Ground states of stealthy hyperuniform potentials: II. Stacked-slider phases G. Zhang and F. H. Stillinger Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA arxiv:1508.04728v1
More informationDesigner Spin Systems via Inverse Statistical Mechanics
Designer Spin Systems via Inverse Statistical Mechanics Princeton University, Department of Chemistry, Princeton, New Jersey, USA IPAM Workshop at UCLA February 27, 2015 The Inverse Problem: Preliminaries
More informationUnexpected Density Fluctuations in Jammed Disordered Sphere Packings. Abstract
Unexpected Density Fluctuations in Jammed Disordered Sphere Packings Aleksandar Donev, 1, 2 Frank H. Stillinger, 3 1, 2, 3, and Salvatore Torquato 1 Program in Applied and Computational Mathematics, Princeton
More informationPhysics 541: Condensed Matter Physics
Physics 541: Condensed Matter Physics In-class Midterm Exam Wednesday, October 26, 2011 / 14:00 15:20 / CCIS 4-285 Student s Name: Instructions There are 23 questions. You should attempt all of them. Mark
More informationThe Phases of Hard Sphere Systems
The Phases of Hard Sphere Systems Tutorial for Spring 2015 ICERM workshop Crystals, Quasicrystals and Random Networks Veit Elser Cornell Department of Physics outline: order by disorder constant temperature
More informationHyperuniformity on the Sphere
Hyperuniformity on the Sphere Johann S. Brauchart j.brauchart@tugraz.at 14. June 2018 Week 2 Workshop MATRIX Research Program "On the Frontiers of High Dimensional Computation" MATRIX, Creswick [ 1 / 56
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 3, 3 March 2006 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationarxiv: v2 [cond-mat.stat-mech] 24 Aug 2014
Hyperuniformity of critical absorbing states Daniel Hexner and Dov Levine, Department of Physics, Technion, Haifa, Israel Initiative for the Theoretical Sciences - CUNY Graduate Center 65 Fifth Avenue,
More informationThe Second Virial Coefficient & van der Waals Equation
V.C The Second Virial Coefficient & van der Waals Equation Let us study the second virial coefficient B, for a typical gas using eq.v.33). As discussed before, the two-body potential is characterized by
More information1 Fluctuations of the number of particles in a Bose-Einstein condensate
Exam of Quantum Fluids M1 ICFP 217-218 Alice Sinatra and Alexander Evrard The exam consists of two independant exercises. The duration is 3 hours. 1 Fluctuations of the number of particles in a Bose-Einstein
More informationSelf-assembly of the simple cubic lattice with an isotropic potential
Self-assembly of the simple cubic lattice with an isotropic potential Mikael C. Rechtsman, 1 Frank H. Stillinger, 2 and Salvatore Torquato 2,3,4,5, * 1 Department of Physics, Princeton University, Princeton,
More informationLecture 11 - Phonons II - Thermal Prop. Continued
Phonons II - hermal Properties - Continued (Kittel Ch. 5) Low High Outline Anharmonicity Crucial for hermal expansion other changes with pressure temperature Gruneisen Constant hermal Heat ransport Phonon
More informationG : Statistical Mechanics
G5.651: Statistical Mechanics Notes for Lecture 8 I. DISTRIBUTION FUNCTIONS AND PERTURBATION THEORY A. General formulation Recall the expression for the configurational partition function: Z N = dr 1 dr
More informationBreakdown of Elasticity Theory for Jammed Hard-Particle Packings: Conical Nonlinear Constitutive Theory
Breakdown of Elasticity Theory for Jammed Hard-Particle Packings: Conical Nonlinear Constitutive Theory S. Torquato, 1,2 A. Donev 2,3, and F. H. Stillinger 1 Department of Chemistry, 1 Princeton Materials
More informationarxiv: v2 [cond-mat.dis-nn] 20 Aug 2018
Absence of hyperuniformity in amorphous hard-sphere packings of nonvanishing complexity M. J. Godfrey and M. A. Moore School of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK (Dated:
More informationT. Egami. Model System of Dense Random Packing (DRP)
Introduction to Metallic Glasses: How they are different/similar to other glasses T. Egami Model System of Dense Random Packing (DRP) Hard Sphere vs. Soft Sphere Glass transition Universal behavior History:
More informationV.C The Second Virial Coefficient & van der Waals Equation
V.C The Second Virial Coefficient & van der Waals Equation Let us study the second virial coefficient B, for a typical gas using eq.(v.33). As discussed before, the two-body potential is characterized
More informationcollisions inelastic, energy is dissipated and not conserved
LECTURE 1 - Introduction to Granular Materials Jamming, Random Close Packing, The Isostatic State Granular materials particles only interact when they touch hard cores: rigid, incompressible, particles
More informationLecture 4: Basic elements of band theory
Phys 769 Selected Topics in Condensed Matter Physics Summer 010 Lecture 4: Basic elements of band theory Lecturer: Anthony J. Leggett TA: Bill Coish 1 Introduction Most matter, in particular most insulating
More informationSolid State Physics 460- Lecture 5 Diffraction and the Reciprocal Lattice Continued (Kittel Ch. 2)
Solid State Physics 460- Lecture 5 Diffraction and the Reciprocal Lattice Continued (Kittel Ch. 2) Ewald Construction 2θ k out k in G Physics 460 F 2006 Lect 5 1 Recall from previous lectures Definition
More informationSupplemental Material Controlled Self-Assembly of Periodic and Aperiodic Cluster Crystals
Supplemental Material Controlled Self-Assembly of Periodic and Aperiodic Cluster Crystals Kobi Barkan, 1 Michael Engel, 2 and Ron Lifshitz 1, 3 1 Raymond and Beverly Sackler School of Physics and Astronomy,
More informationSupplemental Material Controlled Self-Assembly of Periodic and Aperiodic Cluster Crystals
Supplemental Material Controlled Self-Assembly of Periodic and Aperiodic Cluster Crystals Kobi Barkan, 1 Michael Engel, 2 and Ron Lifshitz 1, 3 1 Raymond and Beverly Sackler School of Physics and Astronomy,
More informationReview of Last Class 1
Review of Last Class 1 X-Ray diffraction of crystals: the Bragg formulation Condition of diffraction peak: 2dd sin θθ = nnλλ Review of Last Class 2 X-Ray diffraction of crystals: the Von Laue formulation
More informationPRINCIPLES OF PHYSICS. \Hp. Ni Jun TSINGHUA. Physics. From Quantum Field Theory. to Classical Mechanics. World Scientific. Vol.2. Report and Review in
LONDON BEIJING HONG TSINGHUA Report and Review in Physics Vol2 PRINCIPLES OF PHYSICS From Quantum Field Theory to Classical Mechanics Ni Jun Tsinghua University, China NEW JERSEY \Hp SINGAPORE World Scientific
More informationComplex entangled states of quantum matter, not adiabatically connected to independent particle states. Compressible quantum matter
Complex entangled states of quantum matter, not adiabatically connected to independent particle states Gapped quantum matter Z2 Spin liquids, quantum Hall states Conformal quantum matter Graphene, ultracold
More informationκ = f (r 0 ) k µ µ k ν = κk ν (5)
1. Horizon regularity and surface gravity Consider a static, spherically symmetric metric of the form where f(r) vanishes at r = r 0 linearly, and g(r 0 ) 0. Show that near r = r 0 the metric is approximately
More informationarxiv: v1 [cond-mat.soft] 20 Jun 2008
Accurate determination of crystal structures based on averaged local bond order parameters Wolfgang Lechner and Christoph Dellago Faculty of hysics, University of Vienna, Boltzmanngasse, 19 Vienna, Austria
More informationPhysics 212: Statistical mechanics II Lecture XI
Physics 212: Statistical mechanics II Lecture XI The main result of the last lecture was a calculation of the averaged magnetization in mean-field theory in Fourier space when the spin at the origin is
More informationFirst Problem Set for Physics 847 (Statistical Physics II)
First Problem Set for Physics 847 (Statistical Physics II) Important dates: Feb 0 0:30am-:8pm midterm exam, Mar 6 9:30am-:8am final exam Due date: Tuesday, Jan 3. Review 0 points Let us start by reviewing
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More information(# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble
Recall from before: Internal energy (or Entropy): &, *, - (# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble & = /01Ω maximized Ω: fundamental statistical quantity
More informationCrystal Structure. Dr Bindu Krishnan
Solid State Physics-1 Crystal Structure Dr Bindu Krishnan CRYSTAL LATTICE What is crystal (space) lattice? In crystallography, only the geometrical properties of the crystal are of interest, therefore
More informationGalaxies 626. Lecture 3: From the CMBR to the first star
Galaxies 626 Lecture 3: From the CMBR to the first star Galaxies 626 Firstly, some very brief cosmology for background and notation: Summary: Foundations of Cosmology 1. Universe is homogenous and isotropic
More informationNucleation rate (m -3 s -1 ) Radius of water nano droplet (Å) 1e+00 1e-64 1e-128 1e-192 1e-256
Supplementary Figures Nucleation rate (m -3 s -1 ) 1e+00 1e-64 1e-128 1e-192 1e-256 Calculated R in bulk water Calculated R in droplet Modified CNT 20 30 40 50 60 70 Radius of water nano droplet (Å) Supplementary
More informationStructure and Dynamics : An Atomic View of Materials
Structure and Dynamics : An Atomic View of Materials MARTIN T. DOVE Department ofearth Sciences University of Cambridge OXFORD UNIVERSITY PRESS Contents 1 Introduction 1 1.1 Observations 1 1.1.1 Microscopic
More informationSince their introduction in 1987, photonic band gap (PBG)
Designer disordered materials with large, complete photonic band gaps Marian Florescu a, Salvatore Torquato b,c, and Paul J. Steinhardt a,1 Departments of a Physics and b Chemistry, and c Princeton Center
More information4 Evolution of density perturbations
Spring term 2014: Dark Matter lecture 3/9 Torsten Bringmann (torsten.bringmann@fys.uio.no) reading: Weinberg, chapters 5-8 4 Evolution of density perturbations 4.1 Statistical description The cosmological
More informationHeterogenous Nucleation in Hard Spheres Systems
University of Luxembourg, Softmatter Theory Group May, 2012 Table of contents 1 2 3 Motivation Event Driven Molecular Dynamics Time Driven MD simulation vs. Event Driven MD simulation V(r) = { if r < σ
More informationSIMULATIONAL ANALYSIS OF GLASSES PREPARED VIA DIFFERENT INTERATOMIC POTENTIALS
Journal of Optoelectronics and Advanced Materials Vol. 7, No. 4, August 2005, p. 1915-1922 SIMULATIONAL ANALYSIS OF GLASSES PREPARED VIA DIFFERENT INTERATOMIC POTENTIALS Y. Sano, J. K ga, F. Yonezawa Department
More informationStructure of the First and Second Neighbor Shells of Water: Quantitative Relation with Translational and Orientational Order.
Structure of the First and Second Neighbor Shells of Water: Quantitative Relation with Translational and Orientational Order Zhenyu Yan, Sergey V. Buldyrev,, Pradeep Kumar, Nicolas Giovambattista 3, Pablo
More informationChapter 1. Crystal structure. 1.1 Crystal lattices
Chapter 1 Crystal structure 1.1 Crystal lattices We will concentrate as stated in the introduction, on perfect crystals, i.e. on arrays of atoms, where a given arrangement is repeated forming a periodic
More informationStatistical Mechanics
Franz Schwabl Statistical Mechanics Translated by William Brewer Second Edition With 202 Figures, 26 Tables, and 195 Problems 4u Springer Table of Contents 1. Basic Principles 1 1.1 Introduction 1 1.2
More informationPhys 412 Solid State Physics. Lecturer: Réka Albert
Phys 412 Solid State Physics Lecturer: Réka Albert What is a solid? A material that keeps its shape Can be deformed by stress Returns to original shape if it is not strained too much Solid structure
More informationSpin Superfluidity and Graphene in a Strong Magnetic Field
Spin Superfluidity and Graphene in a Strong Magnetic Field by B. I. Halperin Nano-QT 2016 Kyiv October 11, 2016 Based on work with So Takei (CUNY), Yaroslav Tserkovnyak (UCLA), and Amir Yacoby (Harvard)
More informationPhonons I - Crystal Vibrations (Kittel Ch. 4)
Phonons I - Crystal Vibrations (Kittel Ch. 4) Displacements of Atoms Positions of atoms in their perfect lattice positions are given by: R 0 (n 1, n 2, n 3 ) = n 10 x + n 20 y + n 30 z For simplicity here
More informationSuggestions for Further Reading
Contents Preface viii 1 From Microscopic to Macroscopic Behavior 1 1.1 Introduction........................................ 1 1.2 Some Qualitative Observations............................. 2 1.3 Doing
More informationGeneral formula for the efficiency of quantum-mechanical analog of the Carnot engine
General formula for the efficiency of quantum-mechanical analog of the Carnot engine Sumiyoshi Abe Department of Physical Engineering, Mie University, Mie 514-8507, Japan Abstract: An analog of the Carnot
More informationarxiv: v1 [cond-mat.soft] 31 Dec 2018
Spherical structure factor and classification of hyperuniformity on the sphere Anže Lošdorfer Božič Department of Theoretical Physics, Jožef Stefan Institute, SI-1000 Ljubljana, Slovenia arxiv:181.11763v1
More information16. GAUGE THEORY AND THE CREATION OF PHOTONS
6. GAUGE THEORY AD THE CREATIO OF PHOTOS In the previous chapter the existence of a gauge theory allowed the electromagnetic field to be described in an invariant manner. Although the existence of this
More informationSUPPLEMENTARY INFORMATION
A phase diagram for jammed matter Chaoming Song, Ping Wang, Hernán A. Makse Here, we describe the details of the calculations leading to the phase diagram of jammed matter. The main theoretical result
More informationMetropolis Monte Carlo simulation of the Ising Model
Metropolis Monte Carlo simulation of the Ising Model Krishna Shrinivas (CH10B026) Swaroop Ramaswamy (CH10B068) May 10, 2013 Modelling and Simulation of Particulate Processes (CH5012) Introduction The Ising
More informationClass 29: Reciprocal Space 3: Ewald sphere, Simple Cubic, FCC and BCC in Reciprocal Space
Class 29: Reciprocal Space 3: Ewald sphere, Simple Cubic, FCC and BCC in Reciprocal Space We have seen that diffraction occurs when, in reciprocal space, Let us now plot this information. Let us designate
More informationPolymer Solution Thermodynamics:
Polymer Solution Thermodynamics: 3. Dilute Solutions with Volume Interactions Brownian particle Polymer coil Self-Avoiding Walk Models While the Gaussian coil model is useful for describing polymer solutions
More informationLinear Theory of Evolution to an Unstable State
Chapter 2 Linear Theory of Evolution to an Unstable State c 2012 by William Klein, Harvey Gould, and Jan Tobochnik 1 October 2012 2.1 Introduction The simple theory of nucleation that we introduced in
More informationPHYS 771, Quantum Mechanics, Final Exam, Fall 2011 Instructor: Dr. A. G. Petukhov. Solutions
PHYS 771, Quantum Mechanics, Final Exam, Fall 11 Instructor: Dr. A. G. Petukhov Solutions 1. Apply WKB approximation to a particle moving in a potential 1 V x) = mω x x > otherwise Find eigenfunctions,
More informationMolecular Dynamics Simulations
Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline Long-range Interactions Ewald Sum Fast Multipole Method Spherically Truncated Coulombic Potential Speeding up Calculations SPaSM
More informationRenormalization Group for the Two-Dimensional Ising Model
Chapter 8 Renormalization Group for the Two-Dimensional Ising Model The two-dimensional (2D) Ising model is arguably the most important in statistical physics. This special status is due to Lars Onsager
More informationMP464: Solid State Physics Problem Sheet
MP464: Solid State Physics Problem Sheet 1 Write down primitive lattice vectors for the -dimensional rectangular lattice, with sides a and b in the x and y-directions respectively, and a face-centred rectangular
More informationPHASE TRANSITIONS IN SOFT MATTER SYSTEMS
OUTLINE: Topic D. PHASE TRANSITIONS IN SOFT MATTER SYSTEMS Definition of a phase Classification of phase transitions Thermodynamics of mixing (gases, polymers, etc.) Mean-field approaches in the spirit
More informationLattice spin models: Crash course
Chapter 1 Lattice spin models: Crash course 1.1 Basic setup Here we will discuss the basic setup of the models to which we will direct our attention throughout this course. The basic ingredients are as
More informationThe Effects of Inhomogeneities on the Universe Today. Antonio Riotto INFN, Padova
The Effects of Inhomogeneities on the Universe Today Antonio Riotto INFN, Padova Frascati, November the 19th 2004 Plan of the talk Short introduction to Inflation Short introduction to cosmological perturbations
More informationThe goal of equilibrium statistical mechanics is to calculate the diagonal elements of ˆρ eq so we can evaluate average observables < A >= Tr{Â ˆρ eq
Chapter. The microcanonical ensemble The goal of equilibrium statistical mechanics is to calculate the diagonal elements of ˆρ eq so we can evaluate average observables < A >= Tr{Â ˆρ eq } = A that give
More informationPhysics 211B : Problem Set #0
Physics 211B : Problem Set #0 These problems provide a cross section of the sort of exercises I would have assigned had I taught 211A. Please take a look at all the problems, and turn in problems 1, 4,
More informationPoint Sets and Dynamical Systems in the Autocorrelation Topology
Point Sets and Dynamical Systems in the Autocorrelation Topology Robert V. Moody and Nicolae Strungaru Department of Mathematical and Statistical Sciences University of Alberta, Edmonton Canada, T6G 2G1
More informationPhysics 562: Statistical Mechanics Spring 2002, James P. Sethna Prelim, due Wednesday, March 13 Latest revision: March 22, 2002, 10:9
Physics 562: Statistical Mechanics Spring 2002, James P. Sethna Prelim, due Wednesday, March 13 Latest revision: March 22, 2002, 10:9 Open Book Exam Work on your own for this exam. You may consult your
More informationTheoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 11
WiSe 22..23 Prof. Dr. A-S. Smith Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg Problem. Theoretische Physik 2:
More informationNon-magnetic states. The Néel states are product states; φ N a. , E ij = 3J ij /4 2 The Néel states have higher energy (expectations; not eigenstates)
Non-magnetic states Two spins, i and j, in isolation, H ij = J ijsi S j = J ij [Si z Sj z + 1 2 (S+ i S j + S i S+ j )] For Jij>0 the ground state is the singlet; φ s ij = i j i j, E ij = 3J ij /4 2 The
More informationIntroduction. Statistical physics: microscopic foundation of thermodynamics degrees of freedom 2 3 state variables!
Introduction Thermodynamics: phenomenological description of equilibrium bulk properties of matter in terms of only a few state variables and thermodynamical laws. Statistical physics: microscopic foundation
More informationPacking hyperspheres in high-dimensional Euclidean spaces
Packing hyperspheres in high-dimensional Euclidean spaces Monica Skoge, 1 Aleksandar Donev, 2,3 Frank H. Stillinger, 4 and Salvatore Torquato 2,3,4,5, * 1 Department of Physics, Princeton University, Princeton,
More informationMP464: Solid State Physics Problem Sheet
MP464: Solid State Physics Problem Sheet 1) Write down primitive lattice vectors for the -dimensional rectangular lattice, with sides a and b in the x and y-directions respectively, and a face-centred
More informationChapter 6. Phase transitions. 6.1 Concept of phase
Chapter 6 hase transitions 6.1 Concept of phase hases are states of matter characterized by distinct macroscopic properties. ypical phases we will discuss in this chapter are liquid, solid and gas. Other
More informationA variational approach to Ising spin glasses in finite dimensions
. Phys. A: Math. Gen. 31 1998) 4127 4140. Printed in the UK PII: S0305-447098)89176-2 A variational approach to Ising spin glasses in finite dimensions R Baviera, M Pasquini and M Serva Dipartimento di
More informationPhase transitions and critical phenomena
Phase transitions and critical phenomena Classification of phase transitions. Discontinous (st order) transitions Summary week -5 st derivatives of thermodynamic potentials jump discontinously, e.g. (
More informationMeasuring Ellipsoids 1
Measuring Ellipsoids 1 Igor Rivin Temple University 2 What is an ellipsoid? E = {x E n Ax = 1}, where A is a non-singular linear transformation of E n. Remark that Ax = Ax, Ax = x, A t Ax. The matrix Q
More informationA Variational Analysis of a Gauged Nonlinear Schrödinger Equation
A Variational Analysis of a Gauged Nonlinear Schrödinger Equation Alessio Pomponio, joint work with David Ruiz Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Variational and Topological
More informationBasic Principles of Quasicrystallography
Chapter 2 Basic Principles of Quasicrystallography Conventional crystals are constructed by a periodic repetition of a single unit-cell. This periodic long-range order only agrees with two-, three-, four-
More informationGeneral Formula for the Efficiency of Quantum-Mechanical Analog of the Carnot Engine
Entropy 2013, 15, 1408-1415; doi:103390/e15041408 Article OPEN ACCESS entropy ISSN 1099-4300 wwwmdpicom/journal/entropy General Formula for the Efficiency of Quantum-Mechanical Analog of the Carnot Engine
More informationSPACETIME FROM ENTANGLEMENT - journal club notes -
SPACETIME FROM ENTANGLEMENT - journal club notes - Chris Heinrich 1 Outline 1. Introduction Big picture: Want a quantum theory of gravity Best understanding of quantum gravity so far arises through AdS/CFT
More informationThe existence of Burnett coefficients in the periodic Lorentz gas
The existence of Burnett coefficients in the periodic Lorentz gas N. I. Chernov and C. P. Dettmann September 14, 2006 Abstract The linear super-burnett coefficient gives corrections to the diffusion equation
More informationExam TFY4230 Statistical Physics kl Wednesday 01. June 2016
TFY423 1. June 216 Side 1 av 5 Exam TFY423 Statistical Physics l 9. - 13. Wednesday 1. June 216 Problem 1. Ising ring (Points: 1+1+1 = 3) A system of Ising spins σ i = ±1 on a ring with periodic boundary
More information8.333: Statistical Mechanics I Problem Set # 5 Due: 11/22/13 Interacting particles & Quantum ensembles
8.333: Statistical Mechanics I Problem Set # 5 Due: 11/22/13 Interacting particles & Quantum ensembles 1. Surfactant condensation: N surfactant molecules are added to the surface of water over an area
More informationAs in my QM3 and QM4, it is useful to replace the phase by its gradient: In the present, nonlinear case the hydrodynamic equations take a form:
1 Lecture 3 BEC What happens if we subject the condensate to an action of a time dependent potential? In the mean field approximation we may expect that the single particle wave function will satisfy the
More informationScaling Theory. Roger Herrigel Advisor: Helmut Katzgraber
Scaling Theory Roger Herrigel Advisor: Helmut Katzgraber 7.4.2007 Outline The scaling hypothesis Critical exponents The scaling hypothesis Derivation of the scaling relations Heuristic explanation Kadanoff
More informationComputer simulation methods (1) Dr. Vania Calandrini
Computer simulation methods (1) Dr. Vania Calandrini Why computational methods To understand and predict the properties of complex systems (many degrees of freedom): liquids, solids, adsorption of molecules
More information424 Index. Eigenvalue in quantum mechanics, 174 eigenvector in quantum mechanics, 174 Einstein equation, 334, 342, 393
Index After-effect function, 368, 369 anthropic principle, 232 assumptions nature of, 242 autocorrelation function, 292 average, 18 definition of, 17 ensemble, see ensemble average ideal,23 operational,
More information5. Systems in contact with a thermal bath
5. Systems in contact with a thermal bath So far, isolated systems (micro-canonical methods) 5.1 Constant number of particles:kittel&kroemer Chap. 3 Boltzmann factor Partition function (canonical methods)
More informationStatistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany
Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical
More informationWhat is a Quasicrystal?
July 23, 2013 Rotational symmetry An object with rotational symmetry is an object that looks the same after a certain amount of rotation. Rotational symmetry An object with rotational symmetry is an object
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More informationOn Schrödinger equations with inverse-square singular potentials
On Schrödinger equations with inverse-square singular potentials Veronica Felli Dipartimento di Statistica University of Milano Bicocca veronica.felli@unimib.it joint work with Elsa M. Marchini and Susanna
More informationFrom the ideal gas to an ideal glass: a thermodynamic route to random close packing.
Leslie V. Woodcock From the ideal gas to an ideal glass: a thermodynamic route to random close packing. Department of Chemical and Biomolecular Engineering National University of Singapore, Singapore 117576
More informationRandom closest packing in a 2D lattice model
J. Phys. A: Math. Gen. 33 (2000) 1729 1734. Printed in the UK PII: S0305-4470(00)09254-4 Random closest packing in a 2D lattice model E Eisenberg and A Baram Department of Physics, Bar-Ilan University,
More informationCollective behavior, from particles to fields
978-0-51-87341-3 - Statistical Physics of Fields 1 Collective behavior, from particles to fields 1.1 Introduction One of the most successful aspects of physics in the twentieth century was revealing the
More informationMarkov Chains and Stochastic Sampling
Part I Markov Chains and Stochastic Sampling 1 Markov Chains and Random Walks on Graphs 1.1 Structure of Finite Markov Chains We shall only consider Markov chains with a finite, but usually very large,
More informationPhase transitions and coarse graining for a system of particles in the continuum
Phase transitions and coarse graining for a system of particles in the continuum Elena Pulvirenti (joint work with E. Presutti and D. Tsagkarogiannis) Leiden University April 7, 2014 Symposium on Statistical
More informationPreliminaries. Growth of Structure. Today s measured power spectrum, P(k) Simple 1-D example of today s P(k) Growth in roughness: δρ/ρ. !(r) =!!
Growth of Structure Notes based on Teaching Company lectures, and associated undergraduate text with some additional material added. For a more detailed discussion, see the article by Peacock taken from
More informationDuality of multiparameter Hardy spaces H p on spaces of homogeneous type
Duality of multiparameter Hardy spaces H p on spaces of homogeneous type Yongsheng Han, Ji Li, and Guozhen Lu Department of Mathematics Vanderbilt University Nashville, TN Internet Analysis Seminar 2012
More information