Source author record

Sourav Chatterjee

Sourav Chatterjee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

88works

30topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2023arXiv

A computationally efficient and mechanically compatible multi-phase-field model applied to coherently stressed three-phase solids

Engineering alloys generally exhibit multi-phase microstructures. For simulating their microstructure evolution during solid-state phase transformation, CALPHAD-guided multi-phase-field models coupled with micro-mechanics have proven to be a reliable simulation tool. Nevertheless, their efficiency and accuracy still depend on the homogenization scheme used to interpolate the elastic properties in the interfacial regions. In this paper, we present a phase-field model for multi-phase and multi-component solids using a partial rank-one homogenization scheme that enforces static and kinematic compatibilities in the interfacial regions. To this end, we first extend the rank-one homogenization scheme to multi-phase systems. Moreover, for computational efficiency, we analytically solve the static compatibility equations for linear elastic three-phase solids. For quantitative accuracy, a coupling technique is used to extract the prerequisite thermodynamic and kinetic properties from CALPHAD databases. The model is solved numerically in an open source finite-element framework. As numerical applications, the microstructure of two elastically stressed intermetallic-containing three-phase alloys: Ni-Al and Al-Cr-Ni, are simulated. The accuracy of the model is verified against analytically obtained solutions for planar and concentric ring interfaces. We show that the simulation results remain unaltered with varying interface width. Except for one simulation, all cases show better or nearly equal convergence using the partial rank-one scheme compared to the Voigt-Taylor scheme.

preprint2023arXiv

An efficient and quantitative phase-field model for elastically heterogeneous two-phase solids based on a partial rank-one homogenization scheme

This paper presents an efficient and quantitative phase-field model for elastically heterogeneous alloys that ensures the two mechanical compatibilities$\unicode{x2014}$static and kinematic, in conjunction with chemical equilibrium within the interfacial region. Our model contrasts with existing phase-field models that either violate static compatibility or interfacial chemical equilibrium or are computationally costly. For computational efficiency, the partial rank-one homogenization (PRH) scheme is employed to enforce both static and kinematic compatibilities at the interface. Moreover, interfacial chemical equilibrium is ensured by replacing the composition field with the diffusion potential field as the independent variable of the model. Its performance is demonstrated by simulating four single-particle and one multi-particle cases for two binary two-phase alloys: Ni-Al $γ^{\prime}/γ$ and UO$_2$/void. Its accuracy is then investigated against analytical solutions. For the single-particle $γ^{\prime}/γ$ alloy, we find that the accuracy of the phase-field results remains unaffected for both planar and non-planar geometries when the PRH scheme is employed. Fortuitously, in the UO$_2$/void simulations, despite a strong elastic heterogeneity$\unicode{x2014}$the ratio of Young's modulus of the void phase to that of the UO$_2$ phase is $10^{-4}\unicode{x2014}$we find that the PRH scheme shows significantly better convergence compared to the Voigt-Taylor scheme (VTS) for both planar and non-planar geometries. Nevertheless, for the same interface width range as in the $γ^{'}/γ$ case, the interface migration in these simulations shows dependence on interface width. Contrary to the $γ^{\prime}/γ$ simulations, we also find that the simulated elastic fields show deviations from the analytical solution in the non-planar UO$_2$/void case using the PRH scheme.

preprint2022arXiv

A random walk on the Rado graph

The Rado graph, also known as the random graph $G(\infty, p)$, is a classical limit object for finite graphs. We study natural ball walks as a way of understanding the geometry of this graph. For the walk started at $i$, we show that order $\log_2^*i$ steps are sufficient, and for infinitely many $i$, necessary for convergence to stationarity. The proof involves an application of Hardy's inequality for trees.

preprint2022arXiv

Compact Object Modeling in the Globular Cluster 47 Tucanae

The globular cluster 47~Tucanae (47~Tuc) is one of the most massive star clusters in the Milky Way and is exceptionally rich in exotic stellar populations. For several decades it has been a favorite target of observers, and yet it is computationally very challenging to model because of its large number of stars ($N\gtrsim 10^6$) and high density. Here we present detailed and self-consistent 47~Tuc models computed with the \texttt{Cluster Monte Carlo} code (\texttt{CMC}). The models include all relevant dynamical interactions coupled to stellar and binary evolution, and reproduce various observations, including the surface brightness and velocity dispersion profiles, pulsar accelerations, and numbers of compact objects. We show that the present properties of 47~Tuc are best reproduced by adopting an initial stellar mass function that is both bottom-heavy and top-light relative to standard assumptions \citep[as in, e.g.,][]{Kroupa2001}, and an initial Elson profile \citep{Elson1987} that is overfilling the cluster's tidal radius. We include new prescriptions in \texttt{CMC} for the formation of binaries through giant star collisions and tidal captures, and we show that these mechanisms play a crucial role in the formation of neutron star binaries and millisecond pulsars in 47~Tuc; our best-fit model contains $\sim 50$ millisecond pulsars, $70\%$ of which are formed through giant collisions and tidal captures. Our models also suggest that 47~Tuc presently contains up to $\sim 200$ stellar-mass black holes, $\sim 5$ binary black holes, $\sim 15$ low-mass X-ray binaries, and $\sim 300$ cataclysmic variables.

preprint2022arXiv

Existence of stationary ballistic deposition on the infinite lattice

Ballistic deposition is one of the many models of interface growth that are believed to be in the KPZ universality class, but have so far proved to be largely intractable mathematically. In this model, blocks of size one fall independently as Poisson processes at each site on the $d$-dimensional lattice, and either attach themselves to the column growing at that site, or to the side of an adjacent column, whichever comes first. It is not hard to see that if we subtract off the height of the column at the origin from the heights of the other columns, the resulting interface process is Markovian. The main result of this article is that this Markov process has at least one invariant probability measure. We conjecture that the invariant measure is not unique, and provide some partial evidence.

preprint2022arXiv

Globular Clusters UVIT Legacy Survey (GlobULeS) I. FUV-optical Color-Magnitude Diagrams for Eight Globular Clusters

We present the first results of eight Globular Clusters (GCs) from the AstroSat/UVIT Legacy Survey program GlobULeS based on the observations carried out in two FUV filters (F148W and F169M). The FUV-optical and FUV-FUV color-magnitude diagrams (CMDs) of GCs with the proper motion membership were constructed by combining the UVIT data with HST UV Globular Cluster Survey (HUGS) data for inner regions and Gaia Early Data Release (EDR3) for regions outside the HST's field. We detect sources as faint as F148W $\sim$ 23.5~mag which are classified based on their locations in CMDs by overlaying stellar evolutionary models. The CMDs of 8 GCs are combined with the previous UVIT studies of 3 GCs to create stacked FUV-optical CMDs to highlight the features/peculiarities found in the different evolutionary sequences. The FUV (F148W) detected stellar populations of 11 GCs comprises 2,816 Horizontal Branch (HB) stars (190 Extreme HB candidates), 46 post-HB (pHB), 221 Blue Straggler Stars (BSS), and 107 White Dwarf (WD) candidates. We note that the blue HB color extension obtained from F148W$-$G color and the number of FUV detected EHB candidates are strongly correlated with the maximum internal Helium (He) variation within each GC, suggesting that the FUV-optical plane is the most sensitive to He abundance variations in the HB. We discuss the potential science cases that will be addressed using these catalogues including HB morphologies, BSSs, pHB, and, WD stars.

preprint2022arXiv

Isomorphisms between random graphs

Consider two independent Erdős-Rényi $G(N,1/2)$ graphs. We show that with probability tending to $1$ as $N\to\infty$, the largest induced isomorphic subgraph has size either $\lfloor x_N-\varepsilon_N\rfloor$ or $\lfloor x_N+\varepsilon_N \rfloor$, where $x_N=4\log_2 N -2 \log_2 \log_2 N - 2\log_2(4/e)+1$ and $\varepsilon_N = (4\log_2 N)^{-1/2}$. Using similar techniques, we also show that if $Γ_1$ and $Γ_2$ are independent $G(n,1/2)$ and $G(N,1/2)$ random graphs, then $Γ_2$ contains an isomorphic copy of $Γ_1$ as an induced subgraph with high probability if $n\le \lfloor y_N - \varepsilon_N \rfloor$ and does not contain an isomorphic copy of $Γ_1$ as an induced subgraph with high probability if $n>\lfloor y_N+\varepsilon_N \rfloor$, where $y_N=2\log_2 N+1$ and $\varepsilon_N$ is as above.

preprint2022arXiv

Matrix completion with data-dependent missingness probabilities

The problem of completing a large matrix with lots of missing entries has received widespread attention in the last couple of decades. Two popular approaches to the matrix completion problem are based on singular value thresholding and nuclear norm minimization. Most of the past works on this subject assume that there is a single number $p$ such that each entry of the matrix is available independently with probability $p$ and missing otherwise. This assumption may not be realistic for many applications. In this work, we replace it with the assumption that the probability that an entry is available is an unknown function $f$ of the entry itself. For example, if the entry is the rating given to a movie by a viewer, then it seems plausible that high value entries have greater probability of being available than low value entries. We propose two new estimators, based on singular value thresholding and nuclear norm minimization, to recover the matrix under this assumption. The estimators involve no tuning parameters, and are shown to be consistent under a low rank assumption. We also provide a consistent estimator of the unknown function $f$.

preprint2022arXiv

MOSPAT: AutoML based Model Selection and Parameter Tuning for Time Series Anomaly Detection

Organizations leverage anomaly and changepoint detection algorithms to detect changes in user behavior or service availability and performance. Many off-the-shelf detection algorithms, though effective, cannot readily be used in large organizations where thousands of users monitor millions of use cases and metrics with varied time series characteristics and anomaly patterns. The selection of algorithm and parameters needs to be precise for each use case: manual tuning does not scale, and automated tuning requires ground truth, which is rarely available. In this paper, we explore MOSPAT, an end-to-end automated machine learning based approach for model and parameter selection, combined with a generative model to produce labeled data. Our scalable end-to-end system allows individual users in large organizations to tailor time-series monitoring to their specific use case and data characteristics, without expert knowledge of anomaly detection algorithms or laborious manual labeling. Our extensive experiments on real and synthetic data demonstrate that this method consistently outperforms using any single algorithm.

preprint2022arXiv

Polarization correction towards satellite-based QKD without an active feedback

