Catalog footprint

What is connected

54works
25topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

54 published item(s)

preprint2020arXiv

Age-structured impact of social distancing on the COVID-19 epidemic in India

The outbreak of the novel coronavirus, COVID-19, has been declared a pandemic by the WHO. The structures of social contact critically determine the spread of the infection and, in the absence of vaccines, the control of these structures through large-scale social distancing measures appears to be the most effective means of mitigation. Here we use an age-structured SIR model with social contact matrices obtained from surveys and Bayesian imputation to study the progress of the COVID-19 epidemic in India. The basic reproductive ratio R0 and its time-dependent generalization are computed based on case data, age distribution and social contact structure. The impact of social distancing measures - workplace non-attendance, school closure, lockdown - and their efficacy with durations are then investigated. A three-week lockdown is found insufficient to prevent a resurgence and, instead, protocols of sustained lockdown with periodic relaxation are suggested. Forecasts are provided for the reduction in age-structured morbidity and mortality as a result of these measures. Our study underlines the importance of age and social contact structures in assessing the country-specific impact of mitigatory social distancing.

preprint2020arXiv

Controlled Optofluidic Crystallization of Colloids Tethered at Interfaces

We report experiments that show rapid crystallization of colloids tethered to an oil-water interface in response to laser illumination. This light-induced transition is due to a combination of long-ranged thermophoretic pumping and local optical binding. We show that the flow-induced force on the colloids can be described as the gradient of a potential. The nonequilibrium steady state due to local heating thus admits an effective equilibrium description. The optofluidic manipulation explored in this work opens novel ways to manipulate and assemble colloidal particles

preprint2020arXiv

Inference, prediction and optimization of non-pharmaceutical interventions using compartment models: the PyRoss library

PyRoss is an open-source Python library that offers an integrated platform for inference, prediction and optimisation of NPIs in age- and contact-structured epidemiological compartment models. This report outlines the rationale and functionality of the PyRoss library, with various illustrations and examples focusing on well-mixed, age-structured populations. The PyRoss library supports arbitrary structured models formulated stochastically (as master equations) or deterministically (as ODEs) and allows mid-run transitioning from one to the other. By supporting additional compartmental subdivision ad libitum, PyRoss can emulate time-since-infection models and allows medical stages such as hospitalization or quarantine to be modelled and forecast. The PyRoss library enables fitting to epidemiological data, as available, using Bayesian parameter inference, so that competing models can be weighed by their evidence. PyRoss allows fully Bayesian forecasts of the impact of idealized NPIs by convolving uncertainties arising from epidemiological data, model choice, parameters, and intrinsic stochasticity. Algorithms to optimize time-dependent NPI scenarios against user-defined cost functions are included. PyRoss's current age-structured compartment framework for well-mixed populations will in future reports be extended to include compartments structured by location, occupation, use of travel networks and other attributes relevant to assessing disease spread and the impact of NPIs. We argue that such compartment models, by allowing social data of arbitrary granularity to be combined with Bayesian parameter estimation for poorly-known disease variables, could enable more powerful and robust prediction than other approaches to detailed epidemic modelling. We invite others to use the PyRoss library for research to address today's COVID-19 crisis, and to plan for future pandemics.

preprint2020arXiv

Periodic Orbits of Active Particles induced by Hydrodynamic Monopoles

Terrestrial experiments on active particles, such as Volvox, involve gravitational forces, torques and accompanying monopolar fluid flows. Taking these into account, we analyse the dynamics of a pair of self-propelling, self-spinning active particles between widely separated parallel planes. Neglecting flow reflected by the planes, the dynamics of orientation and horizontal separation is symplectic, with a Hamiltonian exactly determining limit cycle oscillations. Near the bottom plane, gravitational torque damps and reflected flow excites this oscillator, sustaining a second limit cycle that can be perturbatively related to the first. Our work provides a theory for dancing Volvox and highlights the importance of monopolar flow in active matter.

preprint2020arXiv

PyStokes: phoresis and Stokesian hydrodynamics in Python

We present a modular Python library for computing many-body hydrodynamic and phoretic interactions between spherical active particles in suspension, when these are given by solutions of the Stokes and Laplace equations. Underpinning the library is a grid-free methodology that combines dimensionality reduction, spectral expansion, and Ritz-Galerkin discretization, thereby reducing the computation to the solution of a linear system. The system can be solved analytically as a series expansion or numerically at a cost quadratic in the number of particles. Suspension-scale quantities like fluid flow, entropy production, and rheological response are obtained at a small additional cost. The library is agnostic to boundary conditions and includes, amongst others, confinement by plane walls or liquid-liquid interfaces. The use of the library is demonstrated with six fully coded examples simulating active phenomena of current experimental interest.

preprint2020arXiv

Transients generate memory and break hyperbolicity in stochastic enzymatic networks

The hyperbolic dependence of catalytic rate on substrate concentration is a classical result in enzyme kinetics, quantified by the celebrated Michaelis-Menten equation. The ubiquity of this relation in diverse chemical and biological contexts has recently been rationalized by a graph-theoretic analysis of deterministic reaction networks. Experiments, however, have revealed that "molecular noise" - intrinsic stochasticity at the molecular scale - leads to significant deviations from classical results and to unexpected effects like "molecular memory", i.e., the breakdown of statistical independence between turnover events. Here we show, through a new method of analysis, that memory and non-hyperbolicity have a common source in an initial, and observably long, transient peculiar to stochastic reaction networks of multiple enzymes. Networks of single enzymes do not admit such transients. The transient yields, asymptotically, to a steady-state in which memory vanishes and hyperbolicity is recovered. We propose new statistical measures, defined in terms of turnover times, to distinguish between the transient and steady states and apply these to experimental data from a landmark experiment that first observed molecular memory in a single enzyme with multiple binding sites. Our study shows that catalysis at the molecular level with more than one enzyme always contains a non-classical regime and provides insight on how the classical limit is attained.

preprint2018arXiv

Magnetotransport in phase-separated (Ga,Fe)N with $γ$'-Ga$_y$Fe$_{4-y}$N nanocrystals

