Light scattering in random media with wavelength-scale structures : astronomical and industrial applications

Light scattering, or scattering and absorption of electromagnetic waves, is an important tool in all remote-sensing observations. In astronomy, the light scattered or absorbed by a distant object can be the only source of information. In Solar-system studies, the light-scattering methods are employed when interpreting observations of atmosphereless bodies such as asteroids, atmospheres of planets, and cometary or interplanetary dust. Our Earth is constantly monitored from artificial satellites at different wavelengths. With remote sensing of Earth the light-scattering methods are not the only source of information: there is always the possibility to make in situ measurements. The satellite-based remote sensing is, however, superior in the sense of speed and coverage if only the scattered signal can be reliably interpreted. The optical properties of many industrial products play a key role in their quality. Especially for products such as paint and paper, the ability to obscure the background and to reflect light is of utmost importance. High-grade papers are evaluated based on their brightness, opacity, color, and gloss. In product development, there is a need for computer-based simulation methods that could predict the optical properties and, therefore, could be used in optimizing the quality while reducing the material costs. With paper, for instance, pilot experiments with an actual paper machine can be very time- and resource-consuming. The light-scattering methods presented in this thesis solve rigorously the interaction of light and material with wavelength-scale structures. These methods are computationally demanding, thus the speed and accuracy of the methods play a key role. Different implementations of the discrete-dipole approximation are compared in the thesis and the results provide practical guidelines in choosing a suitable code. In addition, a novel method is presented for the numerical computations of orientation-averaged light-scattering properties of a particle, and the method is compared against existing techniques. Simulation of light scattering for various targets and the possible problems arising from the finite size of the model target are discussed in the thesis. Scattering by single particles and small clusters is considered, as well as scattering in particulate media, and scattering in continuous media with porosity or surface roughness. Various techniques for modeling the scattering media are presented and the results are applied to optimizing the structure of paper. However, the same methods can be applied in light-scattering studies of Solar-system regoliths or cometary dust, or in any remote-sensing problem involving light scattering in random media with wavelength-scale structures.

Interfacial shear rheology of highly confined glassy polymers

We present results on interfacial shear rheology measurements on Langmuir monolayers of two different polymers, poly(vinyl acetate) and poly(methyl methacrylate) as a function of surface concentration and temperature. While for the high glass transition poly(methyl methacrylate) polymer we find a systematic transition from a viscous dominated regime to an elastic dominated regime as surface concentration is increased, monolayers of the low glass transition polymer, poly(vinyl acetate), remain viscous even at very high surface concentrations. We further interpret the results in terms of the soft glassy rheology model of Sollich et al. P. Sollich, F. C. Lequeux, P. Hebraud and M. E. Cates, Phys. Rev. Lett., 1997, 78, 2020-2023] and provide evidence of possible reduction in glass transition temperatures in both poly(methyl methacrylate) and poly(vinyl acetate) monolayers due to finite size effects.

Role of multicritical points in addressing some key problems concerning condensed matter science

We illustrate the potential of using higher order critical points in the deeper understanding of several interesting problems of condensed matter science, e.g. critical adsorption, finite size effects, morphology of critical fluctuations, reversible aggregation of colloids, dynamics of the ordering process, etc.

The collective dynamics of self-propelled particles

We propose a method for the dynamic simulation of a collection of self-propelled particles in a viscous Newtonian fluid. We restrict attention to particles whose size and velocity are small enough that the fluid motion is in the creeping flow regime. We propose a simple model for a self-propelled particle, and extended the Stokesian Dynamics method to conduct dynamic simulations of a collection of such particles. In our description, each particle is treated as a sphere with an orientation vector p, whose locomotion is driven by the action of a force dipole Sp of constant magnitude S0 at a point slightly displaced from its centre. To simplify the calculation, we place the dipole at the centre of the particle, and introduce a virtual propulsion force Fp to effect propulsion. The magnitude F0 of this force is proportional to S0. The directions of Sp and Fp are determined by p. In isolation, a self-propelled particle moves at a constant velocity u0 p, with the speed u0 determined by S0. When it coexists with many such particles, its hydrodynamic interaction with the other particles alters its velocity and, more importantly, its orientation. As a result, the motion of the particle is chaotic. Our simulations are not restricted to low particle concentration, as we implement the full hydrodynamic interactions between the particles, but we restrict the motion of particles to two dimensions to reduce computation. We have studied the statistical properties of a suspension of self-propelled particles for a range of the particle concentration, quantified by the area fraction φa. We find several interesting features in the microstructure and statistics. We find that particles tend to swim in clusters wherein they are in close proximity. Consequently, incorporating the finite size of the particles and the near-field hydrodynamic interactions is of the essence. There is a continuous process of breakage and formation of the clusters. We find that the distributions of particle velocity at low and high φa are qualitatively different; it is close to the normal distribution at high φa, in agreement with experimental measurements. The motion of the particles is diffusive at long time, and the self-diffusivity decreases with increasing φa. The pair correlation function shows a large anisotropic build-up near contact, which decays rapidly with separation. There is also an anisotropic orientation correlation near contact, which decays more slowly with separation. Movies are available with the online version of the paper.