Quantum key distribution (QKD) is a cryptographic protocol to enable two parties to share a secure key string, which can be used in one-time pad cryptosystem. There has been an ongoing surge of interest in implementing long-haul photonic-implementation of QKD protocols. However, the endeavour is challenging in many aspects. In particular, one of the major challenges is the polarization degree of freedom of single-photons getting affected while transmission through optical fibres, or atmospheric turbulence. Conventionally, an active feedback-based mechanism is employed to achieve real-time polarization tracking. In this work, we propose an alternative approach where we first perform a state tomography to reconstruct the output density matrix. We then evaluate the optimal measurement bases at Bob's end that leads to the maximum (anti-)correlation in the measurement outcomes of both parties. As a proof-of-principle demonstration, we implement an in-lab BBM92 protocol -- a particular variant of a QKD protocol using quantum entanglement as a resource -- to exemplify the performance of our technique. We experimentally generate polarization-entangled photon pairs having $94\%$ fidelity with $\ketψ_1 = 1/\sqrt{2}\,(\ket{HV}+\ket{VH})$ state and a concurrence of $0.92$. By considering a representative 1 ns coincidence window span, we are able to achieve a quantum-bit-error-rate (QBER) of $\approx 5\%$, and a key rate of $\approx 35$ Kbps. The protocol performance is independent of local polarization rotations through optical fibres. We also develop an algorithmic approach to optimize the trade-off between the key rate and QBER. Our approach obviates the need for active polarization tracking. Our method is also applicable to entanglement-based QKD demonstrations using partially mixed as well as non-maximally entangled states, and extends to single-photon implementations over fibre channels.

preprint2022arXiv

Superconcentration in surface growth

Height functions of growing random surfaces are often conjectured to be superconcentrated, meaning that their variances grow sublinearly in time. This article introduces a new concept, called subroughness, meaning that there exist two distinct points such that the expected squared difference between the heights at these points grows sublinearly in time. The main result of the paper is that superconcentration is equivalent to subroughness in a class of growing random surfaces. The result is applied to establish superconcentration in a variant of the restricted solid-on-solid (RSOS) model and in a variant of the ballistic deposition model, and give new proofs of superconcentration in directed last-passage percolation and directed polymers.

preprint2022arXiv

The Yang-Mills heat flow with random distributional initial data

We construct local solutions to the Yang-Mills heat flow (in the DeTurck gauge) for a certain class of random distributional initial data, which includes the 3D Gaussian free field. The main idea, which goes back to work of Bourgain as well as work of Da Prato-Debussche, is to decompose the solution into a rougher linear part and a smoother nonlinear part, and to control the latter by probabilistic arguments. In a companion work, we use the main results of this paper to propose a way towards the construction of 3D Yang-Mills measures.

preprint2022arXiv

Weak convergence of directed polymers to deterministic KPZ at high temperature

It is shown that when $d\ge 3$, the growing random surface generated by the $(d+1)$-dimensional directed polymer model at sufficiently high temperature, after being smoothed by taking microscopic local averages, converges to a solution of the deterministic KPZ equation in a suitable scaling limit.

preprint2021arXiv

Convergence of deterministic growth models

We prove the uniform in space and time convergence of the scaled heights of large classes of deterministic growth models that are monotone and equivariant under translations by constants. The limits are characterized as the unique (viscosity solutions) of first- or second-order partial differential equations depending on whether the growth models are scaled hyperbolically or parabolically. The results greatly simplify and extend a recent work by the first author to more general surface growth models. The proofs are based on the methodology developed by Barles and the second author to prove convergence of approximation schemes.

preprint2021arXiv

Intermediate-mass Black Holes from High Massive-star Binary Fractions in Young Star Clusters

Black holes formed in dense star clusters, where dynamical interactions are frequent, may have fundamentally different properties than those formed through isolated stellar evolution. Theoretical models for single star evolution predict a gap in the black hole mass spectrum from roughly $40-120\,M_{\odot}$ caused by (pulsational) pair-instability supernovae. Motivated by the recent LIGO/Virgo event GW190521, we investigate whether black holes with masses within or in excess of this "upper-mass gap" can be formed dynamically in young star clusters through strong interactions of massive stars in binaries. We perform a set of $N$-body simulations using the CMC cluster-dynamics code to study the effects of the high-mass binary fraction on the formation and collision histories of the most massive stars and their remnants. We find that typical young star clusters with low metallicities and high binary fractions in massive stars can form several black holes in the upper-mass gap and often form at least one intermediate-mass black hole. These results provide strong evidence that dynamical interactions in young star clusters naturally lead to the formation of more massive black hole remnants.

preprint2021arXiv

Micro-Estimates of Wealth for all Low- and Middle-Income Countries

Many critical policy decisions, from strategic investments to the allocation of humanitarian aid, rely on data about the geographic distribution of wealth and poverty. Yet many poverty maps are out of date or exist only at very coarse levels of granularity. Here we develop the first micro-estimates of wealth and poverty that cover the populated surface of all 135 low and middle-income countries (LMICs) at 2.4km resolution. The estimates are built by applying machine learning algorithms to vast and heterogeneous data from satellites, mobile phone networks, topographic maps, as well as aggregated and de-identified connectivity data from Facebook. We train and calibrate the estimates using nationally-representative household survey data from 56 LMICs, then validate their accuracy using four independent sources of household survey data from 18 countries. We also provide confidence intervals for each micro-estimate to facilitate responsible downstream use. These estimates are provided free for public use in the hope that they enable targeted policy response to the COVID-19 pandemic, provide the foundation for new insights into the causes and consequences of economic development and growth, and promote responsible policymaking in support of the Sustainable Development Goals.

preprint2021arXiv

Modeling Dense Star Clusters in the Milky Way and Beyond with the Cluster Monte Carlo Code

We describe the public release of the Cluster Monte Carlo Code (CMC) a parallel, star-by-star $N$-body code for modeling dense star clusters. CMC treats collisional stellar dynamics using Hénon's method, where the cumulative effect of many two-body encounters is statistically reproduced as a single effective encounter between nearest-neighbor particles on a relaxation timescale. The star-by-star approach allows for the inclusion of additional physics, including strong gravitational three- and four-body encounters, two-body tidal and gravitational-wave captures, mass loss in arbitrary galactic tidal fields, and stellar evolution for both single and binary stars. The public release of CMC is pinned directly to the COSMIC population synthesis code, allowing dynamical star cluster simulations and population synthesis studies to be performed using identical assumptions about the stellar physics and initial conditions. As a demonstration, we present two examples of star cluster modeling: first, we perform the largest ($N = 10^8$) star-by-star $N$-body simulation of a Plummer sphere evolving to core collapse, reproducing the expected self-similar density profile over more than 15 orders of magnitude; second, we generate realistic models for typical globular clusters, and we show that their dynamical evolution can produce significant numbers of black hole mergers with masses greater than those produced from isolated binary evolution (such as GW190521, a recently reported merger with component masses in the pulsational pair-instability mass gap).

preprint2021arXiv

The Observed Rate of Binary Black Hole Mergers can be Entirely Explained by Globular Clusters

Since the first signal in 2015, the gravitational-wave detections of merging binary black holes (BBHs) by the LIGO and Virgo collaborations (LVC) have completely transformed our understanding of the lives and deaths of compact object binaries, and have motivated an enormous amount of theoretical work on the astrophysical origin of these objects. We show that the phenomenological fit to the redshift-dependent merger rate of BBHs from Abbott et al. (2020) is consistent with a purely dynamical origin for these objects, and that the current merger rate of BBHs from the LVC could be explained entirely with globular clusters alone. While this does not prove that globular clusters are the dominant formation channel, we emphasize that many formation scenarios could contribute a significant fraction of the current LVC rate, and that any analysis that assumes a single (or dominant) mechanism for producing BBH mergers is implicitly using a specious astrophysical prior.

preprint2020arXiv

A Dynamical Survey of Stellar-Mass Black Holes in 50 Milky Way Globular Clusters

Recent numerical simulations of globular clusters (GCs) have shown that stellar-mass black holes (BHs) play a fundamental role in driving cluster evolution and shaping their present-day structure. Rapidly mass-segregating to the center of GCs, BHs act as a dynamical energy source via repeated super-elastic scattering, delaying onset of core collapse and limiting mass segregation for visible stars. While recent discoveries of BH candidates in Galactic and extragalactic GCs have further piqued interest in BH-mediated cluster dynamics, numerical models show that even if significant BH populations remain in today's GCs, they are typically in configurations that are not directly detectable. We demonstrated in Weatherford et al. (2018) that an anti-correlation between a suitable measure of mass segregation ($Δ$) in observable stellar populations and the number of retained BHs in GC models can be applied to indirectly probe BH populations in real GCs. Here, we estimate the number and total mass of BHs presently retained in 50 Milky Way GCs from the ACS Globular Cluster Survey by measuring $Δ$ between populations of main sequence stars, using correlations found between $Δ$ and BH retention in the CMC Cluster Catalog models. We demonstrate that the range in $Δ$'s distribution from our models matches that for observed GCs to a remarkable degree. Our results further provide the narrowest constraints to-date on the retained BH populations in the GCs analyzed. Of these 50 GCs, we identify NGCs 2808, 5927, 5986, 6101, and 6205 to presently contain especially large BH populations, each with total BH mass exceeding $10^3\,\rm{M_{\odot}}$.

preprint2020arXiv

A new coefficient of correlation

Is it possible to define a coefficient of correlation which is (a) as simple as the classical coefficients like Pearson's correlation or Spearman's correlation, and yet (b) consistently estimates some simple and interpretable measure of the degree of dependence between the variables, which is 0 if and only if the variables are independent and 1 if and only if one is a measurable function of the other, and (c) has a simple asymptotic theory under the hypothesis of independence, like the classical coefficients? This article answers this question in the affirmative, by producing such a coefficient. No assumptions are needed on the distributions of the variables. There are several coefficients in the literature that converge to 0 if and only if the variables are independent, but none that satisfy any of the other properties mentioned above.

preprint2020arXiv

Demographics of triple systems in dense star clusters

Depending on the stellar type, more than $\sim 50$\% and $\sim 15\%$ of stars in the field have at least one and two stellar companions, respectively. Hierarchical systems can be assembled dynamically in dense star clusters, as a result of few-body encounters among stars and/or compact remnants in the cluster core. In this paper, we present the demographics of stellar and compact-object triples formed via binary--binary encounters in the \texttt{CMC Cluster Catalog}, a suite of cluster simulations with present-day properties representative of the globular clusters (GCs) observed in the Milky Way. We show how the initial properties of the host cluster set the typical orbital parameters and formation times of the formed triples. We find that a cluster typically assembles hundreds of triples with at least one black hole (BH) in the inner binary, while only clusters with sufficiently small virial radii are efficient in producing triples with no BHs, as a result of the BH-burning process. We show that a typical GC is expected to host tens of triples with at least one luminous component at present day. We discuss how the Lidov-Kozai mechanism can drive the inner binary of the formed triples to high eccentricities, whenever it takes place before the triple is dynamically reprocessed by encountering another cluster member. Some of these systems can reach sufficiently large eccentricities to form a variety of transients and sources, such as blue stragglers, X-ray binaries, Type Ia Supernovae, Thorne-Zytkow objects, and LIGO/Virgo sources.

