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.

On the criteria for reverse transition in a two-dimensional boundary layer flow

Experiments on reverse transition were conducted in two-dimensional accelerated incompressible turbulent boundary layers. Mean velocity profiles, longitudinal velocity fluctuations $\tilde{u}^{\prime}(=(\overline{u^{\prime 2}})^{\frac{1}{2}})$ and the wall-shearing stress (TW) were measured. The mean velocity profiles show that the wall region adjusts itself to laminar conditions earlier than the outer region. During the reverse transition process, increases in the shape parameter (H) are accompanied by a decrease in the skin friction coefficient (Cf). Profiles of turbulent intensity (u’2) exhibit near similarity in the turbulence decay region. The breakdown of the law of the wall is characterized by the parameter \[ \Delta_p (=\nu[dP/dx]/\rho U^{*3}) = - 0.02, \] where U* is the friction velocity. Downstream of this region the decay of $\tilde{u}^{\prime}$ fluctuations occurred when the momentum thickness Reynolds number (R) decreased roughly below 400.

Stress induced phase transition in a monomolecular perfluroalkylsilane film self assembled on aluminium surface

We control the stiffnesses of two dual double cantelevers placed in series to control penetration into a perflurooctyltrichlorosilane monolayer self assembled on aluminium and silicon substrates. The top cantilever which carries the probe is displaced with respect to the bottom cantilever which carries the substrate, the difference in displacement recorded using capacitors gives penetration. We further modulate the input displacement sinusoidally to deconvolute the viscoelastic properties of the monolayer. When the intervention is limited to the terminal end of the molecule there is a strong viscous response in consonance with the ability of the molecule to dissipate energy by the generation of gauche defects freely. When the intervention reaches the backbone, at a contact mean pressure of 0.2GPa the damping disappears abruptly and the molecule registers a steep rise in elastic modulus and relaxation time constant, with increasing contact pressure. We offer a physical explanation of the process and describe this change as due to a phase transition from a liquid like to a solid like state.

Particle dynamics in a turbulent particle–gas suspension at high Stokes number. Part 1. Velocity and acceleration distributions

The effect of fluid velocity fluctuations on the dynamics of the particles in a turbulent gas–solid suspension is analysed in the low-Reynolds-number and high Stokes number limits, where the particle relaxation time is long compared with the correlation time for the fluid velocity fluctuations, and the drag force on the particles due to the fluid can be expressed by the modified Stokes law. The direct numerical simulation procedure is used for solving the Navier–Stokes equations for the fluid, the particles are modelled as hard spheres which undergo elastic collisions and a one-way coupling algorithm is used where the force exerted by the fluid on the particles is incorporated, but not the reverse force exerted by the particles on the fluid. The particle mean and root-mean-square (RMS) fluctuating velocities, as well as the probability distribution function for the particle velocity fluctuations and the distribution of acceleration of the particles in the central region of the Couette (where the velocity profile is linear and the RMS velocities are nearly constant), are examined. It is found that the distribution of particle velocities is very different from a Gaussian, especially in the spanwise and wall-normal directions. However, the distribution of the acceleration fluctuation on the particles is found to be close to a Gaussian, though the distribution is highly anisotropic and there is a correlation between the fluctuations in the flow and gradient directions. The non-Gaussian nature of the particle velocity fluctuations is found to be due to inter-particle collisions induced by the large particle velocity fluctuations in the flow direction. It is also found that the acceleration distribution on the particles is in very good agreement with the distribution that is calculated from the velocity fluctuations in the fluid, using the Stokes drag law, indicating that there is very little correlation between the fluid velocity fluctuations and the particle velocity fluctuations in the presence of one-way coupling. All of these results indicate that the effect of the turbulent fluid velocity fluctuations can be accurately represented by an anisotropic Gaussian white noise.

Structural Relaxation and Crystallization in a Pd40Cu30Ni10P20 Bulk Metallic Glass

The stability of a Pd40Cu30Ni10P20 bulk metallic glass (BMG) against structural relaxation is investigated by isothermal and isochronal annealing heat treatments below and above its glass transition temperature, Tg, for varying periods. Differential scanning calorimetry (DSC) of the annealed samples shows an excess endotherm at Tg, irrespective of the annealing temperature. This recovery peak evolves exponentially with annealing time and is due to the destruction of anneal-induced compositional short range ordering. The alloy exhibits a high resistance to crystallization on annealing below Tg and complex Pd- and Ni-phosphides evolve on annealing above Tg.

Interface states of laser ablated BaTiO3 and Ba0.9Ca0.1TiO3 thin films in MFS structure determined by DLTS and C – V technique