The magnetotransport in phase-separated (Ga,Fe)N containing $γ$'-Ga$_y$Fe$_{4-y}$N (0\,$<$\,y\,$<$1) nanocrystals (NCs) is studied in the temperature range between 2\,K and 300\,K. The evolution of the resistivity and of the magnetoresistance (MR) as a function of temperature points at two conduction mechanisms: namely a conventional Arrhenius-type one down to 50\,K, and Mott variable range hopping at lower temperatures, where the spin-polarized current is transported between NCs in a regime in which phonon-scattering effects are not dominant. Below 25\,K, the MR shows a hysteretic contribution at magnetic fields $<$1\,T and proportional to the coercive field. Anisotropic magnetoresistance with values one order of magnitude greater than those previously reported for $γ$'-Fe$_4$N thin films over the whole considered temperature range, confirms that the observed MR in these layers is determined by the embedded nanocrystals.

preprint2016arXiv

Rashba semiconductor as spin Hall material: Experimental demonstration of spin pumping in wurtzite $n$-GaN:Si

Pure spin currents in semiconductors are essential for implementation in the next generation of spintronic elements. Heterostructures of III- nitride semiconductors are currently employed as central building-blocks for lighting and high-power devices. Moreover, the long relaxation times and the spin-orbit coupling (SOC) in these materials indicate them as privileged hosts for spin currents and related phenomena. Spin pumping is an efficient mechanism for the inception of spin current and its conversion into charge current in non-magnetic metals and semiconductors with Rashba SOC $via$ spin Hall effects. We report on the generation in $n$-GaN:Si\,--\,at room temperature and through spin pumping\,--\,of pure spin current, fundamental for the understanding of the spin dynamics in these non-centrosymmetric Rashba systems. We find for $n$-GaN:Si a spin Hall angle $θ_{\mathrm{SH}}$=$3.03\times10^{-3}$, exceeding by one order of magnitude those reported for other semiconductors, pointing at III-nitrides as particularly efficient spin current generators.

preprint2016arXiv

Stretching magnetism with an electric field in a nitride semiconductor

By direct magnetization measurements, performed employing a new detection scheme, we demonstrate an electrical control of magnetization in wurtzite (Ga,Mn)N. In this dilute magnetic insulator the Fermi energy is pinned by Mn ions in the mid-gap region, and the Mn3+ ions show strong single-ion anisotropy. We establish that (Ga,Mn)N sustains an electric field up to at least 5 MV/cm, indicating that Mn doping turns GaN into a worthwhile semi-insulating material. Under these conditions, the magnetoelectric coupling may be driven by the inverse piezoelectric effect that stretches the elementary cell along the c axis and, thus, affects the magnitude of magnetic anisotropy. We develop a corresponding theory and show that it describes the experimentally determined dependence of magnetization on the electric field quantitatively with no adjustable parameters as a function of the magnetic field and temperature. In this way, our work bridges two research domains developed so far independently: piezoelectricity of wurtzite semiconductors and electrical control of magnetization in hybrid and composite magnetic structures containing piezoelectric components.

preprint2016arXiv

Universal hydrodynamic mechanisms for crystallization in active colloidal suspensions

The lack of detailed balance in active colloidal suspensions allows dissipation to determine stationary states. Here we show that slow viscous flow produced by polar or apolar active colloids near plane walls mediates attractive hydrodynamic forces that drive crystallization. Hydrodynamically mediated torques tend to destabilize the crystal but stability can be regained through critical amounts of bottom-heaviness or chiral activity. Numerical simulations show that crystallization is not nucleational, as in equilibrium, but is preceded by a spinodal-like instability. Harmonic excitations of the active crystal relax diffusively but the normal modes are distinct from an equilibrium colloidal crystal. The hydrodynamic mechanisms presented here are universal and rationalize recent experiments on the crystallization of active colloids.

preprint2015arXiv

Brownian microhydrodynamics of active filaments

Slender bodies capable of spontaneous motion in the absence of external actuation in an otherwise quiescent fluid are common in biological, physical and technological contexts. The interplay between the spontaneous fluid flow, Brownian motion, and the elasticity of the body presents a challenging fluid-structure interaction problem. Here, we model this problem by approximating the slender body as an elastic filament that can impose non-equilibrium velocities or stresses at the fluid-structure interface. We derive equations of motion for such an active filament by enforcing momentum conservation in the fluid-structure interaction and assuming slow viscous flow in the fluid. The fluid-structure interaction is obtained, to any desired degree of accuracy, through the solution of an integral equation. A simplified form of the equations of motion, that allows for efficient numerical solutions, is obtained by applying the Kirkwood-Riseman superposition approximation to the integral equation. We use this form of the equation of motion to study dynamical steady states in free and hinged minimally active filaments. Our model provides the foundation to study collective phenomena in momentum-conserving, Brownian, active filament suspensions.

preprint2015arXiv

Controlling a three dimensional electron slab of graded Al$_{x}$Ga$_{1-x}$N

Polarization induced degenerate $n$-type doping with electron concentrations up to $\sim$10$^{20}$\,cm$^{-3}$ is achieved in graded Al$_{x}$Ga$_{1-x}$N layers ($x$: 0\%$\rightarrow$37\%) grown on unintentionally doped and on $n$-doped GaN:Si buffer/reservoir layers by metal organic vapor phase epitaxy. High resolution x-ray diffraction, transmission electron microscopy and electron energy loss spectroscopy confirm the gradient in the composition of the Al$_{x}$Ga$_{1-x}$N layers, while magnetotransport studies reveal the formation of a three dimensional electron slab, whose conductivity can be adjusted through the GaN(:Si) buffer/reservoir.

preprint2015arXiv

Flow-induced nonequilibrium self-assembly in suspensions of stiff, apolar, active filaments

Active bodies in viscous fluids interact hydrodynamically through self-generated flows. Here we study spontaneous aggregation induced by hydrodynamic flow in a suspension of stiff, apolar, active filaments. Lateral hydrodynamic attractions in extensile filaments lead, independent of volume fraction, to anisotropic aggregates which translate and rotate ballistically. Lateral hydrodynamic repulsion in contractile filaments lead, with increasing volume fractions, to microstructured states of asters, clusters, and incipient gels where, in each case, filament motion is diffusive. Our results demonstrate that the interplay of active hydrodynamic flows and anisotropic excluded volume interactions provides a generic nonequilibrium mechanism for hierarchical self-assembly of active soft matter.

preprint2015arXiv

Many-body microhydrodynamics of colloidal particles with active boundary layers