Study of Low Temperature Magnetic Properties of a Single Chain Magnet with Alternate Isotropic and Non-collinear Anisotropic Units

Here we study thermodynamic properties of an important class of single-chain magnets (SCMs), where alternate units are isotropic and anisotropic with anisotropy axes being non-collinear. This class of SCMs shows slow relaxation at low temperatures which results from the interplay of two different relaxation mechanisms, namely dynamical and thermal. Here anisotropy is assumed to be large and negative, as a result, anisotropic units behave like canted spins at low temperatures; but even then simple Ising-type model does not capture the essential physics of the system due to quantum mechanical nature of the isotropic units. We here show how statistical behavior of this class of SCMs can be studied using a transfer matrix (TM) method. We also, for the first time, discuss in detail how weak inter-chain interactions can be treated by a TM method. The finite size effect is also discussed which becomes important for low temperature dynamics. At the end of this paper, we apply this technique to study a real helical chain magnet.

Fidelity susceptibility of one-dimensional models with twisted boundary conditions

Recently it has been shown that the fidelity of the ground state of a quantum many-body system can be used todetect its quantum critical points (QCPs). If g denotes the parameter in the Hamiltonian with respect to which the fidelity is computed, we find that for one-dimensional models with large but finite size, the fidelity susceptibility chi(F) can detect a QCP provided that the correlation length exponent satisfies nu < 2. We then show that chi(F) can be used to locate a QCP even if nu >= 2 if we introduce boundary conditions labeled by a twist angle N theta, where N is the system size. If the QCP lies at g = 0, we find that if N is kept constant, chi(F) has a scaling form given by chi(F) similar to theta(-2/nu) f (g/theta(1/nu)) if theta << 2 pi/N. We illustrate this both in a tight-binding model of fermions with a spatially varying chemical potential with amplitude h and period 2q in which nu = q, and in a XY spin-1/2 chain in which nu = 2. Finally we show that when q is very large, the model has two additional QCPs at h = +/- 2 which cannot be detected by studying the energy spectrum but are clearly detected by chi(F). The peak value and width of chi(F) seem to scale as nontrivial powers of q at these QCPs. We argue that these QCPs mark a transition between extended and localized states at the Fermi energy. DOI: 10.1103/PhysRevB.86.245424

Cavitation in materials with distributed weak zones: Implications on the origin of brittle fracture in metallic glasses

Several experimental studies have shown that fracture surfaces in brittle metallic glasses (MGs) generally exhibit nanoscale corrugations which may be attributed to the nucleation and coalescence of nanovoids during crack propagation. Recent atomistic simulations suggest that this phenomenon is due to large spatial fluctuations in material properties in a brittle MG, which leads to void nucleation in regions of low atomic density and then catastrophic fracture through void coalescence. To explain this behavior, we propose a model of a heterogeneous solid containing a distribution of weak zones to represent a brittle MG. Plane strain continuum finite element analysis of cavitation in such an elastic-plastic solid is performed with the weak zones idealized as periodically distributed regions having lower yield strength than the background material. It is found that the presence of weak zones can significantly reduce the critical hydrostatic stress for the onset of cavitation which is controlled uniquely by the local yield properties of these zones. Also, the presence of weak zones diminishes the sensitivity of the cavitation stress to the volume fraction of a preexisting void. These results provide plausible explanations for the observations reported in recent atomistic simulations of brittle MGs. An analytical solution for a composite, incompressible elastic-plastic solid with a weak inner core is used to investigate the effect of volume fraction and yield strength of the core on the nature of cavitation bifurcation. It is shown that snap-cavitation may occur, giving rise to sudden formation of voids with finite size, which does not happen in a homogeneous plastic solid. (c) 2012 Elsevier Ltd. All rights reserved.

Turbulence in the two-dimensional Fourier-truncated Gross-Pitaevskii equation

We undertake a systematic, direct numerical simulation of the twodimensional, Fourier-truncated, Gross-Pitaevskii equation to study the turbulent evolutions of its solutions for a variety of initial conditions and a wide range of parameters. We find that the time evolution of this system can be classified into four regimes with qualitatively different statistical properties. Firstly, there are transients that depend on the initial conditions. In the second regime, powerlaw scaling regions, in the energy and the occupation-number spectra, appear and start to develop; the exponents of these power laws and the extents of the scaling regions change with time and depend on the initial condition. In the third regime, the spectra drop rapidly for modes with wave numbers k > kc and partial thermalization takes place for modes with k < kc; the self-truncation wave number kc(t) depends on the initial conditions and it grows either as a power of t or as log t. Finally, in the fourth regime, complete thermalization is achieved and, if we account for finite-size effects carefully, correlation functions and spectra are consistent with their nontrivial Berezinskii-Kosterlitz-Thouless forms. Our work is a natural generalization of recent studies of thermalization in the Euler and other hydrodynamical equations; it combines ideas from fluid dynamics and turbulence, on the one hand, and equilibrium and nonequilibrium statistical mechanics on the other.

Supersolid in a one-dimensional model of hard-core bosons