BaTiO3 and Ba0.9Ca0.1TiO3 thin films were deposited on the p – type Si substrate by pulsed excimer laser ablation technique. The Capacitance – Voltage (C-V) measurement measured at 1 MHz exhibited a clockwise rotating hysteresis loop with a wide memory window for the Metal – Ferroelectric – Semiconductor (MFS) capacitor confirming the ferroelectric nature. The low frequency C – V measurements exhibited the response of the minority carriers in the inversion region while at 1 MHz the C – V is of a high frequency type with minimum capacitance in the inversion region. The interface states of both the MFS structures were calculated from the Castagne – Vaipaille method (High – low frequency C – V curve). Deep Level Transient Spectroscopy (DLTS) was used to analyze the interface traps and capture cross section present in the MFS capacitor. There were distinct peaks present in the DLTS spectrum and these peaks were attributed to the presence of the discrete interface states present at the semiconductor – ferroelectric interface. The distribution of calculated interface states were mapped with the silicon energy band gap for both the undoped and Ca doped BaTiO3 thin films using both the C – V and DLTS method. The interface states of the Ca doped BaTiO3 thin films were found to be higher than the pure BaTiO3 thin films.

Hydrogen Passivation of Shallow Dopants in InP Studied by Photoluminescence Spectroscopy

The effect of hydrogenation on the photoluminescence (PL) of InP : Mg, InP : Zn and undoped n-InP is presented. An increase in the near band edge pl intensity due to passivation of non-radiative centers was observed in all the samples. A donor - acceptor pair transition was observed before hydrogenation in the InP : Mg sample and after hydrogenation in the InP : Zn sample due to the acceptor deactivation. In n-InP the enhancement of donor bound exciton after hydrogenation points to the absence of donor passivation.

On the solution of the problem of scattering of surface water waves by a sharp discontinuity in the surface boundary conditions

Closed-form analytical expressions are derived for the reflection and transmission coefficients for the problem of scattering of surface water waves by a sharp discontinuity in the surface-boundary-conditions, for the case of deep water. The method involves the use of the Havelock-type expansion of the velocity potential along with an analysis to solve a Carleman-type singular integral equation over a semi-infinite range. This method of solution is an alternative to the Wiener-Hopf technique used previously.

Axisymmetric turbulent boundary layer along a circular cylinder at constant pressure

A constant-pressure axisymmetric turbulent boundary layer along a circular cylinder of radius a is studied at large values of the frictional Reynolds number a+ (based upon a) with the boundary-layer thickness δ of order a. Using the equations of mean motion and the method of matched asymptotic expansions, it is shown that the flow can be described by the same two limit processes (inner and outer) as are used in two-dimensional flow. The condition that the two expansions match requires the existence, at the lowest order, of a log region in the usual two-dimensional co-ordinates (u+, y+). Examination of available experimental data shows that substantial log regions do in fact exist but that the intercept is possibly not a universal constant. Similarly, the solution in the outer layer leads to a defect law of the same form as in two-dimensional flow; experiment shows that the intercept in the defect law depends on δ/a. It is concluded that, except in those extreme situations where a+ is small (in which case the boundary layer may not anyway be in a fully developed turbulent state), the simplest analysis of axisymmetric flow will be to use the two-dimensional laws with parameters that now depend on a+ or δ/a as appropriate.

Growth and Characterization of GaAs Epitaxial Layers by MOCVD

technique, on both semi-insulating and semi-conducting CraAs substrates with (100) orientation, offset by 2° towards (110) direction. Systematic variation of As/Ga was performed to gain an understanding of growth process, type of formation and other related physical properties. The films were characterized by using the variety of techniques, such as SEM, EDAX, HRTEM, XRD, and PL. Optical and electrical properties of undoped CyaAs epilayers are presented with reference to the growth conditions and AsH3/TMGa ratio. Photoluminescence measurements of GaAs epilayers were recorded at 4.2K and shows the emission of free exciton and confirmed their high purity. The dominant residual impurities in GaAs are presented by using PL. Finally, electrochemical depth profiling exhibited almost homogeneous background carrier distribution and excellent abruptness between the thin GaAs epilayer and substrate.

Nanopatterning on Reconstructed Ceramic Surfaces

The use of reconstructed ceramic surfaces as templates for nanopatterning has been demonstrated recently. This technique differs from the surface decoration by Au on stepped surfaces of alkali halides which has been a topic of intense research in the past. Some of the intriguing aspects related to the physical origin of the phenomena have been considered here. Based on heterogeneous nucleation of Pt vapor on wedged alumina surfaces, it has been shown that the valley sites are the preferred sites for nucleation. However, the hill sites are decorated by the particles in the present study pointing out to a different physical origin for the formation of the nanoparticles. The role of electrostatic energy reduction on the formation of such nanopatterns is discussed.