Colloidal particles with active boundary layers - regions surrounding the particles where nonequilibrium processes produce large velocity gradients - are common in many physical, chemical and biological contexts. The velocity or stress at the edge of the boundary layer determines the exterior fluid flow and, hence, the many-body interparticle hydrodynamic interaction. Here, we present a method to compute the many-body hydrodynamic interaction between $N$ spherical active particles induced by their exterior microhydrodynamic flow. First, we use a boundary integral representation of the Stokes equation to eliminate bulk fluid degrees of freedom. Then, we expand the boundary velocities and tractions of the integral representation in an infinite-dimensional basis of tensorial spherical harmonics and, on enforcing boundary conditions in a weak sense on the surface of each particle, obtain a system of linear algebraic equations for the unknown expansion coefficients. The truncation of the infinite series, fixed by the degree of accuracy required, yields a finite linear system that can be solved accurately and efficiently by iterative methods. The solution linearly relates the unknown rigid body motion to the known values of the expansion coefficients, motivating the introduction of propulsion matrices. These matrices completely characterize hydrodynamic interactions in active suspensions just as mobility matrices completely characterize hydrodynamic interactions in passive suspensions. The reduction in the dimensionality of the problem, from a three-dimensional partial differential equation to a two-dimensional integral equation, allows for dynamic simulations of hundreds of thousands of active particles on multi-core computational architectures.

preprint2014arXiv

Experimental determination of Rashba spin-orbit coupling in wurtzite $n$-GaN:Si

Millikelvin magnetotransport studies are carried out on heavily $n$-doped wurtzite GaN:Si films grown on semi-insulating GaN:Mn buffer layers by metal-organic vapor phase epitaxy. The dependency of the conductivity on magnetic field and temperature is interpreted in terms of theories that take into account disorder-induced quantum interference of one-electron and many-electron self-crossing trajectories. The Rashba parameter $α_{\text{R}}\,=\,(4.5 \pm 1)$ meV$Å$ is determined, and it is shown that in the previous studies of electrons adjacent to GaN/(Al,Ga)N interfaces, bulk inversion asymmetry was dominant over structural inversion asymmetry. The comparison of experimental and theoretical values of $α_{\text{R}}$ across a series of wurtzite semiconductors is presented as a test of current relativistic ab initio computation schemes. It is found that electron-electron scattering with small energy transfer accounts for low temperature decoherence in these systems.

preprint2014arXiv

Irreducible Representations Of Oscillatory And Swirling Flows In Active Soft Matter

Recent experiments imaging fluid flow around swimming microorganisms have revealed complex time-dependent velocity fields that differ qualitatively from the stresslet flow commonly employed in theoretical descriptions of active matter. Here we obtain the most general flow around a finite sized active particle by expanding the surface stress in irreducible Cartesian tensors. This expansion, whose first term is the stresslet, must include, respectively, third-rank polar and axial tensors to minimally capture crucial features of the active oscillatory flow around translating Chlamydomonas and the active swirling flow around rotating Volvox. The representation provides explicit expressions for the irreducible symmetric, antisymmetric and isotropic parts of the continuum active stress. Antisymmetric active stresses do not conserve orbital angular momentum and our work thus shows that spin angular momentum is necessary to restore angular momentum conservation in continuum hydrodynamic descriptions of active soft matter.

preprint2014arXiv

Upper bound for the s-d exchange integral in n-(Ga,Mn)N:Si from magnetotransport studies

A series of recent magnetooptical studies pointed to contradicting values of the s-d exchange energy N0α in Mn-doped GaAs and GaN as well as in Fe-doped GaN. Here, a strong sensitivity of weak-localization phenomena to symmetry breaking perturbations (such as spin-splitting and spin-disorder scattering) is exploited to evaluate the magnitude of N0α for n-type wurtzite (Ga,Mn)N:Si films grown by metalorganic vapor phase epitaxy. Millikelvin magnetoresistance studies and their quantitative interpretation point to N0α < 40 meV, a value at least 5 times smaller than the one found with similar measurements on, e.g., $n$-(Zn,Mn)O. It is shown that this striking difference in the values of the s-d coupling between $n$-type III-V and II-VI dilute magnetic semiconductors can be explained by a theory that takes into account the acceptor character of Mn in III-V compounds.

preprint2013arXiv

A First Search for coincident Gravitational Waves and High Energy Neutrinos using LIGO, Virgo and ANTARES data from 2007

We present the results of the first search for gravitational wave bursts associated with high energy neutrinos. Together, these messengers could reveal new, hidden sources that are not observed by conventional photon astronomy, particularly at high energy. Our search uses neutrinos detected by the underwater neutrino telescope ANTARES in its 5 line configuration during the period January - September 2007, which coincided with the fifth and first science runs of LIGO and Virgo, respectively. The LIGO-Virgo data were analysed for candidate gravitational-wave signals coincident in time and direction with the neutrino events. No significant coincident events were observed. We place limits on the density of joint high energy neutrino - gravitational wave emission events in the local universe, and compare them with densities of merger and core-collapse events.

preprint2013arXiv

Hydrodynamic Instabilities Provide A Generic Route To Spontaneous Biomimetic Oscillations In Chemomechanically Active Filaments

Non-equilibrium processes which convert chemical energy into mechanical motion enable the motility of organisms. Bundles of inextensible filaments driven by energy transduction of molecular motors form essential components of micron-scale motility engines like cilia and flagella. The mimicry of cilia-like motion in recent experiments on synthetic active filaments supports the idea that generic physical mechanisms may be sufficient to generate such motion. Here we show, theoretically, that the competition between the destabilising effect of hydrodynamic interactions induced by force-free and torque-free chemomechanically active flows, and the stabilising effect of nonlinear elasticity, provides a generic route to spontaneous oscillations in active filaments. These oscillations, reminiscent of prokaryotic and eukaryotic flagellar motion, are obtained without having to invoke structural complexity or biochemical regulation. This minimality implies that biomimetic oscillations, previously observed only in complex bundles of active filaments, can be replicated in simple chains of generic chemomechanically active beads.

preprint2013arXiv

Parameter estimation for compact binary coalescence signals with the first generation gravitational-wave detector network