preprint2020arXiv

GW190412 as a Third-Generation Black Hole Merger from a Super Star Cluster

We explore the possibility that GW190412, a binary black hole merger with a non-equal-mass ratio and significantly spinning primary, was formed through repeated black hole mergers in a dense super star cluster. Using a combination of semi-analytic prescriptions for the remnant spin and recoil kick of black hole mergers, we show that the mass ratio and spin of GW190412 are consistent with a binary black hole whose primary component has undergone two successive mergers from a population of $\sim 10M_{\odot}$ black holes in a high-metallicity environment. We then explore the production of GW190412-like analogs in the CMC Cluster Catalog, a grid of 148 $N$-body star cluster models, as well as a new model, behemoth, with nearly $10^7$ particles and initial conditions taken from a cosmological MHD simulation of galaxy formation. We show that the production of binaries with GW190412-like masses and spins is dominated by massive super star clusters with high metallicities and large central escape speeds. While many are observed in the local universe, our results suggest that a careful treatment of these massive clusters, many of which may have been disrupted before the present day, is necessary to characterize the production of unique gravitational-wave events produced through dynamics.

preprint2020arXiv

Impact of Local Stiffness on Entropy Driven Microscopic Dynamics of Polythiophene

We exploited the high temporal and spatial resolution of neutron spin echo spectroscopy to investigate the large-scale dynamics of semiflexible conjugated polymer chains in solutions. We obtained the first experimental demonstration of earlier predicted single chain glassy state. We used a generalized approach of the well-established Zimm model of flexible polymers to describe the relaxation mode spectra of locally stiff polythiophene chains. The Zimm mode analysis confirms the existence of beads with a finite length that corresponds to a reduced number of segmental modes in semiflexible chains. Irrespective of the temperature and the molecular weight of the conjugated polymer, we witness a universal behavior of the local chain stiffness and invariability of the beads length. Our experimental findings indicate possibly minor role of the change in π-electron conjugation length (and therefore conjugated backbone planar to non-planar conformational transition) in the observed thermochromic behavior of polythiophene but instead point on the major role of chain dynamics in this phenomenon.

preprint2020arXiv

Modeling Dense Star Clusters in the Milky Way and Beyond with the $\texttt{CMC}$ Cluster Catalog

We present a set of 148 independent $N$-body simulations of globular clusters (GCs) computed using the code $\texttt{CMC}$ ($\texttt{Cluster Monte Carlo}$). At an age of $\sim10-13\,$Gyr, the resulting models cover nearly the full range of cluster properties exhibited by the Milky Way GCs, including total mass, core and half-light radii, metallicity, and galactocentric distance. We use our models to investigate the role that stellar-mass black holes play in the process of core collapse. Furthermore, we study how dynamical interactions affect the formation and evolution of several important types of sources in GCs, including low-mass X-ray binaries, millisecond pulsars, blue stragglers, cataclysmic variables, Type Ia supernovae, calcium-rich transients, and merging compact binaries. While our focus here is on old, low-metallicity GCs, our $\texttt{CMC}$ simulations follow the evolution of clusters over a Hubble time, and they include a wide range of metallicities (up to solar), so that our results can also be used to study younger and higher-metallicity star clusters.

preprint2020arXiv

Speeding up Markov chains with deterministic jumps

We show that the convergence of finite state space Markov chains to stationarity can often be considerably speeded up by alternating every step of the chain with a deterministic move. Under fairly general conditions, we show that not only do such schemes exist, they are numerous.

preprint2020arXiv

Wilson loops in Ising lattice gauge theory

Wilson loop expectation in 4D $\mathbb{Z}_2$ lattice gauge theory is computed to leading order in the weak coupling regime. This is the first example of a rigorous theoretical calculation of Wilson loop expectation in the weak coupling regime of a 4D lattice gauge theory. All prior results are either inequalities or strong coupling expansions.

preprint2019arXiv

On the Rate of Neutron Star Binary Mergers from Globular Clusters

The first detection of gravitational waves from a neutron star - neutron star (NS-NS) merger, GW170817, and the increasing number of observations of short gamma-ray bursts (SGRBs) have greatly motivated studies of the origins of NS-NS and neutron star - black hole (NS-BH) binaries. We calculate the merger rates of NS-NS and NS-BH binaries from globular clusters (GCs) using realistic GC simulations with the \texttt{CMC} cluster catalog. We use a large sample of models with a range of initial numbers of stars, metallicities, virial radii and galactocentric distances, representative of the present-day Milky Way GCs, to quantify the inspiral times and volumetric merger rates as a function of redshift, both inside and ejected from clusters. We find that over the complete lifetime of most GCs, stellar BHs dominate the cluster cores and prevent the mass segregation of NSs, thereby reducing the dynamical interaction rates of NSs so that at most a few NS binary mergers are ever produced. We estimate the merger rate in the local universe to be $\sim\rm{0.02\,Gpc^{-3}\,yr^{-1}}$ for both NS-NS and NS-BH binaries, or a total of $\sim 0.04$~Gpc$^{-3}$~yr$^{-1}$ for both populations. These rates are about 5 orders of magnitude below the current empirical merger rate from LIGO/Virgo. We conclude that dynamical interactions in GCs do not play a significant role in enhancing the NS-NS and NS-BH merger rates.

preprint2019arXiv

qkdSim: An experimenter's simulation toolkit for QKD with imperfections, and its performance analysis with a demonstration of the B92 protocol using heralded photon

Quantum Key Distribution (QKD) is one of the most important aspects of quantum cryptography. Using laws of quantum mechanics as the basis for security, the key distribution process is made information theoretically secure in QKD. With the advancement and commercialization of QKD, an end-to-end QKD simulation software is required that can include experimental imperfections. Software of this kind will ensure that resources are invested only after prior performance analysis, and is faithful to experimental capacities and limitations. In this work, we introduce our QKD simulation toolkit qkdSim, which is ultimately aimed at being developed into such a software package that can precisely model and analyse any generic QKD protocol. We present the design, implementation and testing of a prototype of qkdSim that can accurately simulate our own experimental demonstration of the B92 protocol. The simulation results match well with experiment; a representative key rate and QBER from experiment is $51 \pm 0.5$ Kbits/sec and $4.79\% \pm 0.01\%$ respectively, wherein the simulation yields $52.83 \pm 0.36$ Kbits/sec and $4.79\% \pm 0.01\%$ respectively.

preprint2019arXiv

The Role of "Black Hole Burning" in the Evolution of Dense Star Clusters

As self-gravitating systems, dense star clusters exhibit a natural diffusion of energy from their innermost to outermost regions, which leads to a slow and steady contraction of the core until it ultimately collapses under gravity. However, in spite of the natural tendency toward so-called "core collapse," the globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the core-collapsed and non-core-collapsed clusters. This suggests an internal energy source is at work, delaying the onset of core collapse in many clusters. Primordial binary stars have been thought for a long time to provide this energy source, but recent analyses have cast doubt upon the corresponding "binary-burning" mechanism as a viable explanation. Over the past decade, a large amount of both observational and theoretical work has suggested that many stellar-mass black holes (BHs) are retained in typical clusters today and that they play a dynamically-significant role in these clusters throughout their entire lifetimes. Here we review our latest understanding of the formation and evolution of BH populations in GCs and demonstrate that, through their dynamical interaction with their host cluster, BHs can naturally explain the distinction between core-collapsed and non-core-collapsed clusters through a process we call "black hole burning."

preprint2016arXiv

A central limit theorem for a new statistic on permutations

This paper does three things: It proves a central limit theorem for novel permutation statistics (for example, the number of descents plus the number of descents in the inverse). It provides a clear illustration of a new approach to proving central limit theorems more generally. It gives us an opportunity to acknowledge the work of our teacher and friend B. V. Rao.

preprint2016arXiv

A note about the uniform distribution on the intersection of a simplex and a sphere

Uniform probability distributions on $\ell_p$ balls and spheres have been studied extensively and are known to behave like product measures in high dimensions. In this note we consider the uniform distribution on the intersection of a simplex and a sphere. Certain new and interesting features, such as phase transitions and localization phenomena emerge.

preprint2016arXiv

An introduction to large deviations for random graphs

This article gives an overview of the emerging literature on large deviations for random graphs. Written for the general mathematical audience, the article begins with a short introduction to the theory of large deviations. This is followed by a description of some large deviation questions about random graphs, and an outline of the recent progress on this topic. A more elaborate discussion follows, with a brief account of graph limit theory and its application in constructing a large deviation theory for dense random graphs. The role of Szemerédi's regularity lemma is explained, together with a sketch of the proof of the main large deviation result and some examples. Applications to exponential random graph models are briefly touched upon. The remainder of the paper is devoted to large deviations for sparse graphs. Since the regularity lemma is not applicable in the sparse regime, new tools are needed. Fortunately, there have been several new breakthroughs that managed to achieve the goal by an indirect method. These are discussed, together with an exposition of the underlying theory. The last section contains a list of open problems.

preprint2016arXiv

Binary Black Hole Mergers from Globular Clusters: Masses, Merger Rates, and the Impact of Stellar Evolution

The recent discovery of GW150914, the binary black hole merger detected by Advanced LIGO, has the potential to revolutionize observational astrophysics. But to fully utilize this new window into the universe, we must compare these new observations to detailed models of binary black hole formation throughout cosmic time. Expanding upon our previous work (Rodriguez et al., 2015), we study merging binary black holes formed in globular clusters using our Monte Carlo approach to stellar dynamics. We have created a new set of 52 cluster models with different masses, metallicities, and radii to fully characterize the binary black hole merger rate. These models include all the relevant dynamical processes (such as two-body relaxation, strong encounters, and three-body binary formation) and agree well with detailed direct N-body simulations. In addition, we have enhanced our stellar evolution algorithms with updated metallicity-dependent stellar wind and supernova prescriptions, allowing us to compare our results directly to the most recent population synthesis predictions for merger rates from isolated binary evolution. We explore the relationship between a cluster's global properties and the population of binary black holes that it produces. In particular, we derive a numerically calibrated relationship between the merger times of ejected black hole binaries and a cluster's mass and radius. We explore the masses and mass ratios of these binaries as a function of redshift, and find a merger rate of ~5 Gpc$^{-3}$ yr$^{-1}$ in the local universe, with 80% of sources having total masses from $32M_{\odot}$ to $64M_{\odot}$. Under standard assumptions, approximately 1 out of every 7 binary black hole mergers in the local universe will have originated in a globular cluster, but we also explore the sensitivity of this result to different assumptions for binary stellar evolution. [Abridged]

