Two Modes of Solid State Nucleation - Ferrites, Martensites and Isothermal Transformation Curves

When a crystalline solid such as iron is cooled across a structural transition, its final microstructure depends sensitively on the cooling rate. For instance, an adiabatic cooling across the transition results in an equilibrium `ferrite', while a rapid cooling gives rise to a metastable twinned `martensite'. There exists no theoretical framework to understand the dynamics and conditions under which both these microstructures obtain. Existing theories of martensite dynamics describe this transformation in terms of elastic strain, without any explanation for the occurence of the ferrite. Here we provide evidence for the crucial role played by non-elastic variables, {\it viz.}, dynamically generated interfacial defects. A molecular dynamics (MD) simulation of a model 2-dimensional (2d) solid-state transformation reveals two distinct modes of nucleation depending on the temperature of quench. At high temperatures, defects generated at the nucleation front relax quickly giving rise to an isotropically growing `ferrite'. At low temperatures, the defects relax extremely slowly, forcing a coordinated motion of atoms along specific directions. This results in a twinned critical nucleus which grows rapidly at speeds comparable to that of sound. Based on our MD results, we propose a solid-state nucleation theory involving the elastic strain and non-elastic defects, which successfully describes the transformation to both a ferrite and a martensite. Our work provides useful insights on how to formulate a general dynamics of solid state transformations.Comment: 3 pages, 4 B/W + 2 color figure

Detection of a Series of X-ray Dips Associated with a Radio Flare in GRS 1915+105

We report the detection of a series of X-ray dips in the Galactic black hole candidate GRS 1915+105 during 1999 June 6-17 from observations carried out with the Pointed Proportional Counters of the Indian X-ray Astronomy Experiment on board the Indian satellite IRS-P3. The observations were made after the source made a transition from a steady low-hard state to a chaotic state which occuered within a few hours. Dips of about 20-160 seconds duration are observed on most of the days. The X-ray emission outside the dips shows a QPO at ~ 4 Hz which has characteristics similar to the ubiquitous 0.5 - 10 Hz QPO seen during the low-hard state of the source. During the onset of dips this QPO is absent and also the energy spectrum is soft and the variability is low compared to the non-dip periods. These features gradually re-appear as the dip recovers. The onset of the occurrence of a large number of such dips followed the start of a huge radio flare of strength 0.48 Jy (at 2.25 GHz). We interpret these dips as the cause for mass ejection due to the evacuation of matter from an accretion disk around the black hole. We propose that a super-position of a large number of such dip events produces a huge radio jet in GRS 1915+105.Comment: 18 pages, 7 figures, Accepted for publication in Ap

Effect of nonnegativity on estimation errors in one-qubit state tomography with finite data

We analyze the behavior of estimation errors evaluated by two loss functions, the Hilbert-Schmidt distance and infidelity, in one-qubit state tomography with finite data. We show numerically that there can be a large gap between the estimation errors and those predicted by an asymptotic analysis. The origin of this discrepancy is the existence of the boundary in the state space imposed by the requirement that density matrices be nonnegative (positive semidefinite). We derive an explicit form of a function reproducing the behavior of the estimation errors with high accuracy by introducing two approximations: a Gaussian approximation of the multinomial distributions of outcomes, and linearizing the boundary. This function gives us an intuition for the behavior of the expected losses for finite data sets. We show that this function can be used to determine the amount of data necessary for the estimation to be treated reliably with the asymptotic theory. We give an explicit expression for this amount, which exhibits strong sensitivity to the true quantum state as well as the choice of measurement.Comment: 9 pages, 4 figures, One figure (FIG. 1) is added to the previous version, and some typos are correcte

The Analysis of Large Order Bessel Functions in Gravitational Wave Signals from Pulsars

In this work, we present the analytic treatment of the large order Bessel functions that arise in the Fourier Transform (FT) of the Gravitational Wave (GW) signal from a pulsar. We outline several strategies which employ asymptotic expansions in evaluation of such Bessel functions which also happen to have large argument. Large order Bessel functions also arise in the Peters-Mathews model of binary inspiralling stars emitting GW and several problems in potential scattering theory. Other applications also arise in a variety of problems in Applied Mathematics as well as in the Natural Sciences and present a challenge for High Performance Computing(HPC).Comment: 8 pages, Uses IEEE style files: Ieee.cls, Ieee.clo and floatsty.sty. Accepted for publication in High Performance Computing Symposium, May 15-18 (HPCS 2005) Guelph, Ontario, Canad

Lattice Green's function for crystals containing a planar interface

Flexible boundary condition methods couple an isolated defect to a harmonically responding medium through the bulk lattice Green's function; in the case of an interface, interfacial lattice Green's functions. We present a method to compute the lattice Green's function for a planar interface with arbitrary atomic interactions suited for the study of line defect/interface interactions. The interface is coupled to two different semi-infinite bulk regions, and the Green's function for interface-interface, bulk-interface and bulk-bulk interactions are computed individually. The elastic bicrystal Green's function and the bulk lattice Green's function give the interaction between bulk regions. We make use of partial Fourier transforms to treat in-plane periodicity. Direct inversion of the force constant matrix in the partial Fourier space provides the interface terms. The general method makes no assumptions about the atomic interactions or crystal orientations. We simulate a screw dislocation interacting with a $(10\bar{1}2)$ twin boundary in Ti using flexible boundary conditions and compare with traditional fixed boundary conditions results. Flexible boundary conditions give the correct core structure with significantly less atoms required to relax by energy minimization. This highlights the applicability of flexible boundary conditions methods to modeling defect/interface interactions by \textit{ab initio} methods