Compact binary systems with neutron stars or black holes are one of the most promising sources for ground-based gravitational wave detectors. Gravitational radiation encodes rich information about source physics; thus parameter estimation and model selection are crucial analysis steps for any detection candidate events. Detailed models of the anticipated waveforms enable inference on several parameters, such as component masses, spins, sky location and distance that are essential for new astrophysical studies of these sources. However, accurate measurements of these parameters and discrimination of models describing the underlying physics are complicated by artifacts in the data, uncertainties in the waveform models and in the calibration of the detectors. Here we report such measurements on a selection of simulated signals added either in hardware or software to the data collected by the two LIGO instruments and the Virgo detector during their most recent joint science run, including a "blind injection" where the signal was not initially revealed to the collaboration. We exemplify the ability to extract information about the source physics on signals that cover the neutron star and black hole parameter space over the individual mass range 1 Msun - 25 Msun and the full range of spin parameters. The cases reported in this study provide a snap-shot of the status of parameter estimation in preparation for the operation of advanced detectors.

preprint2013arXiv

Search for Gravitational Waves from Binary Black Hole Inspiral, Merger and Ringdown in LIGO-Virgo Data from 2009-2010

We report a search for gravitational waves from the inspiral, merger and ringdown of binary black holes (BBH) with total mass between 25 and 100 solar masses, in data taken at the LIGO and Virgo observatories between July 7, 2009 and October 20, 2010. The maximum sensitive distance of the detectors over this period for a (20,20) Msun coalescence was 300 Mpc. No gravitational wave signals were found. We thus report upper limits on the astrophysical coalescence rates of BBH as a function of the component masses for non-spinning components, and also evaluate the dependence of the search sensitivity on component spins aligned with the orbital angular momentum. We find an upper limit at 90% confidence on the coalescence rate of BBH with non-spinning components of mass between 19 and 28 Msun of 3.3 \times 10^-7 mergers /Mpc^3 /yr.

preprint2012arXiv

All-sky search for gravitational-wave bursts in the second joint LIGO-Virgo run

We present results from a search for gravitational-wave bursts in the data collected by the LIGO and Virgo detectors between July 7, 2009 and October 20, 2010: data are analyzed when at least two of the three LIGO-Virgo detectors are in coincident operation, with a total observation time of 207 days. The analysis searches for transients of duration < 1 s over the frequency band 64-5000 Hz, without other assumptions on the signal waveform, polarization, direction or occurrence time. All identified events are consistent with the expected accidental background. We set frequentist upper limits on the rate of gravitational-wave bursts by combining this search with the previous LIGO-Virgo search on the data collected between November 2005 and October 2007. The upper limit on the rate of strong gravitational-wave bursts at the Earth is 1.3 events per year at 90% confidence. We also present upper limits on source rate density per year and Mpc^3 for sample populations of standard-candle sources. As in the previous joint run, typical sensitivities of the search in terms of the root-sum-squared strain amplitude for these waveforms lie in the range 5 10^-22 Hz^-1/2 to 1 10^-20 Hz^-1/2. The combination of the two joint runs entails the most sensitive all-sky search for generic gravitational-wave bursts and synthesizes the results achieved by the initial generation of interferometric detectors.

preprint2012arXiv

Autonomous Motility of Active Filaments due to Spontaneous Flow-Symmetry Breaking

We simulate the nonlocal Stokesian hydrodynamics of an elastic filament which is active due a permanent distribution of stresslets along its contour. A bending instability of an initially straight filament spontaneously breaks flow symmetry and leads to autonomous filament motion which, depending on conformational symmetry, can be translational or rotational. At high ratios of activity to elasticity, the linear instability develops into nonlinear fluctuating states with large amplitude deformations. The dynamics of these states can be qualitatively understood as a superposition of translational and rotational motion associated with filament conformational modes of opposite symmetry. Our results can be tested in molecular-motor filament mixtures, synthetic chains of autocatalytic particles, or other linearly connected systems where chemical energy is converted to mechanical energy in a fluid environment.

preprint2012arXiv

Einstein@Home all-sky search for periodic gravitational waves in LIGO S5 data

This paper presents results of an all-sky searches for periodic gravitational waves in the frequency range [50, 1190] Hz and with frequency derivative ranges of [-2 x 10^-9, 1.1 x 10^-10] Hz/s for the fifth LIGO science run (S5). The novelty of the search lies in the use of a non-coherent technique based on the Hough-transform to combine the information from coherent searches on timescales of about one day. Because these searches are very computationally intensive, they have been deployed on the Einstein@Home distributed computing project infrastructure. The search presented here is about a factor 3 more sensitive than the previous Einstein@Home search in early S5 LIGO data. The post-processing has left us with eight surviving candidates. We show that deeper follow-up studies rule each of them out. Hence, since no statistically significant gravitational wave signals have been detected, we report upper limits on the intrinsic gravitational wave amplitude h0. For example, in the 0.5 Hz-wide band at 152.5 Hz, we can exclude the presence of signals with h0 greater than 7.6 x 10^-25 with a 90% confidence level.

preprint2012arXiv

Implementation and testing of the first prompt search for gravitational wave transients with electromagnetic counterparts

Aims. A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize GW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify GW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency GW pipeline's ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with ~50% or better probability with a few pointings of wide-field telescopes.

preprint2012arXiv

Implications For The Origin Of GRB 051103 From LIGO Observations

We present the results of a LIGO search for gravitational waves (GWs) associated with GRB 051103, a short-duration hard-spectrum gamma-ray burst (GRB) whose electromagnetically determined sky position is coincident with the spiral galaxy M81, which is 3.6 Mpc from Earth. Possible progenitors for short-hard GRBs include compact object mergers and soft gamma repeater (SGR) giant flares. A merger progenitor would produce a characteristic GW signal that should be detectable at the distance of M81, while GW emission from an SGR is not expected to be detectable at that distance. We found no evidence of a GW signal associated with GRB 051103. Assuming weakly beamed gamma-ray emission with a jet semi-angle of 30 deg we exclude a binary neutron star merger in M81 as the progenitor with a confidence of 98%. Neutron star-black hole mergers are excluded with > 99% confidence. If the event occurred in M81 our findings support the the hypothesis that GRB 051103 was due to an SGR giant flare, making it the most distant extragalactic magnetar observed to date.

preprint2012arXiv

Search for gravitational waves associated with gamma-ray bursts during LIGO science run 6 and Virgo science runs 2 and 3