preprint2016arXiv

Broadcasting of Quantum Correlations: Possibilities & Impossibilities

In this work, we extensively study the problem of broadcasting of quantum correlations. This includes broadcasting of quantum entanglement as well as correlations that go beyond the notion of entanglement. It is quite well known from the "No-Broadcasting theorem" that perfect broadcasting of quantum correlation is not possible. However it does not rule out the possibility of partial broadcasting of correlations where we can get lesser correlated states from a given correlated state. In order to have a holistic view of broadcasting, we investigate this problem by starting with most general representation of two qubit mixed states in terms of the Bloch vectors. As a cloning transformation we have used universal symmetric optimal Buzek-Hillery (B-H) cloner both locally and nonlocally. More specifically, we obtain a set of ranges in terms of Bloch vectors for which broadcasting of entanglement will be possible. In addition to the idea of broadcasting of entanglement for general two qubit mixed states, we explore broadcasting of quantum correlations that go beyond entanglement with the help of both local and nonlocal cloners. Remarkably, we find that it is impossible to even partially broadcast such correlations by using both local and nonlocal quantum copying machines. Taking two different types Buzek-Hillery quantum cloners (state dependent and state independent) we analytically prove the impossibility of even partial broadcasting of correlation. Lastly, we generalize this impossibility result for any symmetric or asymmetric cloning machines as well. This result brings out a fundamental difference between the correlation defined from the perspective of entanglement and the correlation measure which claims to go beyond entanglement.

preprint2016arXiv

Complementarity in Quantum Information Processing Tasks

Complementarity have been an intriguing feature of physical systems for a long time. In this work we establish a new kind of complimentary relations in the frame work of quantum information processing tasks. In broadcasting of entanglement we create many pairs of less entangled states from a given entangled state both by local and non local cloning operations. These entangled states can be used in various information processing tasks like teleportation and superdense coding. Since these states are less entangled states it is quite intuitive that these states are not going to be as powerful resource as the initial states. In this work we study the usefulness of these states in tasks like teleportation and super dense coding. More precisely, we found out bounds of their capabilities in terms of several complimentary relations involving fidelity of broadcasting. In principle we have considered general mixed as a resource also separately providing different examples like a) Werner like states, b) Bell diagonal states. Here we have used both local and non local cloning operations as means of broadcasting. In the later part of our work, we extend this result by obtaining bounds in form of complimentary relations in a situation where we have used $1-N$ cloning transformations instead of $1-2$ cloning transformations.

preprint2016arXiv

Concentration inequalities with exchangeable pairs (Ph.D. thesis)

The purpose of this dissertation is to introduce a version of Stein's method of exchangeable pairs to solve problems in measure concentration. We specifically target systems of dependent random variables, since that is where the power of Stein's method is fully realized. Because the theory is quite abstract, we have tried to put in as many examples as possible. Some of the highlighted applications are as follows: (a) We shall find an easily verifiable condition under which a popular heuristic technique originating from physics, known as the "mean field equations" method, is valid. No such condition is currently known. (b) We shall present a way of using couplings to derive concentration inequalities. Although couplings are routinely used for proving decay of correlations, no method for using couplings to derive concentration bounds is available in the literature. This will be used to obtain (c) concentration inequalities with explicit constants under Dobrushin's condition of weak dependence. (d) We shall give a method for obtaining concentration of Haar measures using convergence rates of related random walks on groups. Using this technique and one of the numerous available results about rates of convergence of random walks, we will then prove (e) a quantitative version of Voiculescu's celebrated connection between random matrix theory and free probability.

preprint2016arXiv

Dynamical Formation of the GW150914 Binary Black Hole

We explore the possibility that GW150914, the binary black hole merger recently detected by Advanced LIGO, was formed by gravitational interactions in the core of a dense star cluster. Using models of globular clusters with detailed $N$-body dynamics and stellar evolution, we show that a typical cluster with a mass of $3\times10^5M_{\odot}$ to $6\times10^5M_{\odot}$ is optimal for forming GW150914-like binary black holes that will merge in the local universe. We identify the most likely dynamical processes for forming GW150914 in such a cluster, and we show that the detection of GW150914 is consistent with the masses and merger rates expected for binary black hole mergers from globular clusters. Our results show that dynamical processes provide a significant and well-understood pathway for forming binary black hole mergers in the local universe. Understanding the contribution of dynamics to the binary black hole merger problem is a critical step in unlocking the full potential of gravitational-wave astronomy.

preprint2016arXiv

Million-Body Star Cluster Simulations: Comparisons between Monte Carlo and Direct $N$-body

We present the first detailed comparison between million-body globular cluster simulations computed with a Hénon-type Monte Carlo code, CMC, and a direct $N$-body code, NBODY6++GPU. Both simulations start from an identical cluster model with $10^6$ particles, and include all of the relevant physics needed to treat the system in a highly realistic way. With the two codes "frozen" (no fine-tuning of any free parameters or internal algorithms of the codes) we find excellent agreement in the overall evolution of the two models. Furthermore, we find that in both models, large numbers of stellar-mass black holes (> 1000) are retained for 12 Gyr. Thus, the very accurate direct $N$-body approach confirms recent predictions that black holes can be retained in present-day, old globular clusters. We find only minor disagreements between the two models and attribute these to the small-$N$ dynamics driving the evolution of the cluster core for which the Monte Carlo assumptions are less ideal. Based on the overwhelming general agreement between the two models computed using these vastly different techniques, we conclude that our Monte Carlo approach, which is more approximate, but dramatically faster compared to the direct $N$-body, is capable of producing a very accurate description of the long-term evolution of massive globular clusters even when the clusters contain large populations of stellar-mass black holes.

preprint2016arXiv

Minimal spanning trees and Stein's method

Kesten and Lee [36] proved that the total length of a minimal spanning tree on certain random point configurations in $\mathbb{R}^d$ satisfies a central limit theorem. They also raised the question: how to make these results quantitative? However, techniques employed to tackle the same problem for other functionals studied in geometric probability do not apply directly to the minimal spanning tree. Thus the problem of determining the convergence rate in the central limit theorem for Euclidean minimal spanning trees has remained open. In this work, we establish bounds on the convergence rate for the Poissonized version of this problem by using a variation of Stein's method. We also derive bounds on the convergence rate for the analogous problem in the setup of the lattice $\mathbb{Z}^d$. The contribution of this paper is twofold. First, we develop a general technique to compute convergence rates in central limit theorems satisfied by minimal spanning trees on sequence of weighted graphs which includes minimal spanning trees on Poisson points. Secondly, we present a way of quantifying the Burton-Keane argument for the uniqueness of the infinite open cluster. The latter is interesting in its own right and based on a generalization of our technique, Duminil-Copin, Ioffe and Velenik [28] have recently obtained bounds on probability of two-arm events in a broad class of translation-invariant percolation models.

preprint2016arXiv

Nonlinear large deviations

We present a general technique for computing large deviations of nonlinear functions of independent Bernoulli random variables. The method is applied to compute the large deviation rate functions for subgraph counts in sparse random graphs. Previous technology, based on Szemeredi's regularity lemma, works only for dense graphs. Applications are also made to exponential random graphs and three-term arithmetic progressions in random sets of integers.

preprint2016arXiv

Prediction error of cross-validated Lasso

In spite of the wealth of literature on the theoretical properties of the Lasso, there is very little known when the value of the tuning parameter is chosen using the data, even though this is what actually happens in practice. We give a general upper bound on the prediction error of Lasso when the tuning parameter is chosen using a variant of 2-fold cross-validation. No special assumption is made about the structure of the design matrix, and the tuning parameter is allowed to be optimized over an arbitrary data-dependent set of values. The proof is based on a general principle that may extend to other kinds of cross-validation as well as to other penalized regression methods. Based on this result, we propose a new estimate for error variance in high dimensional regression and prove that it has good properties under minimal assumptions.

preprint2016arXiv

Sequential Quantum Secret Sharing in a Noisy Environment aided with Weak Measurements

In this work we give a $(n,n)$-threshold protocol for sequential secret sharing of quantum information for the first time. By sequential secret sharing we refer to a situation where the dealer is not having all the secrets at the same time, at the beginning of the protocol; however if the dealer wishes to share secrets at subsequent phases she/he can realize it with the help of our protocol. First of all we present our protocol for three parties and later we generalize it for the situation where we have $(n>3)$ parties. Further in a much more realistic situation, we consider the sharing of qubits through two kinds of noisy channels, namely the phase damping channel (PDC) and the amplitude damping channel (ADC). When we carry out the sequential secret sharing in the presence of noise we observe that the fidelity of secret sharing at the $k^{th}$ iteration is independent of the effect of noise at the $(k-1)^{th}$ iteration. In case of ADC we have seen that the average fidelity of secret sharing drops down to $\frac{1}{2}$ which is equivalent to a random guess of the quantum secret. Interestingly, we find that by applying weak measurements one can enhance the average fidelity. This increase of the average fidelity can be achieved with certain trade off with the success probability of the weak measurements.

preprint2016arXiv

The $1/N$ expansion for SO(N) lattice gauge theory at strong coupling

The $1/N$ expansion is an asymptotic series expansion for certain quantities in large-$N$ lattice gauge theories. This article gives a rigorous formulation and proof of the $1/N$ expansion for Wilson loop expectations in SO(N) lattice gauge theory in the strong coupling regime in any dimension. The terms in the expansion are expressed as sums over trajectories of strings in a lattice string theory, establishing an explicit gauge-string duality. The trajectories trace out surfaces of genus zero for the first term in the expansion, and surfaces of higher genus for the higher terms.

preprint2016arXiv

The leading term of the Yang-Mills free energy

This article gives an explicit formula for the leading term of the free energy of three-dimensional $U(N)$ lattice gauge theory for any $N$, as the lattice spacing tends to zero. The proof is based on a novel technique that avoids phase cell renormalization. The technique also yields a similar formula for the four-dimensional theory, but only in the weak coupling~limit.

preprint2015arXiv

An Overview of Inside-Out Planet Formation

