Researcher profile

Sourav Chatterjee

Sourav Chatterjee contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
29works
0followers
25topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

29 published item(s)

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.

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.

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."