We present the results of a search for gravitational waves associated with 154 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments in 2009-2010, during the sixth LIGO science run and the second and third Virgo science runs. We perform two distinct searches: a modeled search for coalescences of either two neutron stars or a neutron star and black hole; and a search for generic, unmodeled gravitational-wave bursts. We find no evidence for gravitational-wave counterparts, either with any individual GRB in this sample or with the population as a whole. For all GRBs we place lower bounds on the distance to the progenitor, under the optimistic assumption of a gravitational-wave emission energy of 10^-2 M c^2 at 150 Hz, with a median limit of 17 Mpc. For short hard GRBs we place exclusion distances on binary neutron star and neutron star-black hole progenitors, using astrophysically motivated priors on the source parameters, with median values of 16 Mpc and 28 Mpc respectively. These distance limits, while significantly larger than for a search that is not aided by GRB satellite observations, are not large enough to expect a coincidence with a GRB. However, projecting these exclusions to the sensitivities of Advanced LIGO and Virgo, which should begin operation in 2015, we find that the detection of gravitational waves associated with GRBs will become quite possible.

preprint2012arXiv

Search for Gravitational Waves from Intermediate Mass Binary Black Holes

We present the results of a weakly modeled burst search for gravitational waves from mergers of non-spinning intermediate mass black holes (IMBH) in the total mass range 100--450 solar masses and with the component mass ratios between 1:1 and 4:1. The search was conducted on data collected by the LIGO and Virgo detectors between November of 2005 and October of 2007. No plausible signals were observed by the search which constrains the astrophysical rates of the IMBH mergers as a function of the component masses. In the most efficiently detected bin centered on 88+88 solar masses, for non-spinning sources, the rate density upper limit is 0.13 per Mpc^3 per Myr at the 90% confidence level.

preprint2012arXiv

Search for Gravitational Waves from Low Mass Compact Binary Coalescence in LIGO's Sixth Science Run and Virgo's Science Runs 2 and 3

We report on a search for gravitational waves from coalescing compact binaries using LIGO and Virgo observations between July 7, 2009 and October 20, 2010. We searched for signals from binaries with total mass between 2 and 25 solar masses; this includes binary neutron stars, binary black holes, and binaries consisting of a black hole and neutron star. The detectors were sensitive to systems up to 40 Mpc distant for binary neutron stars, and further for higher mass systems. No gravitational-wave signals were detected. We report upper limits on the rate of compact binary coalescence as a function of total mass, including the results from previous LIGO and Virgo observations. The cumulative 90%-confidence rate upper limits of the binary coalescence of binary neutron star, neutron star- black hole and binary black hole systems are 1.3 x 10^{-4}, 3.1 x 10^{-5} and 6.4 x 10^{-6} Mpc^{-3}yr^{-1}, respectively. These upper limits are up to a factor 1.4 lower than previously derived limits. We also report on results from a blind injection challenge.

preprint2012arXiv

Swift follow-up observations of candidate gravitational-wave transient events

We present the first multi-wavelength follow-up observations of two candidate gravitational-wave (GW) transient events recorded by LIGO and Virgo in their 2009-2010 science run. The events were selected with low latency by the network of GW detectors and their candidate sky locations were observed by the Swift observatory. Image transient detection was used to analyze the collected electromagnetic data, which were found to be consistent with background. Off-line analysis of the GW data alone has also established that the selected GW events show no evidence of an astrophysical origin; one of them is consistent with background and the other one was a test, part of a "blind injection challenge". With this work we demonstrate the feasibility of rapid follow-ups of GW transients and establish the sensitivity improvement joint electromagnetic and GW observations could bring. This is a first step toward an electromagnetic follow-up program in the regime of routine detections with the advanced GW instruments expected within this decade. In that regime multi-wavelength observations will play a significant role in completing the astrophysical identification of GW sources. We present the methods and results from this first combined analysis and discuss its implications in terms of sensitivity for the present and future instruments.

preprint2012arXiv

The characterization of Virgo data and its impact on gravitational-wave searches

Between 2007 and 2010 Virgo collected data in coincidence with the LIGO and GEO gravitational-wave (GW) detectors. These data have been searched for GWs emitted by cataclysmic phenomena in the universe, by non-axisymmetric rotating neutron stars or from a stochastic background in the frequency band of the detectors. The sensitivity of GW searches is limited by noise produced by the detector or its environment. It is therefore crucial to characterize the various noise sources in a GW detector. This paper reviews the Virgo detector noise sources, noise propagation, and conversion mechanisms which were identified in the three first Virgo observing runs. In many cases, these investigations allowed us to mitigate noise sources in the detector, or to selectively flag noise events and discard them from the data. We present examples from the joint LIGO-GEO-Virgo GW searches to show how well noise transients and narrow spectral lines have been identified and excluded from the Virgo data. We also discuss how detector characterization can improve the astrophysical reach of gravitational-wave searches.

preprint2012arXiv

Upper limits on a stochastic gravitational-wave background using LIGO and Virgo interferometers at 600-1000 Hz

A stochastic background of gravitational waves is expected to arise from a superposition of many incoherent sources of gravitational waves, of either cosmological or astrophysical origin. This background is a target for the current generation of ground-based detectors. In this article we present the first joint search for a stochastic background using data from the LIGO and Virgo interferometers. In a frequency band of 600-1000 Hz, we obtained a 95% upper limit on the amplitude of $Ω_{\rm GW}(f) = Ω_3 (f/900 \mathrm{Hz})^3$, of $Ω_3 < 0.33$, assuming a value of the Hubble parameter of $h_{100}=0.72$. These new limits are a factor of seven better than the previous best in this frequency band.

preprint2011arXiv

All-sky Search for Periodic Gravitational Waves in the Full S5 LIGO Data

We report on an all-sky search for periodic gravitational waves in the frequency band 50-800 Hz and with the frequency time derivative in the range of 0 through -6e-9 Hz/s. Such a signal could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. After recent improvements in the search program that yielded a 10x increase in computational efficiency, we have searched in two years of data collected during LIGO's fifth science run and have obtained the most sensitive all-sky upper limits on gravitational wave strain to date. Near 150 Hz our upper limit on worst-case linearly polarized strain amplitude $h_0$ is 1e-24, while at the high end of our frequency range we achieve a worst-case upper limit of 3.8e-24 for all polarizations and sky locations. These results constitute a factor of two improvement upon previously published data. A new detection pipeline utilizing a Loosely Coherent algorithm was able to follow up weaker outliers, increasing the volume of space where signals can be detected by a factor of 10, but has not revealed any gravitational wave signals. The pipeline has been tested for robustness with respect to deviations from the model of an isolated neutron star, such as caused by a low-mass or long-period binary companion.

preprint2011arXiv

Beating the spin-down limit on gravitational wave emission from the Vela pulsar