The Kepler-discovered Systems with Tightly-packed Inner Planets (STIPs), typically with several planets of Earth to super-Earth masses on well-aligned, sub-AU orbits may host the most common type of planets, including habitable planets, in the Galaxy. They pose a great challenge for planet formation theories, which fall into two broad classes: (1) formation further out followed by inward migration; (2) formation in situ, in the very inner regions of the protoplanetary disk. We review the pros and cons of these classes, before focusing on a new theory of sequential in situ formation from the inside-out via creation of successive gravitationally unstable rings fed from a continuous stream of small (~cm-m size) "pebbles," drifting inward via gas drag. Pebbles first collect at the pressure trap associated with the transition from a magnetorotational instability (MRI)-inactive ("dead zone") region to an inner, MRI-active zone. A pebble ring builds up that begins to dominate the local mass surface density of the disk and spawns a planet. The planet continues to grow, most likely by pebble accretion, until it becomes massive enough to isolate itself from the accretion flow via gap opening. This reduces the local gas density near the planet, leading to enhanced ionization and a retreat of the dead zone inner boundary. The process repeats with a new pebble ring gathering at the new pressure maximum associated with this boundary. We discuss the theory's predictions for planetary masses, relative mass scalings with orbital radius, and minimum orbital separations, and their comparison with observed systems. Finally, we discuss open questions, including potential causes of diversity of planetary system architectures, i.e., STIPs versus Solar System analogs.

preprint2015arXiv

Binary Black Hole Mergers from Globular Clusters: Implications for Advanced LIGO

The predicted rate of binary black hole mergers from galactic fields can vary over several orders of magnitude and is extremely sensitive to the assumptions of stellar evolution. But in dense stellar environments such as globular clusters, binary black holes form by well-understood gravitational interactions. In this letter, we study the formation of black hole binaries in an extensive collection of realistic globular cluster models. By comparing these models to observed Milky Way and extragalactic globular clusters, we find that the mergers of dynamically-formed binaries could be detected at a rate of ~100 per year, potentially dominating the binary black hole merger rate. We also find that a majority of cluster-formed binaries are more massive than their field-formed counterparts, suggesting that Advanced LIGO could identify certain binaries as originating from dense stellar environments.

preprint2015arXiv

Black hole mergers and blue stragglers from hierarchical triples formed in globular clusters

Hierarchical triple-star systems are expected to form frequently via close binary-binary encounters in the dense cores of globular clusters. In a sufficiently inclined triple, gravitational interactions between the inner and outer binary can cause large-amplitude oscillations in the eccentricity of the inner orbit ("Lidov-Kozai cycles"), which can lead to a collision and merger of the two inner components. In this paper we use Monte Carlo models of dense star clusters to identify all triple systems formed dynamically and we compute their evolution using a highly accurate three-body integrator which incorporates relativistic and tidal effects. We find that a large fraction of these triples evolve through a non-secular dynamical phase which can drive the inner binary to higher eccentricities than predicted by the standard secular perturbation theory (even including octupole-order terms). We place constraints on the importance of Lidov-Kozai-induced mergers for producing: (i) gravitational wave sources detectable by Advanced LIGO (aLIGO), for triples with an inner pair of stellar black holes; and (ii) blue straggler stars, for triples with main-sequence-star components. We find a realistic aLIGO detection rate of black hole mergers due to the Lidov-Kozai mechanism of 1yr^-1, with about 20% of these having a finite eccentricity when they first chirp into the aLIGO frequency band. While rare, these events are likely to dominate among eccentric compact object inspirals that are potentially detectable by aLIGO. For blue stragglers, we find that the Lidov-Kozai mechanism can contribute only up to ~10% of their total numbers in globular clusters.

preprint2015arXiv

High dimensional regression and matrix estimation without tuning parameters

A general theory for Gaussian mean estimation that automatically adapts to unknown sparsity under arbitrary norms is proposed. The theory is applied to produce adaptively minimax rate-optimal estimators in high dimensional regression and matrix estimation that involve no tuning parameters.

preprint2015arXiv

Inside-Out Planet Formation. III. Planet-disk interaction at the dead zone inner boundary

The Kepler mission has discovered more than 4000 exoplanet candidates. Many are in systems with tightly packed inner planets. Inside-Out Planet Formation (IOPF) has been proposed to explain these systems. It involves sequential in situ planet formation at the local pressure maximum of a retreating dead zone inner boundary (DZIB). Pebbles accumulate at this pressure trap, which builds up a ring, and then a planet. The planet is expected to grow until it opens a gap, which helps to both truncate pebble accretion and induce DZIB retreat that sets the location of formation of the next planet. This simple scenario may be modified if the planet migrates significantly from its formation location. Thus planet-disk interactions play a crucial role in the IOPF scenario. We present numerical simulations that first assess migration of planets of various masses that are forming at the DZIB of an active accretion disk, where the effective viscosity rapidly increases in the radially inward direction. We find that the disk's torques on the planet tend to trap the planet at a location very close to the initial pressure maximum where it formed. We then study gap opening by these planets to assess at what mass a significant gap is created. Finally we present a simple model for DZIB retreat due to penetration of X-rays from the star to the disk midplane. Overall, these simulations help to quantify both the mass scale of first,"Vulcan," planet formation and the orbital separation to the location of second planet formation.

preprint2015arXiv

Period Ratio Distribution of Near-Resonant Planets Indicates Planetesimal Scattering

An intriguing trend among Kepler's multi-planet systems is an overabundance of planet pairs with period ratios just wide of mean motion resonances (MMR) and a dearth of systems just narrow of them. In a recently published paper Chatterjee & Ford (2015; henceforth CF15) has proposed that gas-disk migration traps planets in a MMR. After gas dispersal, orbits of these trapped planets are altered through interaction with a residual planetesimal disk. They found that for massive enough disks planet-planetesimal disk interactions can break resonances and naturally create moderate to large positive offsets from the initial period ratio for large ranges of planetesimal disk and planet properties. Divergence from resonance only happens if the mass of planetesimals that interact with the planets is at least a few percent of the total planet mass. This threshold, above which resonances are broken and the offset from resonances can grow, naturally explains why the asymmetric large offsets were not seen in more massive planet pairs found via past radial velocity surveys. In this article we will highlight some of the key findings of CF15. In addition, we report preliminary results from an extension of this study, that investigates the effects of planet-planetesimal disk interactions on initially non-resonant planet pairs. We find that planetesimal scattering typically increases period ratios of non-resonant planets. If the initial period ratios are below and in proximity of a resonance, under certain conditions, this increment in period ratios can create a deficit of systems with period ratios just below the exact integer corresponding to the MMR and an excess just above. From an initially uniform distribution of period ratios just below a 2:1 MMR, planetesimal interactions can create an asymmetric distribution across this MMR similar to what is observed for the Kepler planet pairs.

preprint2015arXiv

Planetesimal Interactions Can Explain the Mysterious Period Ratios of Small Near-Resonant Planets

An intriguing trend among \kepler's multi-planet systems is an overabundance of planet pairs with period ratios just wide of a mean motion resonance (MMR) and a dearth of systems just narrow of them. Traditional planet formation models are at odds with these observations. They are also in contrast with the period ratios of radial-velocity-discovered multi-planet systems which tend to pile up at 2:1 MMR. We propose that gas-disk migration traps planets in a MMR. After gas dispersal, orbits of these trapped planets are altered through interaction with a residual planetesimal disk. We study the effects of planetesimal disk interactions on planet pairs trapped in 2:1 MMR using planets of mass typical of the Kepler planet candidates (KPC) and explore large ranges for the mass, and density profile of the planetesimal disk. We find that planet-planetesimal disk interactions naturally create the observed asymmetry in period-ratio distribution for large ranges of planetesimal disk and planet properties. If the planetesimal disk mass is above a threshold of ~0.2x the planet mass, these interactions typically disrupt MMR. Afterwards, the planets migrate in such a way that the final period-ratio is slightly higher than the integer ratio corresponding to the initial MMR. Below this threshold these interactions typically cannot disrupt the resonance and the period ratio stays close to the integer ratio. The threshold explains why the more massive planet pairs found by RV surveys are still in resonance. We encourage future research to explore how significantly the associated accretion would change the planets' atmospheric and surface properties.

preprint2014arXiv

A new perspective on least squares under convex constraint

Consider the problem of estimating the mean of a Gaussian random vector when the mean vector is assumed to be in a given convex set. The most natural solution is to take the Euclidean projection of the data vector on to this convex set; in other words, performing "least squares under a convex constraint." Many problems in modern statistics and statistical signal processing theory are special cases of this general situation. Examples include the lasso and other high-dimensional regression techniques, function estimation problems, matrix estimation and completion, shape-restricted regression, constrained denoising, linear inverse problems, etc. This paper presents three general results about this problem, namely, (a) an exact computation of the main term in the estimation error by relating it to expected maxima of Gaussian processes (existing results only give upper bounds), (b) a theorem showing that the least squares estimator is always admissible up to a universal constant in any problem of the above kind and (c) a counterexample showing that least squares estimator may not always be minimax rate-optimal. The result from part (a) is then used to compute the error of the least squares estimator in two examples of contemporary interest.

preprint2014arXiv

A short survey of Stein's method

Stein's method is a powerful technique for proving central limit theorems in probability theory when more straightforward approaches cannot be implemented easily. This article begins with a survey of the historical development of Stein's method and some recent advances. This is followed by a description of a "general purpose" variant of Stein's method that may be called the generalized perturbative approach, and an application of this method to minimal spanning trees. The article concludes with the descriptions of some well known open problems that may possibly be solved by the perturbative approach or some other variant of Stein's method.

preprint2014arXiv

Absence of replica symmetry breaking in the random field Ising model

It is shown that replica symmetry is not broken in the random field Ising model in any dimension, at any temperature and field strength, except possibly at a measure-zero set of exceptional temperatures and field strengths.

preprint2014arXiv

Assumptionless consistency of the Lasso

The Lasso is a popular statistical tool invented by Robert Tibshirani for linear regression when the number of covariates is greater than or comparable to the number of observations. The purpose of this note is to highlight the simple fact (noted in a number of earlier papers in various guises) that for the loss function considered in Tibshirani's original paper, the Lasso is consistent under almost no assumptions at all.

preprint2014arXiv

Fluctuations of the Bose-Einstein condensate

This article gives a rigorous analysis of the fluctuations of the Bose-Einstein condensate for a system of non-interacting bosons in an arbitrary potential, assuming that the system is governed by the canonical ensemble. As a result of the analysis, we are able to tell the order of fluctuations of the condensate fraction as well as its limiting distribution upon proper centering and scaling. This yields interesting results. For example, for a system of $n$ bosons in a 3D harmonic trap near the transition temperature, the order of fluctuations of the condensate fraction is $n^{-1/2}$ and the limiting distribution is normal, whereas for the 3D uniform Bose gas, the order of fluctuations is $n^{-1/3}$ and the limiting distribution is an explicit non-normal distribution. For a 2D harmonic trap, the order of fluctuations is $n^{-1/2}(\log n)^{1/2}$, which is larger than $n^{-1/2}$ but the limiting distribution is still normal. All of these results come as easy consequences of a general theorem.

preprint2014arXiv