Lead Stereochemistry in Incommensurate Ferroelectric Perovskites and in Incommensurate Lead Monoxide

We have performed Rietveld refinements on neutron and synchrotron diffraction patterns and density functional calculations on various ferroelectric lead perovskites and on α lead monoxide (litharge). These structural data have allowed to shed some light on lead stereochemistry in these compounds. In particular, we discuss the changing in the lead behaviour between the paraelectric cubic phases and the low temperature anti or ferroelectric phases in Pb2CoWO6 and Pb2MgTeO6 (both incommensurate), in Pb2MgWO6 (antiferroelectric) and in PbMg1/3Nb2/3O3 (relaxor). The possible phase transition mechanisms are reviewed and the bonds are compared to those in the aperiodic structure of α-lead monoxide.

Investigation of Surface Grooves from Migrating Grain Boundaries

Visible-light microscopy (VLM) and atomic-force microscopy (AFM) were used to study the progression of grain-boundary grooving and migration in high-purity alumina (Lucalox™). Groove profiles from the same grain boundaries were revisited using AFM following successive heat-treatments. The grooves measured from migrating grain boundaries were found to have asymmetric partial-angles compared to those measured from boundaries that did not migrate during the experiment. For a moving boundary, the grain with the larger partial-angle was consistently found to grow into the grain with the smaller partial-angle. Migrating boundaries were observed to leave behind remnant thermal grooves. The observations indicate that the boundary may be bowing out during the migration process.

Ellipsometry as a Sensitive Technique to Probe film -Substrate Interfaces: Al2O3 on Si(100)

An attempt has been made to study the film-substrate interface by using a sensitive, non- conventional tool. Because of the prospective use of gate oxide in MOSFET devices, we have chosen to study alumina films grown on silicon. Film-substrate interface of alumina grown by MOCVD on Si(100) was studied systematically using spectroscopic ellipsometry in the range 1.5-5.0 eV, supported by cross-sectional SEM, and SIMS. The (ε1,ε2) versus energy data obtained for films grown at 600°C, 700°C, and 750°C were modeled to fit a substrate/interface/film “sandwich”. The experimental results reveal (as may be expected) that the nature of the substrate -film interface depends strongly on the growth temperature. The simulated (ε1,ε2) patterns are in excellent agreement with observed ellipsometric data. The MOCVD precursors results the presence of carbon in the films. Theoretical simulation was able to account for the ellipsometry data by invoking the presence of “free” carbon in the alumina films.

Tailoring of Structural Morphology of Silver Nanowires in Electrochemical Growth

Noble metal such as Ag normally exists in an fcc crystal structure. However as the size of the material is decreased to nanometer lengthscales, a structural transformation from that of its bulk state can be expected with new atomic arrangements due to competition between internal packing and minimization of surface energy. In many previous studies, it has been shown that silver nanowires (AGNWs) grown inside anodic alumina (AAO) templates by ac or dc electrochemical deposition from silver salts or complexes, adopt fcc structure and below some critical diameter ∼ 20 nm they may acquire hcp structure at low temperature. This is, however, critically dependant on the nature of confinement, as AgNWs grown inside nanotube confinement with subnanometer diameter have been reported to have fcc structure. Hence the question of the crystal structure of metal nanowires under combined influence of confinement, temperature and deposition condition remains open. In this abstract we show that the alternative crystal structures of AGNWs at room temperature can be achieved with electrochemical growth processes under specific conditions determined by the deposition parameters and nature of confinement. We fabricated AgNWs of 4H hexagonal structure with diameters 30 – 80 nm inside polycarbonate (PC) templates with a modified dc electrodeposition technique, where the nanowires were grown at deposition potentials as low as 10 mV in 2 M silver nitrate solution[1]. We call this low-potential electrodeposition (LPED) since the electrodeposition process occurs at potential much less than the standard Nernst potential (770 mV) of silver. Two types of electrodes were used – stainless steel and sputtered thin Pt film, neither of which had any influence on the crystal structure of the nanowires. EDS elemental analysis showed the nanowires to consist only of silver. Although the precise atomic dynamics during the LPED process is unclear at present, we investigated this with HRTEM (high-resolution transmission electron microscopy) characterization of nanowires grown over various deposition times, as well as electrical conductivity measurements. These experiments indicate that nanowire growth does not occur through a three-dimensional diffusion controlled process, as proposed for conventional over-potential deposition, but follow a novel instantaneous linear growth mechanism. Further experiments showed that, (a) conventional electrochemical growth at a small over-potential in a 2 mM AgNO3 solution yields nanowires with expected fcc structure inside the same PC templates, and (2) no nanowire was observed under the LPED conditions inside hard AAO templates, indicating that LPED-growth process, and hcp structure of the corresponding nanowires depend on deposition parameters, as well as nature of confinement.