We present direct upper limits on continuous gravitational wave emission from the Vela pulsar using data from the Virgo detector's second science run. These upper limits have been obtained using three independent methods that assume the gravitational wave emission follows the radio timing. Two of the methods produce frequentist upper limits for an assumed known orientation of the star's spin axis and value of the wave polarization angle of, respectively, $1.9\ee{-24}$ and $2.2\ee{-24}$, with 95% confidence. The third method, under the same hypothesis, produces a Bayesian upper limit of $2.1\ee{-24}$, with 95% degree of belief. These limits are below the indirect {\it spin-down limit} of $3.3\ee{-24}$ for the Vela pulsar, defined by the energy loss rate inferred from observed decrease in Vela's spin frequency, and correspond to a limit on the star ellipticity of $\sim 10^{-3}$. Slightly less stringent results, but still well below the spin-down limit, are obtained assuming the star's spin axis inclination and the wave polarization angles are unknown.

preprint2011arXiv

Effects of Mirror Aberrations on Laguerre-Gaussian Beams in Interferometric Gravitational-Wave Detectors

A fundamental limit to the sensitivity of optical interferometers is imposed by Brownian thermal fluctuations of the mirrors' surfaces. This thermal noise can be reduced by using larger beams which "average out" the random fluctuations of the surfaces. It has been proposed previously that wider, higher-order Laguerre-Gaussian modes can be used to exploit this effect. In this article, we show that susceptibility to spatial imperfections of the mirrors' surfaces limits the effectiveness of this approach in interferometers used for gravitational-wave detection. Possible methods of reducing this susceptibility are also discussed.

preprint2011arXiv

Langevin theory of fluctuations in the discrete Boltzmann equation

The discrete Boltzmann equation for both the ideal and a non-ideal fluid is extended by adding Langevin noise terms in order to incorporate the effects of thermal fluctuations. After casting the fluctuating discrete Boltzmann equation in a form appropriate to the Onsager-Machlup theory of linear fluctuations, the statistical properties of the noise are determined by invoking a fluctuation-dissipation theorem at the kinetic level. By integrating the fluctuating discrete Boltzmann equation, the fluctuating lattice Boltzmann equation is obtained, which provides an efficient way to solve the equations of fluctuating hydrodynamics for ideal and non-ideal fluids. Application of the framework to a generic force-based non-ideal fluid model leads to ideal gas-type thermal noise. Simulation results indicate proper thermalization of all degrees of freedom.

preprint2011arXiv

Non-renewal statistics in the catalytic activity of enzyme molecules at mesoscopic concentrations

Recent fluorescence spectroscopy measurements of single-enzyme kinetics have shown that enzymatic turnovers form a renewal stochastic process in which the inverse of the mean waiting time between turnovers follows the Michaelis-Menten equation. Under typical physiological conditions, however, tens to thousands of enzymes react in catalyzing thousands to millions of substrates. We study enzyme kinetics at these physiologically relevant conditions through a master equation including stochasticity and molecular discreteness. From the exact solution of the master equation we find that the waiting times are neither independent nor are they identically distributed, implying that enzymatic turnovers form a non-renewal stochastic process. The inverse of the mean waiting time shows strong departures from the Michaelis-Menten equation. The waiting times between consecutive turnovers are anti-correlated, where short intervals are more likely to be followed by long intervals and vice versa. Correlations persist beyond consecutive turnovers indicating that multi-scale fluctuations govern enzyme kinetics.

preprint2011arXiv

Search for gravitational waves from binary black hole inspiral, merger and ringdown

We present the first modeled search for gravitational waves using the complete binary black hole gravitational waveform from inspiral through the merger and ringdown for binaries with negligible component spin. We searched approximately 2 years of LIGO data taken between November 2005 and September 2007 for systems with component masses of 1-99 solar masses and total masses of 25-100 solar masses. We did not detect any plausible gravitational-wave signals but we do place upper limits on the merger rate of binary black holes as a function of the component masses in this range. We constrain the rate of mergers for binary black hole systems with component masses between 19 and 28 solar masses and negligible spin to be no more than 2.0 per Mpc^3 per Myr at 90% confidence.

preprint2011arXiv

Technique for high axial shielding factor performance of large-scale, thin, open-ended, cylindrical Metglas magnetic shields

Metglas 2705M is a low-cost commercially-available, high-permeability Cobalt-based magnetic alloy, provided as a 5.08-cm wide and 20.3-$μ$m thick ribbon foil. We present an optimized construction technique for single-shell, large-scale (human-size), thin, open-ended cylindrical Metglas magnetic shields. The measured DC axial and transverse magnetic shielding factors of our 0.61-m diameter and 1.83-m long shields in the Earth's magnetic field were 267 and 1500, for material thicknesses of only 122 $μ$m (i.e., 6 foil layers). The axial shielding performance of our single-shell Metglas magnetic shields, obtained without the use of magnetic shaking techniques, is comparable to the performance of significantly thicker, multiple-shell, open-ended Metglas magnetic shields in comparable-magnitude, low-frequency applied external fields reported previously in the literature.

preprint2010arXiv

A search for gravitational waves associated with the August 2006 timing glitch of the Vela pulsar

The physical mechanisms responsible for pulsar timing glitches are thought to excite quasi-normal mode oscillations in their parent neutron star that couple to gravitational wave emission. In August 2006, a timing glitch was observed in the radio emission of PSR B0833-45, the Vela pulsar. At the time of the glitch, the two co-located Hanford gravitational wave detectors of the Laser Interferometer Gravitational-wave observatory (LIGO) were operational and taking data as part of the fifth LIGO science run (S5). We present the first direct search for the gravitational wave emission associated with oscillations of the fundamental quadrupole mode excited by a pulsar timing glitch. No gravitational wave detection candidate was found. We place Bayesian 90% confidence upper limits of 6.3e-21 to 1.4e-20 on the peak intrinsic strain amplitude of gravitational wave ring-down signals, depending on which spherical harmonic mode is excited. The corresponding range of energy upper limits is 5.0e44 to 1.3e45 erg.

preprint2010arXiv

Calibration of the LIGO Gravitational Wave Detectors in the Fifth Science Run

The Laser Interferometer Gravitational Wave Observatory (LIGO) is a network of three detectors built to detect local perturbations in the space-time metric from astrophysical sources. These detectors, two in Hanford, WA and one in Livingston, LA, are power-recycled Fabry-Perot Michelson interferometers. In their fifth science run (S5), between November 2005 and October 2007, these detectors accumulated one year of triple coincident data while operating at their designed sensitivity. In this paper, we describe the calibration of the instruments in the S5 data set, including measurement techniques and uncertainty estimation.