Matrix estimation by Universal Singular Value Thresholding

Consider the problem of estimating the entries of a large matrix, when the observed entries are noisy versions of a small random fraction of the original entries. This problem has received widespread attention in recent times, especially after the pioneering works of Emmanuel Candès and collaborators. This paper introduces a simple estimation procedure, called Universal Singular Value Thresholding (USVT), that works for any matrix that has "a little bit of structure." Surprisingly, this simple estimator achieves the minimax error rate up to a constant factor. The method is applied to solve problems related to low rank matrix estimation, blockmodels, distance matrix completion, latent space models, positive definite matrix completion, graphon estimation and generalized Bradley--Terry models for pairwise comparison.

preprint2014arXiv

Pebble Delivery for Inside-Out Planet Formation

Inside-Out Planet Formation (IOPF; Chatterjee & Tan 2014, hereafter CT14) is a scenario for sequential in situ planet formation at the pressure traps of retreating dead zone inner boundaries (DZIBs) motivated to explain the many systems with tightly packed inner planets (STIPs) discovered by Kepler. The scenario involves build-up of a pebble-dominated protoplanetary ring, supplied by radial drift of pebbles from the outer disk. It may also involve further build-up of planetary masses to gap-opening scales via continued pebble accretion. Here we study radial drift & growth of pebbles delivered to the DZIB in fiducial IOPF disk models.

preprint2014arXiv

Vulcan Planets: Inside-Out Formation of the Innermost Super-Earths

The compact multi-transiting systems discovered by Kepler challenge traditional planet formation theories. These fall into two broad classes: (1) formation further out followed by migration; (2) formation in situ from a disk of gas and planetesimals. In the former, an abundance of resonant chains is expected, which the Kepler data do not support. In the latter, required disk mass surface densities may be too high. A recently proposed mechanism hypothesizes that planets form in situ at the pressure trap associated with the dead-zone inner boundary (DZIB) where radially drifting "pebbles" accumulate. This scenario predicts planet masses ($M_p$) are set by the gap-opening process that then leads to DZIB retreat, followed by sequential, inside-out planet formation (IOPF). For typical disk accretion rates, IOPF predictions for $M_p$, $M_p$ versus orbital radius $r$, and planet-planet separations are consistent with observed systems. Here we investigate the IOPF prediction for how the masses, $M_{p,1}$, of the innermost ("Vulcan") planets vary with $r$. We show that for fiducial parameters, $M_{p,1}\simeq5.0(r/{\rm{0.1\:AU}})\:M_\oplus$, independent of the disk's accretion rate at time of planet formation. Then, using Monte Carlo sampling of a population of these innermost planets, we test this predicted scaling against observed planet properties, allowing for intrinsic dispersions in planetary densities and Kepler's observational biases. These effects lead to a slightly shallower relation $M_{p,1}\propto{r}^{0.9\pm0.2}$, which is consistent with $M_{p,1}\propto{r}^{0.7\pm0.2}$ of the observed Vulcans. The normalization of the relation constrains the gap-opening process, favoring relatively low viscosities in the inner dead zone.

preprint2013arXiv

Estimating and understanding exponential random graph models

We introduce a method for the theoretical analysis of exponential random graph models. The method is based on a large-deviations approximation to the normalizing constant shown to be consistent using theory developed by Chatterjee and Varadhan [European J. Combin. 32 (2011) 1000-1017]. The theory explains a host of difficulties encountered by applied workers: many distinct models have essentially the same MLE, rendering the problems ``practically'' ill-posed. We give the first rigorous proofs of ``degeneracy'' observed in these models. Here, almost all graphs have essentially no edges or are essentially complete. We supplement recent work of Bhamidi, Bresler and Sly [2008 IEEE 49th Annual IEEE Symposium on Foundations of Computer Science (FOCS) (2008) 803-812 IEEE] showing that for many models, the extra sufficient statistics are useless: most realizations look like the results of a simple Erdős-Rényi model. We also find classes of models where the limiting graphs differ from Erdős-Rényi graphs. A limitation of our approach, inherited from the limitation of graph limit theory, is that it works only for dense graphs.

preprint2013arXiv

Inside-Out Planet Formation

The compact multi-transiting planet systems discovered by Kepler challenge planet formation theories. Formation in situ from disks with radial mass surface density, $Σ$, profiles similar to the minimum mass solar nebula (MMSN) but boosted in normalization by factors $\gtrsim 10$ has been suggested. We propose that a more natural way to create these planets in the inner disk is formation sequentially from the inside-out via creation of successive gravitationally unstable rings fed from a continuous stream of small (~cm--m size) "pebbles", drifting inwards via gas drag. Pebbles collect at the pressure maximum associated with the transition from a magneto-rotational instability (MRI)-inactive ("dead zone") region to an inner MRI-active zone. A pebble ring builds up until it either becomes gravitationally unstable to form an $\sim 1\ M_\oplus$ planet directly or induces gradual planet formation via core accretion. The planet may undergo Type I migration into the active region, allowing a new pebble ring and planet to form behind it. Alternatively if migration is inefficient, the planet may continue to accrete from the disk until it becomes massive enough to isolate itself from the accretion flow. A variety of densities may result depending on the relative importance of residual gas accretion as the planet approaches its isolation mass. The process can repeat with a new pebble ring gathering at the new pressure maximum associated with the retreating dead zone boundary. Our simple analytical model for this scenario of inside-out planet formation yields planetary masses, relative mass scalings with orbital radius, and minimum orbital separations consistent with those seen by Kepler. It provides an explanation of how massive planets can form with tightly-packed and well-aligned system architectures, starting from typical protoplanetary disk properties.

preprint2013arXiv

Large Deviations for Random Matrices