We study a system of hard-core boson on a one-dimensional lattice with frustrated next-nearest-neighbor hopping and nearest-neighbor interaction. At half filling, for equal magnitude of nearest- and next-nearest-neighbor hopping, the ground state of this system exhibits a first-order phase transition from a bond-ordered solid to a charge-density-wave solid as a function of the nearest- neighbor interaction. Moving away from half filling we investigate the system at incommensurate densities, where we find a supersolid phase which has concurrent off-diagonal long-range order and density-wave order which is unusual in a system of hard-core bosons in one dimension. Using the finite-size density-matrix renormalization group method, we obtain the complete phase diagram for this model.

Renyi entropies of free bosons on the torus and holography

We analytically evaluate the Renyi entropies for the two dimensional free boson CFT. The CFT is considered to be compactified on a circle and at finite temperature. The Renyi entropies S-n are evaluated for a single interval using the two point function of bosonic twist fields on a torus. For the case of the compact boson, the sum over the classical saddle points results in the Riemann-Siegel theta function associated with the A(n-1) lattice. We then study the Renyi entropies in the decompactification regime. We show that in the limit when the size of the interval becomes the size of the spatial circle, the entanglement entropy reduces to the thermal entropy of free bosons on a circle. We then set up a systematic high temperature expansion of the Renyi entropies and evaluate the finite size corrections for free bosons. Finally we compare these finite size corrections both for the free boson CFT and the free fermion CFT with the one-loop corrections obtained from bulk three dimensional handlebody spacetimes which have higher genus Riemann surfaces as its boundary. One-loop corrections in these geometries are entirely determined by quantum numbers of the excitations present in the bulk. This implies that the leading finite size corrections contributions from one-loop determinants of the Chern-Simons gauge field and the Dirac field in the dual geometry should reproduce that of the free boson and the free fermion CFT respectively. By evaluating these corrections both in the bulk and in the CFT explicitly we show that this expectation is indeed true.

Short-Time Beta Relaxation in Glass-Forming Liquids Is Cooperative in Nature

Temporal relaxation of density fluctuations in supercooled liquids near the glass transition occurs in multiple steps. Using molecular dynamics simulations for three model glass-forming liquids, we show that the short-time beta relaxation is cooperative in nature. Using finite-size scaling analysis, we extract a growing length scale associated with beta relaxation from the observed dependence of the beta relaxation time on the system size. We find, in qualitative agreement with the prediction of the inhomogeneous mode coupling theory, that the temperature dependence of this length scale is the same as that of the length scale that describes the spatial heterogeneity of local dynamics in the long-time alpha-relaxation regime.

Trans-Scale Mechanics: Looking For The Missing Links Between Continuum And Micro/Nanoscopic Reality

Problems involving coupled multiple space and time scales offer a real challenge for conventional frameworks of either particle or continuum mechanics. In this paper, four cases studies (shear band formation in bulk metallic glasses, spallation resulting from stress wave, interaction between a probe tip and sample, the simulation of nanoindentation with molecular statistical thermodynamics) are provided to illustrate the three levels of trans-scale problems (problems due to various physical mechanisms at macro-level, problems due to micro-structural evolution at macro/micro-level, problems due to the coupling of atoms/molecules and a finite size body at micro/nano-level) and their formulations. Accordingly, non-equilibrium statistical mechanics, coupled trans-scale equations and simultaneous solutions, and trans-scale algorithms based on atomic/molecular interaction are suggested as the three possible modes of trans-scale mechanics.

Local stresses and strains in axially cracked cylindrical shells

A perturbation solution is obtained for the local stress-strain fields in an axially cracked cylindrical shell. The tenth-order differential equations are used that take into account the transverse shear deformation. The perturbation of a curvature parameter, λ, is employed, where . The stress intensity factors for finite size cylindrical shells subjected to bending and internal pressure are evaluated. Sufficient accuracy can be obtained without using fine mesh sizes in regions near the crack tip. Also analyzed are the influence of cylinder diameter and shearing stiffness on bulging.

Study on cracked plates, shells and 3 dimensional bodies

This paper presents a summary of the authors' recent work in following areas: (1) The stress-strain fields at crack tip in Reissner's plate. (2) The calculations of the stress intensity factors in finite size plates. (3) The stress-strain fields at crack tip in Reissner's shell. (4) The calculations of the stress intensity factors and bulging coefficients in finite size spherical shells. (5) The stress-strain fields along crack tip in three dimensional body with surface crack. (6) The calculation of stress intensity factors in a plate with surface crack.

Study on cracked spherical shells

The local-global anatysis method is systematically extended to the fracture analysis of spherical shells. On the basis of the shallow shell theory, which takes into account transverse shear deformations, governing equations for cracked spherical shells expressed in displacement and stress functions f, F and φ are proposed, and then a general solution including Modes, Ⅰ, Ⅱ, Ⅲ for stress-strain fields at crack tip in a spherical shell is obtained, which plays the same role as Williams's expansion in plane elasticity. The numerical results for finite-size spherical shells under different boundary conditions have been obtained. Furthermore, the bulging factors are analyzed with regard to shearing stiffness and an approximate formula is given.