preprint2010arXiv

Endogenous Quasicycles and Stochastic Coherence in a Closed Endemic Model

We study the role of demographic fluctuations in typical endemics as exemplified by the stochastic SIRS model. The birth-death master equation of the model is simulated using exact numerics and analysed within the linear noise approximation. The endemic fixed point is unstable to internal demographic noise, and leads to sustained oscillations. This is ensured when the eigenvalues ($λ$) of the linearised drift matrix are complex, which in turn, is possible only if detailed balance is violated. In the oscillatory state, the phases decorrelate asymptotically, distinguishing such oscillations from those produced by external periodic forcing. These so-called quasicycles are of sufficient strength to be detected reliably only when the ratio $|Im(λ)/Re(λ)|$ is of order unity. The coherence or regularity of these oscillations show a maximum as a function of population size, an effect known variously as stochastic coherence or coherence resonance. We find that stochastic coherence can be simply understood as resulting from a non-monotonic variation of $|Im(λ)/Re(λ)|$ with population size. Thus, within the linear noise approximation, stochastic coherence can be predicted from a purely deterministic analysis. The non-normality of the linearised drift matrix, associated with the violation of detailed balance, leads to enhanced fluctuations in the population amplitudes.

preprint2010arXiv

First search for gravitational waves from the youngest known neutron star

We present a search for periodic gravitational waves from the neutron star in the supernova remnant Cassiopeia A. The search coherently analyzes data in a 12-day interval taken from the fifth science run of the Laser Interferometer Gravitational-Wave Observatory. It searches gravitational wave frequencies from 100 to 300 Hz, and covers a wide range of first and second frequency derivatives appropriate for the age of the remnant and for different spin-down mechanisms. No gravitational wave signal was detected. Within the range of search frequencies, we set 95% confidence upper limits of 0.7--1.2e-24 on the intrinsic gravitational wave strain, 0.4--4e-4 on the equatorial ellipticity of the neutron star, and 0.005--0.14 on the amplitude of r-mode oscillations of the neutron star. These direct upper limits beat indirect limits derived from energy conservation and enter the range of theoretical predictions involving crystalline exotic matter or runaway r-modes. This is the first gravitational wave search to present upper limits on r-modes.

preprint2010arXiv

Fluctuating Nematodynamics using the Stochastic Method of Lines

We construct Langevin equations describing the fluctuations of the tensor order parameter $Q_{αβ}$ in nematic liquid crystals by adding noise terms to time-dependent variational equations that follow from the Ginzburg-Landau-de Gennes free energy. The noise is required to preserve the symmetry and tracelessness of the tensor order parameter and must satisfy a fluctuation-dissipation relation at thermal equilibrium. We construct a noise with these properties in a basis of symmetric traceless matrices and show that the Langevin equations can be solved numerically in this basis using a stochastic version of the method of lines. The numerical method is validated by comparing equilibrium probability distributions, structure factors and dynamic correlations obtained from these numerical solutions with analytic predictions. We demonstrate excellent agreement between numerics and theory. This methodology can be applied to the study of phenomena where fluctuations in both the magnitude and direction of nematic order are important, as for instance in the nematic swarms which produce enhanced opalescence near the isotropic-nematic transition or the problem of nucleation of the nematic from the isotropic phase.

preprint2010arXiv

Predictions for the Rates of Compact Binary Coalescences Observable by Ground-based Gravitational-wave Detectors

We present an up-to-date, comprehensive summary of the rates for all types of compact binary coalescence sources detectable by the Initial and Advanced versions of the ground-based gravitational-wave detectors LIGO and Virgo. Astrophysical estimates for compact-binary coalescence rates depend on a number of assumptions and unknown model parameters, and are still uncertain. The most confident among these estimates are the rate predictions for coalescing binary neutron stars which are based on extrapolations from observed binary pulsars in our Galaxy. These yield a likely coalescence rate of 100 per Myr per Milky Way Equivalent Galaxy (MWEG), although the rate could plausibly range from 1 per Myr per MWEG to 1000 per Myr per MWEG. We convert coalescence rates into detection rates based on data from the LIGO S5 and Virgo VSR2 science runs and projected sensitivities for our Advanced detectors. Using the detector sensitivities derived from these data, we find a likely detection rate of 0.02 per year for Initial LIGO-Virgo interferometers, with a plausible range between 0.0002 and 0.2 per year. The likely binary neutron-star detection rate for the Advanced LIGO-Virgo network increases to 40 events per year, with a range between 0.4 and 400 per year.

preprint2010arXiv

Run-and-tumble particles with hydrodynamics: sedimentation, trapping and upstream swimming

We simulate by lattice Boltzmann the nonequilibrium steady states of run-and-tumble particles (inspired by a minimal model of bacteria), interacting by far-field hydrodynamics, subject to confinement. Under gravity, hydrodynamic interactions barely perturb the steady state found without them, but for particles in a harmonic trap such a state is quite changed if the run length is larger than the confinement length: a self-assembled pump is formed. Particles likewise confined in a narrow channel show a generic upstream flux in Poiseuille flow: chiral swimming is not required.

preprint2010arXiv

Search for gravitational-wave bursts associated with gamma-ray bursts using data from LIGO Science Run 5 and Virgo Science Run 1

We present the results of a search for gravitational-wave bursts associated with 137 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments during the fifth LIGO science run and first Virgo science run. The data used in this analysis were collected from 2005 November 4 to 2007 October 1, and most of the GRB triggers were from the Swift satellite. The search uses a coherent network analysis method that takes into account the different locations and orientations of the interferometers at the three LIGO-Virgo sites. We find no evidence for gravitational-wave burst signals associated with this sample of GRBs. Using simulated short-duration (<1 s) waveforms, we set upper limits on the amplitude of gravitational waves associated with each GRB. We also place lower bounds on the distance to each GRB under the assumption of a fixed energy emission in gravitational waves, with typical limits of D ~ 15 Mpc (E_GW^iso / 0.01 M_o c^2)^1/2 for emission at frequencies around 150 Hz, where the LIGO-Virgo detector network has best sensitivity. We present astrophysical interpretations and implications of these results, and prospects for corresponding searches during future LIGO-Virgo runs.