We prove a large deviation result for a random symmetric n x n matrix with independent identically distributed entries to have a few eigenvalues of size n. If the spectrum S survives when the matrix is rescaled by a factor of n, it can only be the eigenvalues of a Hilbert-Schmidt kernel k(x,y) on [0,1] x [0,1]. The rate function for k is $I(k)=1/2\int h(k(x,y) dxdy$ where h is the Cramer rate function for the common distribution of the entries that is assumed to have a tail decaying faster than any Gaussian. The large deviation for S is then obtained by contraction.

preprint2013arXiv

Mock Observations of Blue Stragglers in Globular Cluster Models

We created artificial color-magnitude diagrams of Monte Carlo dynamical models of globular clusters, and then used observational methods to determine the number of blue stragglers in those clusters. We compared these blue stragglers to various cluster properties, mimicking work that has been done for blue stragglers in Milky Way globular clusters to determine the dominant formation mechanism(s) of this unusual stellar population. We find that a mass-based prescription for selecting blue stragglers will choose approximately twice as many blue stragglers than a selection criterion that was developed for observations of real clusters. However, the two numbers of blue stragglers are well-correlated, so either selection criterion can be used to characterize the blue straggler population of a cluster. We confirm previous results that the simplified prescription for the evolution of a collision or merger product in the BSE code overestimates their lifetimes. We show that our model blue stragglers follow similar trends with cluster properties (core mass, binary fraction, total mass, collision rate) as the true Milky Way blue stragglers, as long as we restrict ourselves to model clusters with an initial binary fraction higher than 5%. We also show that, in contrast to earlier work, the number of blue stragglers in the cluster core does have a weak dependence on the collisional parameter Gamma in both our models and in Milky Way globular clusters.

preprint2013arXiv

On level sets of Gaussian fields

In this short note, we present a theorem concerning certain "additive structure" for the level sets of non-degenerate Gaussian fields, which yields the multiple valley phenomenon for extremal fields with exponentially many valleys.

preprint2013arXiv

Stellar Collisions and Blue Straggler Stars in Dense Globular Clusters

Blue straggler stars (BSS) are abundantly observed in all Galactic globular clusters (GGC) where data exist. However, observations alone cannot reveal the relative importance of various formation channels or the typical formation times for this well studied population of anomalous stars. Using a state-of-the-art Hénon-type Monte Carlo code that includes all relevant physical processes, we create 128 models with properties typical of the observed GGCs. These models include realistic numbers of single and binary stars, use observationally motivated initial conditions, and span large ranges in central density, concentration, binary fraction, and mass. Their properties can be directly compared with those of observed GGCs. We can easily identify the BSSs in our models and determine their formation channels and birth times. We find that for central densities above ~10^3 Msun/pc^3 the dominant formation channel is stellar collisions while for lower density clusters, mass transfer in binaries provides a significant contribution (up to ~ 60% in our models). The majority of these collisions are binary-mediated, occurring during 3-body and 4-body interactions. As a result a strong correlation between the specific frequency of BSSs and the binary fraction in a cluster can be seen in our models. We find that the number of BSSs in the core shows only a weak correlation with the collision rate estimator Γtraditionally used by observers, in agreement with the latest Hubble Space Telescope (ACS) data. Using an idealized "full mixing" prescription for collision products, our models indicate that the BSSs observed today may have formed several Gyrs ago. However, denser clusters tend to have younger (~1 Gyr) BSSs.

preprint2013arXiv

Stochastic solutions of the wave equation

Unlike the heat equation or the Laplace equation, solutions of the wave equation on general domains have no known stochastic representation. This short note gives a simple solution to this well known problem in arbitrary dimensions. The proposed representation has several shortcomings, one of which is that it does not cover all solutions. Still, it is proof that a large class of nontrivial solutions of the wave equation in general dimensions and domains may indeed be represented stochastically.

preprint2012arXiv

Central limit theorem for first-passage percolation time across thin cylinders

We prove that first-passage percolation times across thin cylinders of the form $[0,n]\times [-h_n,h_n]^{d-1}$ obey Gaussian central limit theorems as long as $h_n$ grows slower than $n^{1/(d+1)}$. It is an open question as to what is the fastest that $h_n$ can grow so that a Gaussian CLT still holds. Under the natural but unproven assumption about existence of fluctuation and transversal exponents, and strict convexity of the limiting shape in the direction of $(1,0,...,0)$, we prove that in dimensions 2 and 3 the CLT holds all the way up to the height of the unrestricted geodesic. We also provide some numerical evidence in support of the conjecture in dimension 2.

preprint2012arXiv

Gravitational Slingshot of Young Massive Stars in Orion

The Orion Nebula Cluster (ONC) is the nearest region of massive star formation and thus a crucial testing ground for theoretical models. Of particular interest amongst the ONC's ~1000 members are: θ^1 Ori C, the most massive binary in the cluster with stars of masses 38 and 9 MSun (Kraus et al. 2009); the Becklin-Neugebauer (BN) object, a 30 km/s runaway star of ~8 MSun (Tan 2004); and the Kleinmann-Low (KL) nebula protostar, a highly-obscured, ~15 MSun object still accreting gas while also driving a powerful, apparently "explosive" outflow (Allen & Burton 1993). The unusual behavior of BN and KL is much debated: How did BN acquire its high velocity? How is this related to massive star formation in the KL nebula? Here we report the results of a systematic survey using ~ 10^7 numerical experiments of gravitational interactions of the θ^1C and BN stars. We show that dynamical ejection of BN from this triple system at its observed velocity leaves behind a binary with total energy and eccentricity matching those observed for θ^1C. Five other observed properties of θ^C are also consistent with it having ejected BN and altogether we estimate there is only a <~ 10^{-5} probability that θ^1C has these properties by chance. We conclude that BN was dynamically ejected from the θ^1C system about 4,500 years ago. BN has then plowed through the KL massive-star-forming core within the last 1,000 years causing its recently-enhanced accretion and outflow activity.

preprint2012arXiv

Monte Carlo Simulations of Globular Cluster Evolution. VI. The Influence of an Intermediate Mass Black Hole

We present results of a series of Monte Carlo simulations investigating the imprint of a central intermediate-mass black hole (IMBH) on the structure of a globular cluster. We investigate the three-dimensional and projected density profiles, and stellar disruption rates for idealized as well as realistic cluster models, taking into account a stellar mass spectrum and stellar evolution, and allowing for a larger, more realistic, number of stars than was previously possible with direct N-body methods. We compare our results to other N-body and Fokker-Planck simulations published previously. We find, in general, very good agreement for the overall cluster structure and dynamical evolution between direct N-body simulations and our Monte Carlo simulations. Significant differences exist in the number of stars that are tidally disrupted by the IMBH, which is most likely an effect of the wandering motion of the IMBH, not included in the Monte Carlo scheme. These differences, however, are negligible for the final IMBH masses in realistic cluster models as the disruption rates are generally much lower than for single-mass clusters. As a direct comparison to observations we construct a detailed model for the cluster NGC 5694, which is known to possess a central surface brightness cusp consistent with the presence of an IMBH. We find that not only the inner slope but also the outer part of the surface brightness profile agree well with observations. However, there is only a slight preference for models harboring an IMBH compared to models without.

preprint2012arXiv

Planets in Open Clusters Detectable by Kepler

While hundreds of planets have been discovered around field stars, only a few are known in star clusters. To explain the lack of short-period giant planets in globular clusters (GC), such as 47 Tucane and ω Centauri, it has been suggested that their low metallicities may have prevented planet formation. Alternatively, the high rates of close stellar encounters in these clusters may have influenced the formation and subsequent evolution of planetary systems. How common are planets in clusters around normal main-sequence stars? Here we consider whether this question can be addressed using data from the Kepler mission. The Kepler field of view contains 4 low-density (relative to GCs) open clusters where the metallicities are about solar (or even higher) and stellar encounters are much less frequent than in typical GCs. We provide detailed $N$-body models and show that most planets in Kepler-detectable orbits are not significantly perturbed by stellar encounters in these open clusters. We focus on the most massive cluster, NGC 6791, which has super-solar metallicity, and find that if planets formed in this cluster at the same frequency as observed in the field, Kepler could detect 1 -- 20 transiting planets depending on the planet-size distribution and the duration of data collection. Due to the large distance to NGC 6791 Kepler will have to search relatively faint ($K_p<20$) stars for the full extended mission to achieve such a yield.

preprint2012arXiv

Random Overlap Structures: Properties and Applications to Spin Glasses

Random Overlap Structures (ROSt's) are random elements on the space of probability measures on the unit ball of a Hilbert space, where two measures are identified if they differ by an isometry. In spin glasses, they arise as natural limits of Gibbs measures under the appropriate algebra of functions. We prove that the so called `cavity mapping' on the space of ROSt's is continuous, leading to a proof of the stochastic stability conjecture for the limiting Gibbs measures of a large class of spin glass models. Similar arguments yield the proofs of a number of other properties of ROSt's that may be useful in future attempts at proving the ultrametricity conjecture. Lastly, assuming that the ultrametricity conjecture holds, the setup yields a constructive proof of the Parisi formula for the free energy of the Sherrington-Kirkpatrick model by making rigorous a heuristic of Aizenman, Sims and Starr.

preprint2012arXiv

The universal relation between scaling exponents in first-passage percolation

It has been conjectured in numerous physics papers that in ordinary first-passage percolation on integer lattices, the fluctuation exponent $χ$ and the wandering exponent $ξ$ are related through the universal relation $χ=2ξ-1$, irrespective of the dimension. This is sometimes called the KPZ relation between the two exponents. This article gives a rigorous proof of this conjecture assuming that the exponents exist in a certain sense.

preprint2012arXiv

Understanding the Dynamical State of Globular Clusters: Core-Collapsed vs Non Core-Collapsed

We study the dynamical evolution of globular clusters using our Hénon-type Monte Carlo code for stellar dynamics including all relevant physics such as two-body relaxation, single and binary stellar evolution, Galactic tidal stripping, and strong interactions such as physical collisions and binary mediated scattering. We compute a large database of several hundred models starting from broad ranges of initial conditions guided by observations of young and massive star clusters. We show that these initial conditions very naturally lead to present day clusters with properties including the central density, core radius, half-light radius, half-mass relaxation time, and cluster mass, that match well with those of the old Galactic globular clusters. In particular, we can naturally reproduce the bimodal distribution in observed core radii separating the "core-collapsed" vs the "non core-collapsed" clusters. We see that the core-collapsed clusters are those that have reached or are about to reach the equilibrium "binary burning" phase. The non core-collapsed clusters are still undergoing gravo-thermal contraction.

preprint2011arXiv

Nonnormal approximation by Stein's method of exchangeable pairs with application to the Curie--Weiss model

Let $(W,W')$ be an exchangeable pair. Assume that \[E(W-W'|W)=g(W)+r(W),\] where $g(W)$ is a dominated term and $r(W)$ is negligible. Let $G(t)=\int_0^tg(s)\,ds$ and define $p(t)=c_1e^{-c_0G(t)}$, where $c_0$ is a properly chosen constant and $c_1=1/\int_{-\infty}^{\infty}e^{-c_0G(t)}\,dt$. Let $Y$ be a random variable with the probability density function $p$. It is proved that $W$ converges to $Y$ in distribution when the conditional second moment of $(W-W')$ given $W$ satisfies a law of large numbers. A Berry-Esseen type bound is also given. We use this technique to obtain a Berry-Esseen error bound of order $1/\sqrt{n}$ in the noncentral limit theorem for the magnetization in the Curie-Weiss ferromagnet at the critical temperature. Exponential approximation with application to the spectrum of the Bernoulli-Laplace Markov chain is also discussed.

preprint2011arXiv

Probabilistic methods for discrete nonlinear Schrödinger equations

We show that the thermodynamics of the focusing cubic discrete nonlinear Schrodinger equation are exactly solvable in dimensions three and higher. A number of explicit formulas are derived.

preprint2011arXiv

Random graphs with a given degree sequence

Large graphs are sometimes studied through their degree sequences (power law or regular graphs). We study graphs that are uniformly chosen with a given degree sequence. Under mild conditions, it is shown that sequences of such graphs have graph limits in the sense of Lovász and Szegedy with identifiable limits. This allows simple determination of other features such as the number of triangles. The argument proceeds by studying a natural exponential model having the degree sequence as a sufficient statistic. The maximum likelihood estimate (MLE) of the parameters is shown to be unique and consistent with high probability. Thus $n$ parameters can be consistently estimated based on a sample of size one. A fast, provably convergent, algorithm for the MLE is derived. These ingredients combine to prove the graph limit theorem. Along the way, a continuous version of the Erdős--Gallai characterization of degree sequences is derived.

preprint2011arXiv

Random multiplicative functions in short intervals

We consider random multiplicative functions taking the values $\pm 1$. Using Stein's method for normal approximation, we prove a central limit theorem for the sum of such multiplicative functions in appropriate short intervals.

preprint2011arXiv

Spectral clustering and the high-dimensional stochastic blockmodel

Networks or graphs can easily represent a diverse set of data sources that are characterized by interacting units or actors. Social networks, representing people who communicate with each other, are one example. Communities or clusters of highly connected actors form an essential feature in the structure of several empirical networks. Spectral clustering is a popular and computationally feasible method to discover these communities. The stochastic blockmodel [Social Networks 5 (1983) 109--137] is a social network model with well-defined communities; each node is a member of one community. For a network generated from the Stochastic Blockmodel, we bound the number of nodes "misclustered" by spectral clustering. The asymptotic results in this paper are the first clustering results that allow the number of clusters in the model to grow with the number of nodes, hence the name high-dimensional. In order to study spectral clustering under the stochastic blockmodel, we first show that under the more general latent space model, the eigenvectors of the normalized graph Laplacian asymptotically converge to the eigenvectors of a "population" normalized graph Laplacian. Aside from the implication for spectral clustering, this provides insight into a graph visualization technique. Our method of studying the eigenvectors of random matrices is original.

preprint2011arXiv

The large deviation principle for the Erdős-Rényi random graph

What does an Erdos-Renyi graph look like when a rare event happens? This paper answers this question when p is fixed and n tends to infinity by establishing a large deviation principle under an appropriate topology. The formulation and proof of the main result uses the recent development of the theory of graph limits by Lovasz and coauthors and Szemeredi's regularity lemma from graph theory. As a basic application of the general principle, we work out large deviations for the number of triangles in G(n,p). Surprisingly, even this simple example yields an interesting double phase transition.

preprint2011arXiv

The missing log in large deviations for triangle counts

This paper solves the problem of sharp large deviation estimates for the upper tail of the number of triangles in an Erdos-Renyi random graph, by establishing a logarithmic factor in the exponent that was missing till now. It is possible that the method of proof may extend to general subgraph counts.

preprint2010arXiv

A new approach to strong embeddings

We revisit strong approximation theory from a new perspective, culminating in a proof of the Komlós-Major-Tusnády embedding theorem for the simple random walk. The proof is almost entirely based on a series of soft arguments and easy inequalities. The new technique, inspired by Stein's method of normal approximation, is applicable to any setting where Stein's method works. In particular, one can hope to take it beyond sums of independent random variables.

preprint2010arXiv

Applications of Stein's method for concentration inequalities

Stein's method for concentration inequalities was introduced to prove concentration of measure in problems involving complex dependencies such as random permutations and Gibbs measures. In this paper, we provide some extensions of the theory and three applications: (1) We obtain a concentration inequality for the magnetization in the Curie--Weiss model at critical temperature (where it obeys a nonstandard normalization and super-Gaussian concentration). (2) We derive exact large deviation asymptotics for the number of triangles in the Erdős--Rényi random graph $G(n,p)$ when $p\ge0.31$. Similar results are derived also for general subgraph counts. (3) We obtain some interesting concentration inequalities for the Ising model on lattices that hold at all temperatures.

preprint2010arXiv

How planet-planet scattering can create high-inclination as well as long-period orbits

Recent observations have revealed two new classes of planetary orbits. Rossiter- Mclaughlin (RM) measurements have revealed hot Jupiters in high-obliquity orbits. In addition, direct-imaging has discovered giant planets at large (~ 100 AU) separations via direct-imaging technique. Simple-minded disk-migration scenarios are inconsistent with the high-inclination (and even retrograde) orbits as seen in recent RM measurements. Furthermore, forming giant planets at large semi-major axis (a) may be challenging in the core-accretion paradigm. We perform many N-body simulations to explore the two above-mentioned orbital architectures. Planet-planet scattering in a multi-planet system can naturally excite orbital inclinations. Planets can also get scattered to large distances. Large-a planetary orbits created from planet-planet scattering are expected to have high eccentricities (e). Theoretical models predict that the observed long-period planets, such as Fomalhaut-b have moderate e \approx 0.3. Interestingly, these are also in systems with disks. We find that if a massive-enough outer disk is present, a scattered planet may be circularized at large a via dynamical friction from the disk and repeated scattering of the disk particles.

preprint2010arXiv

Properties of Uniform Doubly Stochastic Matrices

We investigate the properties of uniform doubly stochastic random matrices, that is non-negative matrices conditioned to have their rows and columns sum to 1. The rescaled marginal distributions are shown to converge to exponential distributions and indeed even large sub-matrices of side-length $o(n^{1/2-ε})$ behave like independent exponentials. We determine the limiting empirical distribution of the singular values the the matrix. Finally the mixing time of the associated Markov chains is shown to be exactly 2 with high probability.

preprint2010arXiv

Unstable Planetary Systems Emerging Out Of Gas Disks

The discovery of over 400 extrasolar planets allows us to statistically test our understanding of formation and dynamics of planetary systems via numerical simulations. Traditional N-body simulations of multiple-planet systems without gas disks have successfully reproduced the eccentricity (e) distribution of the observed systems, by assuming that the planetary systems are relatively closely packed when the gas disk dissipates, so that they become dynamically unstable within the stellar lifetime. However, such studies cannot explain the small semi-major axes (a) of extrasolar planetary systems, if planets are formed, as the standard planet formation theory suggests, beyond the ice line. In this paper, we numerically study the evolution of three-planet systems in dissipating gas disks, and constrain the initial conditions that reproduce the observed semi-major axis and eccentricity distributions simultaneously. We adopt the initial conditions that are motivated by the standard planet formation theory, and self-consistently simulate the disk evolution, and planet migration by using a hybrid N-body and 1D gas disk code. We also take account of eccentricity damping, and investigate the effect of saturation of corotation resonances on the evolution of planetary systems. We find that the semi-major axis distribution is largely determined in a gas disk, while the eccentricity distribution is determined after the disk dissipation. We also find that there may be an optimum disk mass which leads to the observed a-e distribution. Our simulations generate a larger fraction of planetary systems trapped in mean-motion resonances (MMRs) than the observations, indicating that the disk's perturbation to the planetary orbits may be important to explain the observed rate of MMRs. We also find much lower occurrence of planets on retrograde orbits than the current observations of close-in planets suggest.

preprint2009arXiv

Central Limit Theorems for the Energy Density in the Sherrington-Kirkpatrick Model

In this paper we consider central limit theorems for various macroscopic observables in the high temperature region of the Sherrington-Kirkpatrick spin glass model. With a particular focus on obtaining a quenched central limit theorem for the energy density of the system with non-zero external field, we show how to combine the mean field cavity method with Stein's method in the quenched regime. The result for the energy density extends the corresponding result of Comets-Neveu.

preprint2009arXiv

Monte Carlo Simulations of Globular Cluster Evolution. V. Binary Stellar Evolution

We study the dynamical evolution of globular clusters containing primordial binaries, including full single and binary stellar evolution using our Monte Carlo cluster evolution code updated with an adaptation of the single and binary stellar evolution codes SSE/BSE from Hurley et. al (2000, 2002). We describe the modifications we have made to the code. We present several test calculations and comparisons with existing studies to illustrate the validity of the code. We show that our code finds very good agreement with direct N-body simulations including primordial binaries and stellar evolution. We find significant differences in the evolution of the global properties of the simulated clusters using stellar evolution compared to simulations without any stellar evolution. In particular, we find that the mass loss from stellar evolution acts as a significant energy production channel simply by reducing the total gravitational binding energy and can significantly prolong the initial core contraction phase before reaching the binary-burning quasi steady state of the cluster evolution as noticed in Paper IV. We simulate a large grid of clusters varying the initial cluster mass, binary fraction, and concentration and compare properties of the simulated clusters with those of the observed Galactic globular clusters (GGCs). We find that our simulated cluster properties agree well with the observed GGC properties. We explore in some detail qualitatively different clusters in different phases of their evolution, and construct synthetic Hertzprung-Russell diagrams for these clusters.

preprint2008arXiv

Multivariate normal approximation using exchangeable pairs

Since the introduction of Stein's method in the early 1970s, much research has been done in extending and strengthening it; however, there does not exist a version of Stein's original method of exchangeable pairs for multivariate normal approximation. The aim of this article is to fill this void. We present three abstract normal approximation theorems using exchangeable pairs in multivariate contexts, one for situations in which the underlying symmetries are discrete, and real and complex versions of a theorem for situations involving continuous symmetry groups. Our main applications are proofs of the approximate normality of rank $k$ projections of Haar measure on the orthogonal and unitary groups, when $k=o(n)$.

Sourav Chatterjee

What is connected

Connect this record

See the researcher in context

Building this map preview

88 published item(s)

A computationally efficient and mechanically compatible multi-phase-field model applied to coherently stressed three-phase solids

An efficient and quantitative phase-field model for elastically heterogeneous two-phase solids based on a partial rank-one homogenization scheme

A random walk on the Rado graph

Compact Object Modeling in the Globular Cluster 47 Tucanae

Existence of stationary ballistic deposition on the infinite lattice

Globular Clusters UVIT Legacy Survey (GlobULeS) I. FUV-optical Color-Magnitude Diagrams for Eight Globular Clusters

Isomorphisms between random graphs

Matrix completion with data-dependent missingness probabilities

MOSPAT: AutoML based Model Selection and Parameter Tuning for Time Series Anomaly Detection

Polarization correction towards satellite-based QKD without an active feedback

Superconcentration in surface growth

The Yang-Mills heat flow with random distributional initial data

Weak convergence of directed polymers to deterministic KPZ at high temperature

Convergence of deterministic growth models

Intermediate-mass Black Holes from High Massive-star Binary Fractions in Young Star Clusters

Micro-Estimates of Wealth for all Low- and Middle-Income Countries

Modeling Dense Star Clusters in the Milky Way and Beyond with the Cluster Monte Carlo Code

The Observed Rate of Binary Black Hole Mergers can be Entirely Explained by Globular Clusters

A Dynamical Survey of Stellar-Mass Black Holes in 50 Milky Way Globular Clusters

A new coefficient of correlation

Demographics of triple systems in dense star clusters

GW190412 as a Third-Generation Black Hole Merger from a Super Star Cluster

Impact of Local Stiffness on Entropy Driven Microscopic Dynamics of Polythiophene

Modeling Dense Star Clusters in the Milky Way and Beyond with the $\texttt{CMC}$ Cluster Catalog

Speeding up Markov chains with deterministic jumps

Wilson loops in Ising lattice gauge theory

On the Rate of Neutron Star Binary Mergers from Globular Clusters

qkdSim: An experimenter's simulation toolkit for QKD with imperfections, and its performance analysis with a demonstration of the B92 protocol using heralded photon

The Role of "Black Hole Burning" in the Evolution of Dense Star Clusters

A central limit theorem for a new statistic on permutations

A note about the uniform distribution on the intersection of a simplex and a sphere

An introduction to large deviations for random graphs

Binary Black Hole Mergers from Globular Clusters: Masses, Merger Rates, and the Impact of Stellar Evolution

Broadcasting of Quantum Correlations: Possibilities & Impossibilities

Complementarity in Quantum Information Processing Tasks

Concentration inequalities with exchangeable pairs (Ph.D. thesis)

Dynamical Formation of the GW150914 Binary Black Hole

Million-Body Star Cluster Simulations: Comparisons between Monte Carlo and Direct $N$-body

Minimal spanning trees and Stein's method

Nonlinear large deviations

Prediction error of cross-validated Lasso

Sequential Quantum Secret Sharing in a Noisy Environment aided with Weak Measurements

The $1/N$ expansion for SO(N) lattice gauge theory at strong coupling

The leading term of the Yang-Mills free energy

An Overview of Inside-Out Planet Formation

Binary Black Hole Mergers from Globular Clusters: Implications for Advanced LIGO

Black hole mergers and blue stragglers from hierarchical triples formed in globular clusters

High dimensional regression and matrix estimation without tuning parameters

Inside-Out Planet Formation. III. Planet-disk interaction at the dead zone inner boundary

Period Ratio Distribution of Near-Resonant Planets Indicates Planetesimal Scattering

Planetesimal Interactions Can Explain the Mysterious Period Ratios of Small Near-Resonant Planets

A new perspective on least squares under convex constraint

A short survey of Stein's method

Absence of replica symmetry breaking in the random field Ising model

Assumptionless consistency of the Lasso

Fluctuations of the Bose-Einstein condensate

Matrix estimation by Universal Singular Value Thresholding

Pebble Delivery for Inside-Out Planet Formation

Vulcan Planets: Inside-Out Formation of the Innermost Super-Earths

Estimating and understanding exponential random graph models

Inside-Out Planet Formation

Large Deviations for Random Matrices

Mock Observations of Blue Stragglers in Globular Cluster Models

On level sets of Gaussian fields

Stellar Collisions and Blue Straggler Stars in Dense Globular Clusters

Stochastic solutions of the wave equation

Central limit theorem for first-passage percolation time across thin cylinders

Gravitational Slingshot of Young Massive Stars in Orion

Monte Carlo Simulations of Globular Cluster Evolution. VI. The Influence of an Intermediate Mass Black Hole

Planets in Open Clusters Detectable by Kepler

Random Overlap Structures: Properties and Applications to Spin Glasses

The universal relation between scaling exponents in first-passage percolation

Understanding the Dynamical State of Globular Clusters: Core-Collapsed vs Non Core-Collapsed

Nonnormal approximation by Stein's method of exchangeable pairs with application to the Curie--Weiss model