preprint2010arXiv

Search for gravitational-wave inspiral signals associated with short Gamma-Ray Bursts during LIGO's fifth and Virgo's first science run

Progenitor scenarios for short gamma-ray bursts (short GRBs) include coalescenses of two neutron stars or a neutron star and black hole, which would necessarily be accompanied by the emission of strong gravitational waves. We present a search for these known gravitational-wave signatures in temporal and directional coincidence with 22 GRBs that had sufficient gravitational-wave data available in multiple instruments during LIGO's fifth science run, S5, and Virgo's first science run, VSR1. We find no statistically significant gravitational-wave candidates within a [-5, +1) s window around the trigger time of any GRB. Using the Wilcoxon-Mann-Whitney U test, we find no evidence for an excess of weak gravitational-wave signals in our sample of GRBs. We exclude neutron star-black hole progenitors to a median 90% CL exclusion distance of 6.7 Mpc.

preprint2010arXiv

Searches for gravitational waves from known pulsars with S5 LIGO data

We present a search for gravitational waves from 116 known millisecond and young pulsars using data from the fifth science run of the LIGO detectors. For this search ephemerides overlapping the run period were obtained for all pulsars using radio and X-ray observations. We demonstrate an updated search method that allows for small uncertainties in the pulsar phase parameters to be included in the search. We report no signal detection from any of the targets and therefore interpret our results as upper limits on the gravitational wave signal strength. The most interesting limits are those for young pulsars. We present updated limits on gravitational radiation from the Crab pulsar, where the measured limit is now a factor of seven below the spin-down limit. This limits the power radiated via gravitational waves to be less than ~2% of the available spin-down power. For the X-ray pulsar J0537-6910 we reach the spin-down limit under the assumption that any gravitational wave signal from it stays phase locked to the X-ray pulses over timing glitches, and for pulsars J1913+1011 and J1952+3252 we are only a factor of a few above the spin-down limit. Of the recycled millisecond pulsars several of the measured upper limits are only about an order of magnitude above their spin-down limits. For these our best (lowest) upper limit on gravitational wave amplitude is 2.3x10^-26 for J1603-7202 and our best (lowest) limit on the inferred pulsar ellipticity is 7.0x10^-8 for J2124-3358.

preprint2010arXiv

Sensitivity to Gravitational Waves from Compact Binary Coalescences Achieved during LIGO's Fifth and Virgo's First Science Run

We summarize the sensitivity achieved by the LIGO and Virgo gravitational wave detectors for compact binary coalescence (CBC) searches during LIGO's fifth science run and Virgo's first science run. We present noise spectral density curves for each of the four detectors that operated during these science runs which are representative of the typical performance achieved by the detectors for CBC searches. These spectra are intended for release to the public as a summary of detector performance for CBC searches during these science runs.

preprint2010arXiv

Thermal fluctuations in the lattice Boltzmann method for non-ideal fluids

We introduce thermal fluctuations in the lattice Boltzmann method for non-ideal fluids. A fluctuation-dissipation theorem is derived within the Langevin framework and applied to a specific lattice Boltzmann model that approximates the linearized fluctuating Navier-Stokes equations for fluids based on square-gradient free energy functionals. The obtained thermal noise is shown to ensure equilibration of all degrees of freedom in a simulation to high accuracy. Furthermore, we demonstrate that satisfactory results for most practical applications of fluctuating hydrodynamics can already be achieved using thermal noise derived in the long wavelength-limit.

preprint2009arXiv

Biaxiality at the Isotropic-Nematic Interface with Planar Anchoring

We revisit the classic problem of the structure of the isotropic-nematic interface within Ginzburg-Landau-de Gennes theory, refining previous analytic treatments of biaxiality at the interface. We compare our analysis with numerical results obtained through a highly accurate spectral collocation scheme for the solution of the Landau-Ginzburg-de Gennes equations. In comparison to earlier work, we obtain improved agreement with numerics for both the uniaxial and biaxial profiles, accurate asymptotic results for the decay of biaxial order on both nematic and isotropic sides of the interface and accurate fits to data from density functional approaches to this problem.

preprint2009arXiv

Statistical analysis of the Indus script using $n$-grams

The Indus script is one of the major undeciphered scripts of the ancient world. The small size of the corpus, the absence of bilingual texts, and the lack of definite knowledge of the underlying language has frustrated efforts at decipherment since the discovery of the remains of the Indus civilisation. Recently, some researchers have questioned the premise that the Indus script encodes spoken language. Building on previous statistical approaches, we apply the tools of statistical language processing, specifically $n$-gram Markov chains, to analyse the Indus script for syntax. Our main results are that the script has well-defined signs which begin and end texts, that there is directionality and strong correlations in the sign order, and that there are groups of signs which appear to have identical syntactic function. All these require no {\it a priori} suppositions regarding the syntactic or semantic content of the signs, but follow directly from the statistical analysis. Using information theoretic measures, we find the information in the script to be intermediate between that of a completely random and a completely fixed ordering of signs. Our study reveals that the Indus script is a structured sign system showing features of a formal language, but, at present, cannot conclusively establish that it encodes {\it natural} language. Our $n$-gram Markov model is useful for predicting signs which are missing or illegible in a corpus of Indus texts. This work forms the basis for the development of a stochastic grammar which can be used to explore the syntax of the Indus script in greater detail.

preprint2009arXiv

The Isotropic-Nematic Interface with an Oblique Anchoring Condition

We present numerical and analytic results for uniaxial and biaxial order at the isotropic-nematic interface within Ginzburg-Landau-de Gennes theory. We study the case where an oblique anchoring condition is imposed asymptotically on the nematic side of the interface, reproducing results of previous work when this condition reduces to planar or homoeotropic anchoring. We construct physically motivated and computationally flexible variational profiles for uniaxial and biaxial order, comparing our variational results to numerical results obtained from a minimization of the Ginzburg-Landau-de Gennes free energy. While spatial variations of the scalar uniaxial and biaxial order parameters are confined to the neighbourhood of the interface, nematic elasticity requires that the director orientation interpolate linearly between either planar or homoeotropic anchoring at the location of the interface and the imposed boundary condition at infinity. The selection of planar or homoeotropic anchoring at the interface is governed by the sign of the Ginzburg-Landau-de Gennes elastic coefficient L_2. Our variational calculations are in close agreement with our numerics and agree qualitatively with results from density functional theory and molecular simulations.