Showing 1–50 of 79 results for author: Bhattacharya, K

Comparative Analysis of $k$-essence and Quintessence Scalar Field Models: A Data Analysis Approach

Authors: Saddam Hussain, Sarath Nelleri, Kaushik Bhattacharya

Abstract: We perform a comparative analysis of quintessence and $k$-essence scalar field models in the data analysis perspective. We study the quintessence field with an exponential potential and the $k$-essence field with an inverse square potential in the present work. Before delving into data analysis, we provide a brief perspective on dynamical evolution on both of the models and obtain the stability co… ▽ More We perform a comparative analysis of quintessence and $k$-essence scalar field models in the data analysis perspective. We study the quintessence field with an exponential potential and the $k$-essence field with an inverse square potential in the present work. Before delving into data analysis, we provide a brief perspective on dynamical evolution on both of the models and obtain the stability constraints on the model parameters. We adopt Bayesian inference procedure to estimate the model parameters that best-fit the data. A comprehensive analysis utilizing Observational Hubble data (OHD) and Pantheon+ compilation of Type Ia supernovae (SNIa) shows that $k$-essence model fits the data slightly better than the quintessence model while the evidence of these models in comparison with the $Λ$CDM model is weak. The value of the Hubble constant predicted by both the models is in close agreement with the value obtained by the Planck2018 collaboration assuming the $Λ$CDM model. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 24 pages, 12 figures, 4 tables

Dielectric microwave resonator with large optical apertures for spin-based quantum devices

Authors: Tatsuki Hamamoto, Amit Bhunia, Rupak Kumar Bhattacharya, Hiroki Takahashi, Yuimaru Kubo

Abstract: Towards a spin-based quantum microwave-optical photon transducer, we demonstrate a low-loss dielectric microwave resonator with an internal quality factor of $2.30\times10^4$ while accommodating optical apertures with a diameter of $8\, \mathrm{mm}$. The two seemingly conflicting requirements, high quality factor and large optical apertures, are satisfied thanks to the large dielectric constant of… ▽ More Towards a spin-based quantum microwave-optical photon transducer, we demonstrate a low-loss dielectric microwave resonator with an internal quality factor of $2.30\times10^4$ while accommodating optical apertures with a diameter of $8\, \mathrm{mm}$. The two seemingly conflicting requirements, high quality factor and large optical apertures, are satisfied thanks to the large dielectric constant of rutile ($\mathrm{TiO_2}$). The quality factor is limited by radiation loss, and we confirmed by numerical simulation that this dielectric resonator can achieve a quality factor exceeding $10^6$ by extending the height of the resonator enclosure. Using this resonator, we performed both continuous-wave (cw) and pulse electron spin resonance (ESR) spectroscopy on 2,2-diphenyl-1-picrylhydrazyl (DPPH) crystalline powder and P1 centers in a diamond crystal in a dilution refrigerator. The cw ESR spectroscopy demonstrated high-cooperativity and strong spin-resonator coupling with the DPPH and P1 centers respectively, while the pulse ESR spectroscopy successfully measured longitudinal and transverse relaxation times. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Iterated learning and multiscale modeling of history-dependent architectured metamaterials

Authors: Yupeng Zhang, Kaushik Bhattacharya

Abstract: Neural network based models have emerged as a powerful tool in multiscale modeling of materials. One promising approach is to use a neural network based model, trained using data generated from repeated solution of an expensive small scale model, as a surrogate for the small scale model in application scale simulations. Such approaches have been shown to have the potential accuracy of concurrent m… ▽ More Neural network based models have emerged as a powerful tool in multiscale modeling of materials. One promising approach is to use a neural network based model, trained using data generated from repeated solution of an expensive small scale model, as a surrogate for the small scale model in application scale simulations. Such approaches have been shown to have the potential accuracy of concurrent multiscale methods like FE2, but at the cost comparable to empirical methods like classical constitutive models or parameter passing. A key question is to understand how much and what kind of data is necessary to obtain an accurate surrogate. This paper examines this question for history dependent elastic-plastic behavior of an architected metamaterial modeled as a truss. We introduce an iterative approach where we use the rich arbitrary class of trajectories to train an initial model, but then iteratively update the class of trajectories with those that arise in large scale simulation and use transfer learning to update the model. We show that such an approach converges to a highly accurate surrogate, and one that is transferable. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Residential clustering and mobility of ethnic groups

Authors: Kunal Bhattacharya, Chandreyee Roy, Tuomas Takko, Anna Rotkirch, Kimmo Kaski

Abstract: We studied residential clustering and mobility of ethnic minorities using a theoretical framework based on null models of spatial distributions and movements of populations. Using microdata from population registers we compared the patterns of clustering amongst various socioethnic groups living in and around the capital region of Finland. Using the models we were able to connect the factors influ… ▽ More We studied residential clustering and mobility of ethnic minorities using a theoretical framework based on null models of spatial distributions and movements of populations. Using microdata from population registers we compared the patterns of clustering amongst various socioethnic groups living in and around the capital region of Finland. Using the models we were able to connect the factors influencing intraurban migration to the spatial patterns that have been developed over time. We could also demonstrate the interrelationship of the movement and clustering with fertility. The observed clustering seems to be a combined effect of fertility and the tendency to migrate locally. The models also highlight the importance of factors like proximity to the city-centre, average neighbourhood income, and similarity of socioeconomic profiles. △ Less

Submitted 9 November, 2023; originally announced November 2023.

Comments: 11 pages, 4 figures, supporting information

Optical penetration depth and periodic motion of a photomechanical strip

Authors: Ameneh Maghsoodi, Kaushik Bhattacharya

Abstract: Liquid crystal elastomers (LCEs) containing light-sensitive molecules exhibit large reversible deformation when subjected to illumination. Here, we investigate the role of optical penetration depth on this photomechanical response. We present a model of the photomechanical behavior of photoactive LCE strips under illumination that goes beyond the common assumption of shallow penetration. This mode… ▽ More Liquid crystal elastomers (LCEs) containing light-sensitive molecules exhibit large reversible deformation when subjected to illumination. Here, we investigate the role of optical penetration depth on this photomechanical response. We present a model of the photomechanical behavior of photoactive LCE strips under illumination that goes beyond the common assumption of shallow penetration. This model reveals how the optical penetration depth and the consequent photomechanically induced deformation can depend on the concentration of photoactive molecules, their absorption cross-sections, and the intensity of illumination. Through a series of examples, we show that the penetration depth can quantitatively and qualitatively affect the photomechanical response of a strip. Shallow illumination leads to monotone curvature change while deep penetration can lead to non-monotone response with illumination duration. Further, the flapping behavior (a cyclic wave-like motion) of doubly clamped and buckled strips that has been proposed for locomotion can reverse direction with sufficiently large penetration depth. This opens the possibility of creating wireless light-driven photomechanical actuators and swimmers whose direction of motion can be controlled by light intensity and frequency. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Learning the Effective Adhesive Properties of Heterogeneous Substrates

Authors: Maximo Cravero Baraja, Kaushik Bhattacharya

Abstract: Adhesion is a fundamental phenomenon that plays a role in many engineering and biological applications. This paper concerns the use of machine learning to characterize the effective adhesive properties when a thin film is peeled from a heterogeneous substrate. There has been recent interest in the use of machine learning in multiscale modeling where macroscale constitutive relations are learnt fro… ▽ More Adhesion is a fundamental phenomenon that plays a role in many engineering and biological applications. This paper concerns the use of machine learning to characterize the effective adhesive properties when a thin film is peeled from a heterogeneous substrate. There has been recent interest in the use of machine learning in multiscale modeling where macroscale constitutive relations are learnt from data gathered from repeated solution of the microscale problem. We extend this approach to peeling; this is challenging because peeling from heterogenous substrates is characterized by pinning where the peel front gets stuck at a heterogeneity followed by an abrupt depinning. This results in a heterogeneity dependent critical force and a singular peel force vs. overall peel rate relationship. We propose a neural architecture that is able to accurately predict both the critical peel force and the singular nature of the peel force vs. overall peel rate relationship from the heterogeneous adhesive pattern. Similar issues arise in other free boundary and free discontinuity problems, and the methods we develop are applicable in those contexts. △ Less

Submitted 31 August, 2023; originally announced September 2023.

Semi-classical understanding of flux quantization in superconductors

Authors: Kolahal Bhattacharya

Abstract: Like electric charge, magnetic flux is also quantised. Theoretically, one can show that the flux quantum must be h/e, as observed in the quantum Hall effect. However, in the superconducting systems, the flux quantum is experimentally observed as h/2e. There is no fundamental explanation for the empirical result. In this article, we argue that this phenomenon is fundamentally linked to the nonlocal… ▽ More Like electric charge, magnetic flux is also quantised. Theoretically, one can show that the flux quantum must be h/e, as observed in the quantum Hall effect. However, in the superconducting systems, the flux quantum is experimentally observed as h/2e. There is no fundamental explanation for the empirical result. In this article, we argue that this phenomenon is fundamentally linked to the nonlocality problem of the Aharonov-Bohm effect and present a new semi-classical explanation for the magnetic flux quantum in superconductivity. This work will also show why the flux quantum should be h/e in the case of the quantum Hall effect. △ Less

Submitted 19 August, 2023; originally announced August 2023.

Comments: 8 pages

Journal ref: Quantum Studies: Mathematics and Foundations, Vol 11, Issue 2, 2024

Semi-classical description of electrostatics and quantization of electric charge

Authors: Kolahal Bhattacharya

Abstract: In this work, we present an explanation of the electric charge quantization based on a semi-classical model of electrostatic fields. We claim that in electrostatics, an electric charge must be equal to a rational multiple of the elementary charge of an electron. However, the charge is quantized if the system has certain boundary conditions that force the wavefunction representing an electric field… ▽ More In this work, we present an explanation of the electric charge quantization based on a semi-classical model of electrostatic fields. We claim that in electrostatics, an electric charge must be equal to a rational multiple of the elementary charge of an electron. However, the charge is quantized if the system has certain boundary conditions that force the wavefunction representing an electric field to vanish at specific surfaces. Next, we develop the corresponding model for the electric displacement vector. It is demonstrated that a number of classical results, e.g. bending of field lines at the interface of two dielectric media, method of images, etc. are all consistent with the predictions of this model. We also present the possible form of Gauss's law or (Poisson's equation), to find the wavefunctions of the field from a source charge distribution, in this model. △ Less

Submitted 4 July, 2023; originally announced July 2023.

Comments: 19 pages, 1 figure, accepted for publication

Journal ref: Physica Scripta, 2023

Modelling exposure between populations using networks of mobility during Covid-19

Authors: Tuomas Takko, Kunal Bhattacharya, Kimmo Kaski

Abstract: The use of mobile phone call detail records and device location data for the calling patterns, movements, and social contacts of individuals, has proven to be valuable for devising models and understanding of their mobility and behaviour patterns. In this study we investigate weighted exposure-networks of human daily activities in the capital region of Finland as a proxy for contacts between posta… ▽ More The use of mobile phone call detail records and device location data for the calling patterns, movements, and social contacts of individuals, has proven to be valuable for devising models and understanding of their mobility and behaviour patterns. In this study we investigate weighted exposure-networks of human daily activities in the capital region of Finland as a proxy for contacts between postal code areas during the pre-pandemic year 2019 and pandemic years 2020, 2021 and early 2022. We investigate the suitability of gravity and radiation type models for reconstructing the exposure-networks based on geo-spatial and population mobility information. For this we use a mobile phone dataset of aggregated daily visits from a postal code area to cellphone grid locations, and treat it as a bipartite network to create weighted one mode projections using a weighted co-occurrence function. We fit a gravitation model and a radiation model to the averaged weekly and yearly projection networks with geo-spatial and socioeconomic variables of the postal code areas and their populations. We also consider an extended gravity type model comprising of additional postal area information such as distance via public transportation and population density. The results show that the co-occurrence of human activities, or exposure, between postal code areas follows both the gravity and radiation type interactions, once fitted to the empirical network. The effects of the pandemic beginning in 2020 can be observed as a decrease of the overall activity as well as of the exposure of the projected networks. In general, the results show that the postal code level networks changed to be more proximity weighted after the pandemic began, following the government imposed non-pharmaceutical interventions, with differences based on the geo-spatial and socioeconomic structure of the areas. △ Less

Submitted 17 April, 2023; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: Supplementary information appended to manuscript. Contents: 15+3 pages, 5+3 figures, 2+1 tables. Submitted to Frontiers in Physics, Social Physics

Light-induced swirling and locomotion

Authors: Ameneh Maghsoodi, Kaushik Bhattacharya

Abstract: Photomechanical liquid crystal elastomers (LCEs) are responsive polymers that can convert light directly into mechanical deformation. This unique feature makes these materials an attractive candidate for soft actuators capable of remote and multi-mode actuation. In this work, we propose a three-dimensional multi-scale model of the nonlinear and nonlocal dynamics of fibers of photomechanical LCEs u… ▽ More Photomechanical liquid crystal elastomers (LCEs) are responsive polymers that can convert light directly into mechanical deformation. This unique feature makes these materials an attractive candidate for soft actuators capable of remote and multi-mode actuation. In this work, we propose a three-dimensional multi-scale model of the nonlinear and nonlocal dynamics of fibers of photomechanical LCEs under illumination. We use the model to show that a pre-stressed helix-like fibers immersed in a fluid can undergo a periodic whirling motion under steady illumination. We analyze the photo-driven spatiotemporal pattern and stability of the whirling deformation, and provide a parametric study. Unlike previous work on photo-driven periodic motion, this whirling motion does not exploit instabilities in the form of snap-through phenomena, or gravity as in rolling. We then show that such motion can be exploited in developing remote controlled bio-inspired microswimmers and novel micromixers. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Aggregation of self-propelled particles with sensitivity to local order

Authors: Kunal Bhattacharya, Abhijit Chakraborty

Abstract: We study a system of self-propelled particles (SPPs) in which individual particles are allowed to switch between a fast aligning and a slow nonaligning state depending upon the degree of the alignment in the neighborhood. The switching is modeled using a threshold for the local order parameter. This additional attribute gives rise to a mixed phase, in contrast to the ordered phases found in clean… ▽ More We study a system of self-propelled particles (SPPs) in which individual particles are allowed to switch between a fast aligning and a slow nonaligning state depending upon the degree of the alignment in the neighborhood. The switching is modeled using a threshold for the local order parameter. This additional attribute gives rise to a mixed phase, in contrast to the ordered phases found in clean SPP systems. As the threshold is increased from zero, we find the sudden appearance of clusters of nonaligners. Clusters of nonaligners coexist with moving clusters of aligners with continual coalescence and fragmentation. The behavior of the system with respect to the clustering of nonaligners appears to be very different for values of low and high global densities. In the low density regime, for an optimal value of the threshold, the largest cluster of nonaligners grows in size up to a maximum that varies logarithmically with the total number of particles. However, on further increasing the threshold the size decreases. In contrast, for the high density regime, an initial abrupt rise is followed by the appearance of a giant cluster of nonaligners. The latter growth can be characterized as a continuous percolation transition. In addition, we find that the speed differences between aligners and nonaligners is necessary for the segregation of aligners and nonaligners. △ Less

Submitted 17 April, 2022; originally announced April 2022.

Comments: 11 pages, 7 figures

Journal ref: Phys. Rev. E 105, 044124 (2022)

A complement to the scalar wave theory of light

Authors: Kolahal Bhattacharya

Abstract: In this paper, we discuss how the concepts of Hamiltonian optics are internally connected to the scalar wave theory of light rays. It is shown that the solutions of the reduced wave equation are similar to Huygen's wavelets, and they can be used to understand Snell's law of refraction. This model can also be used to derive the coefficient of reflection consistent with Fresnel's equation for $s$-po… ▽ More In this paper, we discuss how the concepts of Hamiltonian optics are internally connected to the scalar wave theory of light rays. It is shown that the solutions of the reduced wave equation are similar to Huygen's wavelets, and they can be used to understand Snell's law of refraction. This model can also be used to derive the coefficient of reflection consistent with Fresnel's equation for $s$-polarized light. △ Less

Submitted 23 March, 2022; originally announced March 2022.

Comments: 10 pages, 3 figures, 1 table

Journal ref: European Journal of Physics, Volume 43, No. 3, p 035305. (2022). (IOP Publishing)

Turnover in close friendships: age and gender differences

Authors: Chandreyee Roy, Kunal Bhattacharya, Robin I. M. Dunbar, Kimmo Kaski

Abstract: Humans are social animals and the interpersonal bonds formed between them are crucial for their development and well being in a society. These relationships are usually structured into several layers (Dunbar's layers of friendship) depending on their significance in an individual's life with closest friends and family being the most important ones taking major part of their time and communication… ▽ More Humans are social animals and the interpersonal bonds formed between them are crucial for their development and well being in a society. These relationships are usually structured into several layers (Dunbar's layers of friendship) depending on their significance in an individual's life with closest friends and family being the most important ones taking major part of their time and communication effort. However, we have little idea how the initiation and termination of these relationships occurs across the lifespan. To explore this, we analyse a national cellphone database to determine how and when changes in close relationships occur in the two genders. In general, membership of this inner circle of intimate relationships is extremely stable, at least over a three-year period. However, around 1-4% of alters change every year, with the rate of change being higher among 17-21 year olds than older adults. Young adult females terminate more of their opposite-gender relationships, while older males are more persistent in trying to maintain relationships in decline. These results emphasise the variability in relationship dynamics across age and gender, and remind us that individual differences play an important role in the structure of social networks. Overall, our study provides a holistic understanding of the dynamic nature of relationships during the life-course of humans. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 11 pages, 6 figures

Unexplored aspects of a variational principle in electrostatics

Authors: Kolahal Bhattacharya, Debapriyo Syam

Abstract: The electrostatic field magnitude can play a role in a variational principle similar to the role of the index of refraction in geometrical optics, allowing the determination of the electric field lines. This was identified in the context of the grounded conducting sphere image problem. Assuming the knowledge of the magnitude of the electric field, in this paper we validate this principle for the g… ▽ More The electrostatic field magnitude can play a role in a variational principle similar to the role of the index of refraction in geometrical optics, allowing the determination of the electric field lines. This was identified in the context of the grounded conducting sphere image problem. Assuming the knowledge of the magnitude of the electric field, in this paper we validate this principle for the general case and explore some of its consequences. △ Less

Submitted 8 December, 2021; originally announced January 2022.

Comments: 7 pages (single column), 2 figures; Accepted in American Journal of Physics (AAPT) on 7th December 2021

Fracture Diodes: Directional asymmetry of fracture toughness

Authors: N. R. Brodnik, S. Brach, C. M. Long, G. Ravichandran, B. Bourdin, K. T. Faber, K. Bhattacharya

Abstract: Toughness describes the ability of a material to resist fracture or crack propagation. It is demonstrated here that fracture toughness of a material can be asymmetric, i.e., the resistance of a medium to a crack propagating from right to left can be significantly different from that to a crack propagating from left to right. Such asymmetry is unknown in natural materials, but we show that it can b… ▽ More Toughness describes the ability of a material to resist fracture or crack propagation. It is demonstrated here that fracture toughness of a material can be asymmetric, i.e., the resistance of a medium to a crack propagating from right to left can be significantly different from that to a crack propagating from left to right. Such asymmetry is unknown in natural materials, but we show that it can be built into artificial materials through the proper control of microstructure. This paves the way for control of crack paths and direction, where fracture -- when unavoidable -- can be guided through pre-designed paths to minimize loss of critical components. △ Less

Submitted 9 October, 2021; originally announced October 2021.

Journal ref: Physical Review Letters, 126(2), 025503 (2021)

The Life and Science of Thanu Padmanabhan

Authors: Jasjeet Singh Bagla, Krishnakanta Bhattacharya, Sumanta Chakraborty, Sunu Engineer, Valerio Faraoni, Sanved Kolekar, Dawood Kothawala, Kinjalk Lochan, Sujoy Modak, V. Parameswaran Nair, Aseem Paranjape, Krishnamohan Parattu, Sarada G. Rajeev, Bibhas Ranjan Majhi, Tirthankar Roy Choudhury, Mohammad Sami, Sudipta Sarkar, Sandipan Sengupta, T. R. Seshadri, S. Shankaranarayanan, Suprit Singh, Tejinder P. Singh, L. Sriramkumar, Urjit Yajnik

Abstract: Thanu Padmanabhan was a renowned Indian theoretical physicist known for his research in general relativity, cosmology, and quantum gravity. In an extraordinary career spanning forty-two years, he published more than three hundred research articles, wrote ten highly successful technical and popular books, and mentored nearly thirty graduate students and post-doctoral fellows. He is best known for h… ▽ More Thanu Padmanabhan was a renowned Indian theoretical physicist known for his research in general relativity, cosmology, and quantum gravity. In an extraordinary career spanning forty-two years, he published more than three hundred research articles, wrote ten highly successful technical and popular books, and mentored nearly thirty graduate students and post-doctoral fellows. He is best known for his deep work investigating gravitation as an emergent thermodynamic phenomenon. He was an outstanding teacher, and an indefatigable populariser of science, who travelled very widely to motivate and inspire young students. Paddy, as he was affectionately known, was also a close friend to his students and collaborators, treating them as part of his extended academic family. On September 17, 2021 Paddy passed away very unexpectedly, at the age of sixty-four and at the height of his research career, while serving as a Distinguished Professor at the Inter-University Centre for Astronomy and Astrophysics, Pune. His untimely demise has come as a shock to his family and friends and colleagues. In this article, several of them have come together to pay their tributes and share their fond memories of Paddy. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: In Memoriam of Prof. T. Padmanabhan (1957-2021)

Demystifying the nonlocality problem in Aharonov-Bohm effect

Authors: Kolahal Bhattacharya

Abstract: In this paper, we present a novel semi-classical theory of the electrostatic and magnetostatic fields and explain the nonlocality problem in the context of the Aharonov-Bohm effect [1]. Specifically, we show that the electrostatic and the magnetostatic fields possess a quantum nature that manifests if certain conditions are met. In particular, the wave amplitudes of the fields are seen to exist ev… ▽ More In this paper, we present a novel semi-classical theory of the electrostatic and magnetostatic fields and explain the nonlocality problem in the context of the Aharonov-Bohm effect [1]. Specifically, we show that the electrostatic and the magnetostatic fields possess a quantum nature that manifests if certain conditions are met. In particular, the wave amplitudes of the fields are seen to exist even in the regions where the classical fields vanish and they operate on the electron wave functions locally as unitary phases. This formulation also sheds light on the quantisation of electric charges and magnetic flux. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: 15 pages, 1 figure

Journal ref: Physica Scripta, Published 26 May 2021

Tuning acoustic impedance in load-bearing structures

Authors: Sai Sharan Injeti, Paolo Celli, Kaushik Bhattacharya, Chiara Daraio

Abstract: Acoustic transparency is the capability of a medium to transmit mechanical waves to adjacent media, without scattering. This characteristic can be achieved by carefully engineering the acoustic impedance of the medium -- a combination of wave speed and density, to match that of the surroundings. Owing to the strong correlation between acoustic wave speed and static stiffness, it is challenging to… ▽ More Acoustic transparency is the capability of a medium to transmit mechanical waves to adjacent media, without scattering. This characteristic can be achieved by carefully engineering the acoustic impedance of the medium -- a combination of wave speed and density, to match that of the surroundings. Owing to the strong correlation between acoustic wave speed and static stiffness, it is challenging to design acoustically transparent materials in a fluid, while maintaining their high structural rigidity. In this work, we propose a method to design architected lattices with independent control of the elastic wave speed at a chosen frequency, the mass density, and the static stiffness, along a chosen loading direction. We provide a sensitivity analysis to optimize these properties with respect to design parameters of the structure, that include localized masses at specific positions. We demonstrate the method on five different periodic, three dimensional lattices, to calculate bounds on the longitudinal wave speed as a function of their density and stiffness. We then perform experiments on 3-D printed structures, to validate our numerical simulations. The tools developed in this work can be used to design lightweight and stiff materials with optimized acoustic impedance for a plethora of applications, including ultrasound imaging, wave filtering and waveguiding. △ Less

Submitted 19 June, 2021; originally announced June 2021.

Human-agent coordination in a group formation game

Authors: Tuomas Takko, Kunal Bhattacharya, Daniel Monsivais, Kimmo Kaski

Abstract: Coordination and cooperation between humans and autonomous agents in cooperative games raises interesting questions of human decision making and behaviour changes. Here we report our findings from a group formation game in a small-world network of different mixes of human and agent players, aiming to achieve connected clusters of the same colour by swapping places with neighbouring players using n… ▽ More Coordination and cooperation between humans and autonomous agents in cooperative games raises interesting questions of human decision making and behaviour changes. Here we report our findings from a group formation game in a small-world network of different mixes of human and agent players, aiming to achieve connected clusters of the same colour by swapping places with neighbouring players using non-overlapping information. In the experiments the human players are incentivized by rewarding to prioritize their own cluster while the model of agents' decision making is derived from our previous experiment of purely cooperative game between human players. The experiments were performed by grouping the players in three different setups to investigate the overall effect of having cooperative autonomous agents within teams. We observe that the change in the behavior of human subjects adjusts to playing with autonomous agents by being less risk averse, while keeping the overall performance efficient by splitting the behaviour into selfish and cooperative in the two actions performed during the rounds of the game. Moreover, results from two hybrid human-agent setups suggest that the group composition affects the evolution of clusters. Our findings indicate that in purely or lesser cooperative settings, providing more control to humans could help in maximizing the overall performance of hybrid systems. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: Chosen to be published in Scientific Reports on 24.05.2021 with DOI: 10.1038/s41598-021-90123-8

ACM Class: I.6.0; G.3; J.4

Machine-learned prediction of the electronic fields in a crystal

Authors: Ying Shi Teh, Swarnava Ghosh, Kaushik Bhattacharya

Abstract: We propose an approach for exploiting machine learning to approximate electronic fields in crystalline solids subjected to deformation. Strain engineering is emerging as a widely used method for tuning the properties of materials, and this requires repeated density functional theory calculations of the unit cell subjected to strain. Repeated unit cell calculations are also required for multi-resol… ▽ More We propose an approach for exploiting machine learning to approximate electronic fields in crystalline solids subjected to deformation. Strain engineering is emerging as a widely used method for tuning the properties of materials, and this requires repeated density functional theory calculations of the unit cell subjected to strain. Repeated unit cell calculations are also required for multi-resolution studies of defects in crystalline solids. We propose an approach that uses data from such calculations to train a carefully architected machine learning approximation. We demonstrate the approach on magnesium, a promising light-weight structural material: we show that we can predict the energy and electronic fields to the level of chemical accuracy, and even capture lattice instabilities. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Hierarchical multiscale quantification of material uncertainty

Authors: Burigede Liu, Xingsheng Sun, Kaushik Bhattacharya, Michael Ortiz

Abstract: The macroscopic behavior of many materials is complex and the end result of mechanisms that operate across a broad range of disparate scales. An imperfect knowledge of material behavior across scales is a source of epistemic uncertainty of the overall material behavior. However, assessing this uncertainty is difficult due to the complex nature of material response and the prohibitive computational… ▽ More The macroscopic behavior of many materials is complex and the end result of mechanisms that operate across a broad range of disparate scales. An imperfect knowledge of material behavior across scales is a source of epistemic uncertainty of the overall material behavior. However, assessing this uncertainty is difficult due to the complex nature of material response and the prohibitive computational cost of integral calculations. In this paper, we exploit the multiscale and hierarchical nature of material response to develop an approach to quantify the overall uncertainty of material response without the need for integral calculations. Specifically, we bound the uncertainty at each scale and then combine the partial uncertainties in a way that provides a bound on the overall or integral uncertainty. The bound provides a conservative estimate on the uncertainty. Importantly, this approach does not require integral calculations that are prohibitively expensive. We demonstrate the framework on the problem of ballistic impact of a polycrystalline magnesium plate. Magnesium and its alloys are of current interest as promising light-weight structural and protective materials. Finally, we remark that the approach can also be used to study the sensitivity of the overall response to particular mechanisms at lower scales in a materials-by-design approach. △ Less

Submitted 4 February, 2021; originally announced February 2021.

Spectral quadrature for the first principles study of crystal defects: Application to magnesium

Authors: Swarnava Ghosh, Kaushik Bhattacharya

Abstract: We present an accurate and efficient finite-difference formulation and parallel implementation of Kohn-Sham Density (Operator) Functional Theory (DFT) for non periodic systems embedded in a bulk environment. Specifically, employing non-local pseudopotentials, local reformulation of electrostatics, and truncation of the spatial Kohn-Sham Hamiltonian, and the Linear Scaling Spectral Quadrature metho… ▽ More We present an accurate and efficient finite-difference formulation and parallel implementation of Kohn-Sham Density (Operator) Functional Theory (DFT) for non periodic systems embedded in a bulk environment. Specifically, employing non-local pseudopotentials, local reformulation of electrostatics, and truncation of the spatial Kohn-Sham Hamiltonian, and the Linear Scaling Spectral Quadrature method to solve for the pointwise electronic fields in real-space and the non-local component of the atomic force, we develop a parallel finite difference framework suitable for distributed memory computing architectures to simulate non-periodic systems embedded in a bulk environment. Choosing examples from magnesium-aluminum alloys, we first demonstrate the convergence of energies and forces with respect to spectral quadrature polynomial order, and the width of the spatially truncated Hamiltonian. Next, we demonstrate the parallel scaling of our framework, and show that the computation time and memory scale linearly with respect to the number of atoms. Next, we use the developed framework to simulate isolated point defects and their interactions in magnesium-aluminum alloys. Our findings conclude that the binding energies of divacancies, Al solute-vacancy and two Al solute atoms are anisotropic and are dependent on cell size. Furthermore, the binding is favorable in all three cases. △ Less

Submitted 26 November, 2020; originally announced November 2020.

Understanding the morphotropic phase boundary of perovskite solid solutions as a frustrated state

Authors: Ying Shi Teh, Jiangyu Li, Kaushik Bhattacharya

Abstract: Perovskite solid solutions that have a chemical composition A(C$_x$D$_{1-x})$O$_3$ with transition metals C and D substitutionally occupying the B site of a perovskite lattice are attractive in various applications for their dielectric, piezoelectric and other properties. A remarkable feature of these solid solutions is the \emph{morphotropic phase boundary} (MPB), the composition across which the… ▽ More Perovskite solid solutions that have a chemical composition A(C$_x$D$_{1-x})$O$_3$ with transition metals C and D substitutionally occupying the B site of a perovskite lattice are attractive in various applications for their dielectric, piezoelectric and other properties. A remarkable feature of these solid solutions is the \emph{morphotropic phase boundary} (MPB), the composition across which the crystal symmetry changes. Critically, it has long been observed that the dielectric and piezoelectric as well as the ability to pole a ceramic increases dramatically at the MPB. While this has motivated much study of perovskite MPBs, a number of important questions about the role of disorder remain unanswered. We address these questions using a new approach based on the random-field Ising model with long-range interactions that incorporates the basic elements of the physics at the meso-scale. We show that the MPB emerges naturally in this approach as a frustrated state where stability is exchanged between two well-defined phases. Specifically, long-range interactions suppress the disorder at compositions away from MPB but are unable to do so when there is an exchange of stability. Further, the approach also predicts a number of experimentally observed features like the fragmented domain patterns and superior ability to pole at the MPB. The insights from this model also suggest the possibility of entirely new materials with strong ferroelectric-ferromagnetic coupling using an MPB. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Journal ref: Phys. Rev. B 103, 144201 (2021)

Internal migration and mobile communication patterns among pairs with strong ties

Authors: Mikaela Irene D. Fudolig, Daniel Monsivais, Kunal Bhattacharya, Hang-Hyun Jo, Kimmo Kaski

Abstract: Using large-scale call detail records of anonymised mobile phone service subscribers with demographic and location information, we investigate how a long-distance residential move within the country affects the mobile communication patterns between an ego who moved and a frequently called alter who did not move. By using clustering methods in analysing the call frequency time series, we find that… ▽ More Using large-scale call detail records of anonymised mobile phone service subscribers with demographic and location information, we investigate how a long-distance residential move within the country affects the mobile communication patterns between an ego who moved and a frequently called alter who did not move. By using clustering methods in analysing the call frequency time series, we find that such ego-alter pairs are grouped into two clusters, those with the call frequency increasing and those with the call frequency decreasing after the move of the ego. This indicates that such residential moves are correlated with a change in the communication pattern soon after moving. We find that the pre-move calling behaviour is a relevant predictor for the post-move calling behaviour. While demographic and location information can help in predicting whether the call frequency will rise or decay, they are not relevant in predicting the actual call frequency volume. We also note that at four months after the move, most of these close pairs maintain contact, even if the call frequency is decreased. △ Less

Submitted 5 April, 2021; v1 submitted 1 September, 2020; originally announced September 2020.

Comments: published version in EPJ Data Science

Journal ref: EPJ Data Science 10, 16 (2021)

Can magnetic field be used to reduce cosmic charged particles background to low energy particle detectors?

Authors: Kolahal Bhattacharya

Abstract: The possibility to reduce the background due to cosmic ray charged particles by the use of magnetic field in the ground based low energy particle detectors is explored. The degree of reduction of cosmic rays as a function of the magnetic field strength and its depth is quantified. The possibility to reduce the background due to cosmic ray charged particles by the use of magnetic field in the ground based low energy particle detectors is explored. The degree of reduction of cosmic rays as a function of the magnetic field strength and its depth is quantified. △ Less

Submitted 13 June, 2020; originally announced June 2020.

Comments: 14 pages, 17 figures, 6 tables

Scaling the training of particle classification on simulated MicroBooNE events to multiple GPUs

Authors: Alex Hagen, Eric Church, Jan Strube, Kolahal Bhattacharya, Vinay Amatya

Abstract: Measurements in Liquid Argon Time Projection Chamber (LArTPC) neutrino detectors, such as the MicroBooNE detector at Fermilab, feature large, high fidelity event images. Deep learning techniques have been extremely successful in classification tasks of photographs, but their application to LArTPC event images is challenging, due to the large size of the events. Events in these detectors are typica… ▽ More Measurements in Liquid Argon Time Projection Chamber (LArTPC) neutrino detectors, such as the MicroBooNE detector at Fermilab, feature large, high fidelity event images. Deep learning techniques have been extremely successful in classification tasks of photographs, but their application to LArTPC event images is challenging, due to the large size of the events. Events in these detectors are typically two orders of magnitude larger than images found in classical challenges, like recognition of handwritten digits contained in the MNIST database or object recognition in the ImageNet database. Ideally, training would occur on many instances of the entire event data, instead of many instances of cropped regions of interest from the event data. However, such efforts lead to extremely long training cycles, which slow down the exploration of new network architectures and hyperparameter scans to improve the classification performance. We present studies of scaling a LArTPC classification problem on multiple architectures, spanning multiple nodes. The studies are carried out on simulated events in the MicroBooNE detector. We emphasize that it is beyond the scope of this study to optimize networks or extract the physics from any results here. Institutional computing at Pacific Northwest National Laboratory and the SummitDev machine at Oak Ridge National Laboratory's Leadership Computing Facility have been used. To our knowledge, this is the first use of state-of-the-art Convolutional Neural Networks for particle physics and their attendant compute techniques onto the DOE Leadership Class Facilities. We expect benefits to accrue particularly to the Deep Underground Neutrino Experiment (DUNE) LArTPC program, the flagship US High Energy Physics (HEP) program for the coming decades. △ Less

Submitted 17 April, 2020; originally announced April 2020.

Comments: 6 pages, 4 figures, Accepted for publication in Journal of Physics: Conference Series - Proceedings of the 19th International Workshop on Advanced Computing and Analysis Techniques in Physics Research

Data analytics approach to predict the hardness of the copper matrix composites

Authors: Somesh Kr. Bhattacharya, Ryoji Sahara, Dusan Bozic, Jovana Ruzic

Abstract: Copper matrix composite materials have exhibited a high potential in applications where excellent conductivity and mechanical properties are required. In this study, the machine learning models have been applied to predict the hardness of the copper matrix composite materials produced via powder metallurgy technique. Two particular composites were considered in this work. From experiments, we extr… ▽ More Copper matrix composite materials have exhibited a high potential in applications where excellent conductivity and mechanical properties are required. In this study, the machine learning models have been applied to predict the hardness of the copper matrix composite materials produced via powder metallurgy technique. Two particular composites were considered in this work. From experiments, we extracted the independent variables (features) like the milling time (MT, Hours), dislocation density (DD, 1/m2 ), average particle size (PS,nm), density (gm/cm3 ) and yield stress (MPa) while the Vickers Hardness (MPa) was used as the dependent variable. Feature selection was performed by calculation the Pearson correlation coefficient (PCC) between the independent and dependent variables. We employed six different machine learning regression models to predict the hardness for the two matrix composites. △ Less

Submitted 24 September, 2020; v1 submitted 24 February, 2020; originally announced February 2020.

Comments: 9 pages, 3 figures, 2 tables

Fluid-gravity correspondence in the scalar-tensor theory of gravity: (in)equivalence of Einstein and Jordan frames

Authors: Krishnakanta Bhattacharya, Bibhas Ranjan Majhi, Douglas Singleton

Abstract: The duality of gravitational dynamics (projected on a null hypersurface) and of fluid dynamics is investigated for the scalar tensor (ST) theory of gravity. The description of ST gravity, in both Einstein and Jordan frames, is analyzed from fluid-gravity viewpoint. In the Einstein frame the dynamical equation for the metric leads to the Damour-Navier-Stokes (DNS) equation with an external forcing… ▽ More The duality of gravitational dynamics (projected on a null hypersurface) and of fluid dynamics is investigated for the scalar tensor (ST) theory of gravity. The description of ST gravity, in both Einstein and Jordan frames, is analyzed from fluid-gravity viewpoint. In the Einstein frame the dynamical equation for the metric leads to the Damour-Navier-Stokes (DNS) equation with an external forcing term, coming from the scalar field in ST gravity. In the Jordan frame the situation is more subtle. We observe that finding the DNS equation in this frame can lead to two pictures. In one picture, the usual DNS equation is modified by a Coriolis-like force term, which originates completely from the presence of a non-minimally coupled scalar field ($φ$) on the gravity side. Moreover, the identified fluid variables are no longer conformally equivalent with those in the Einstein frame. However, this picture is consistent with the saturation of Kovtun-Son-Starinets (KSS) bound. In the other picture, we find the standard DNS equation (i.e. without the Coriolis-like force), with the fluid variables conformally equivalent with those in Einstein frame. But, the second picture, may not agree with the KSS bound for some values of $φ$. We conclude by rewriting the Raychaudhuri equation and the tidal force equation in terms of the relevant parameters to demonstrate how the expansion scalar and the shear-tensor evolve in the spacetime. Although, the area law of entropy is broken in ST gravity, we show that the rewritten form of Raychaudhuri's equation correctly results in the generalized second law of black hole thermodynamics. △ Less

Submitted 3 June, 2020; v1 submitted 11 February, 2020; originally announced February 2020.

Comments: New comments added, to appear in JHEP

Journal ref: JHEP07(2020)018

Supervised learning of photoelectron counting in scintillator-based dark matter experiments

Authors: Kolahal Bhattacharya, Christopher Jackson

Abstract: Many scintillator based detectors employ a set of photomultiplier tubes (PMT) to observe the scintillation light from potential signal and background events. It is important to be able to count the number of photoelectrons (PE) in the pulses observed in the PMTs, because the position and energy reconstruction of the events is directly related to how well the spatial distribution of the PEs in the… ▽ More Many scintillator based detectors employ a set of photomultiplier tubes (PMT) to observe the scintillation light from potential signal and background events. It is important to be able to count the number of photoelectrons (PE) in the pulses observed in the PMTs, because the position and energy reconstruction of the events is directly related to how well the spatial distribution of the PEs in the PMTs as well as their total number might be measured. This task is challenging for fast scintillators, since the PEs often overlap each other in time. Standard Bayesian statistics methods are often used and this has been the method employed in analyzing the data from liquid argon experiments such as MiniCLEAN and DEAP. In this work, we show that for the MiniCLEAN detector it is possible to use a multi-layer perceptron to learn the number of PEs using only raw pulse features with better accuracy and precision than existing methods. This can even help to perform position reconstruction with better accuracy and precision, at least in some generic cases. △ Less

Submitted 29 October, 2019; v1 submitted 30 September, 2019; originally announced October 2019.

Comments: Talk presented at the 2019 Meeting of the Division of Particles and Fields of the American Physical Society (DPF2019), July 29 - August 2, 2019, Northeastern University, Boston, C1907293

Dynamical properties of hierarchical networks of Van Der Pol oscillators

Authors: Daniel Monsivais, Kunal Bhattacharya, Rafael A. Barrio, Philip K. Maini, Kimmo K. Kaski

Abstract: Oscillator networks found in social and biological systems are characterized by the presence of wide ranges of coupling strengths and complex organization. Yet robustness and synchronization of oscillations are found to emerge on macro-scales that eventually become key to the functioning of these systems. In order to model this kind of dynamics observed, for example, in systems of circadian oscill… ▽ More Oscillator networks found in social and biological systems are characterized by the presence of wide ranges of coupling strengths and complex organization. Yet robustness and synchronization of oscillations are found to emerge on macro-scales that eventually become key to the functioning of these systems. In order to model this kind of dynamics observed, for example, in systems of circadian oscillators, we study networks of Van der Pol oscillators that are connected with hierarchical couplings. For each isolated oscillator we assume the same fundamental frequency. Using numerical simulations, we show that the coupled system goes to a phase-locked state, with both phase and frequency being the same for every oscillator at each level of the hierarchy. The observed frequency at each level of the hierarchy changes, reaching an asymptotic lowest value at the uppermost level. Notably, the asymptotic frequency can be tuned to any value below the fundamental frequency of an uncoupled Van der Pol oscillator. We compare the numerical results with those of an approximate analytic solution and find them to be in qualitative agreement. △ Less

Submitted 16 December, 2019; v1 submitted 10 September, 2019; originally announced September 2019.

Comments: Updated version of the manuscript submitted to PRE

Link-centric analysis of variation by demographics in mobile phone communication patterns

Authors: Mikaela Irene D. Fudolig, Kunal Bhattacharya, Daniel Monsivais, Hang-Hyun Jo, Kimmo Kaski

Abstract: We present a link-centric approach to study variation in the mobile phone communication patterns of individuals. Unlike most previous research on call detail records that focused on the variation of phone usage across individual users, we examine how the calling and texting patterns obtained from call detail records vary among pairs of users and how these patterns are affected by the nature of rel… ▽ More We present a link-centric approach to study variation in the mobile phone communication patterns of individuals. Unlike most previous research on call detail records that focused on the variation of phone usage across individual users, we examine how the calling and texting patterns obtained from call detail records vary among pairs of users and how these patterns are affected by the nature of relationships between users. To demonstrate this link-centric perspective, we extract factors that contribute to the variation in the mobile phone communication patterns and predict demographics-related quantities for pairs of users. The time of day and the channel of communication (calls or texts) are found to explain most of the variance among pairs that frequently call each other. Furthermore, we find that this variation can be used to predict the relationship between the pairs of users, as inferred from their age and gender, as well as the age of the younger user in a pair. From the classifier performance across different age and gender groups as well as the inherent class overlap suggested by the estimate of the bounds of the Bayes error, we gain insights into the similarity and differences of communication patterns across different relationships. △ Less

Submitted 16 December, 2019; v1 submitted 31 July, 2019; originally announced July 2019.

Journal ref: PLoS ONE 15(1) (2020): e0227037

Large scale ab-initio simulations of dislocations

Authors: Mauricio Ponga, Kaushik Bhattacharya, Michael Ortiz

Abstract: We present a novel methodology to compute relaxed dislocations core configurations, and their energies in crystalline metallic materials using large-scale \emph{ab-intio} simulations. The approach is based on MacroDFT, a coarse-grained density functional theory method that accurately computes the electronic structure but with sub-linear scaling resulting in a tremendous reduction in cost. Due to i… ▽ More We present a novel methodology to compute relaxed dislocations core configurations, and their energies in crystalline metallic materials using large-scale \emph{ab-intio} simulations. The approach is based on MacroDFT, a coarse-grained density functional theory method that accurately computes the electronic structure but with sub-linear scaling resulting in a tremendous reduction in cost. Due to its implementation in \emph{real-space}, MacroDFT has the ability to harness petascale resources to study materials and alloys through accurate \emph{ab-initio} calculations. Thus, the proposed methodology can be used to investigate dislocation cores and other defects where long range elastic defects play an important role, such as in dislocation cores, grain boundaries and near precipitates in crystalline materials. We demonstrate the method by computing the relaxed dislocation cores in prismatic dislocation loops and dislocation segments in magnesium (Mg). We also study the interaction energy with a line of Aluminum (Al) solutes. Our simulations elucidate the essential coupling between the quantum mechanical aspects of the dislocation core and the long range elastic fields that they generate. In particular, our quantum mechanical simulations are able to describe the logarithmic divergence of the energy in the far field as is known from classical elastic theory. In order to reach such scaling, the number of atoms in the simulation cell has to be exceedingly large, and cannot be achieved with the state-of-the-art density functional theory implementations. △ Less

Submitted 4 July, 2019; originally announced July 2019.

Triplet Lifetime in Gaseous Argon

Authors: Michael Akashi-Ronquest, Amanda Bacon, Christopher Benson, Kolahal Bhattacharya, Thomas Caldwell, Joseph A. Formaggio, Dan Gastler, Brianna Grado-White, Jeff Griego, Michael Gold, Andrew Hime, Christopher M. Jackson, Stephen Jaditz, Chris Kachulis, Edward Kearns, Joshua R. Klein, Antonio Ledesma, Steve Linden, Frank Lopez, Sean MacMullin, Andrew Mastbaum, Jocelyn Monroe, James Nikkel, John Oertel, Gabriel D. Orebi Gann , et al. (5 additional authors not shown)

Abstract: MiniCLEAN is a single-phase liquid argon dark matter experiment. During the initial cooling phase, impurities within the cold gas ($<$140 K) were monitored by measuring the scintillation light triplet lifetime, and ultimately a triplet lifetime of 3.480 $\pm$ 0.001 (stat.) $\pm$ 0.064 (sys.) $μ$s was obtained, indicating ultra-pure argon. This is the longest argon triplet time constant ever report… ▽ More MiniCLEAN is a single-phase liquid argon dark matter experiment. During the initial cooling phase, impurities within the cold gas ($<$140 K) were monitored by measuring the scintillation light triplet lifetime, and ultimately a triplet lifetime of 3.480 $\pm$ 0.001 (stat.) $\pm$ 0.064 (sys.) $μ$s was obtained, indicating ultra-pure argon. This is the longest argon triplet time constant ever reported. The effect of quenching of separate components of the scintillation light is also investigated. △ Less

Submitted 29 August, 2019; v1 submitted 15 March, 2019; originally announced March 2019.

Design and construction of the MicroBooNE Cosmic Ray Tagger system

Authors: MicroBooNE collaboration, C. Adams, M. Alrashed, R. An, J. Anthony, J. Asaadi, A. Ashkenazi, M. Auger, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, M. Bass, F. Bay, A. Bhat, K. Bhattacharya, M. Bishai, A. Blake, T. Bolton, L. Camilleri, D. Caratelli, I. Caro Terrazas, R. Carr, R. Castillo Fernandez, F. Cavanna , et al. (149 additional authors not shown)

Abstract: The MicroBooNE detector utilizes a liquid argon time projection chamber (LArTPC) with an 85 t active mass to study neutrino interactions along the Booster Neutrino Beam (BNB) at Fermilab. With a deployment location near ground level, the detector records many cosmic muon tracks in each beam-related detector trigger that can be misidentified as signals of interest. To reduce these cosmogenic backgr… ▽ More The MicroBooNE detector utilizes a liquid argon time projection chamber (LArTPC) with an 85 t active mass to study neutrino interactions along the Booster Neutrino Beam (BNB) at Fermilab. With a deployment location near ground level, the detector records many cosmic muon tracks in each beam-related detector trigger that can be misidentified as signals of interest. To reduce these cosmogenic backgrounds, we have designed and constructed a TPC-external Cosmic Ray Tagger (CRT). This sub-system was developed by the Laboratory for High Energy Physics (LHEP), Albert Einstein center for fundamental physics, University of Bern. The system utilizes plastic scintillation modules to provide precise time and position information for TPC-traversing particles. Successful matching of TPC tracks and CRT data will allow us to reduce cosmogenic background and better characterize the light collection system and LArTPC data using cosmic muons. In this paper we describe the design and installation of the MicroBooNE CRT system and provide an overview of a series of tests done to verify the proper operation of the system and its components during installation, commissioning, and physics data-taking. △ Less

Submitted 15 March, 2019; v1 submitted 9 January, 2019; originally announced January 2019.

Rejecting cosmic background for exclusive neutrino interaction studies with Liquid Argon TPCs; a case study with the MicroBooNE detector

Authors: MicroBooNE collaboration, C. Adams, M. Alrashed, R. An, J. Anthony, J. Asaadi, A. Ashkenazi, M. Auger, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, M. Bass, F. Bay, A. Bhat, K. Bhattacharya, M. Bishai, A. Blake, T. Bolton, L. Camilleri, D. Caratelli, I. Caro Terrazas, R. Carr, R. Castillo Fernandez, F. Cavanna , et al. (150 additional authors not shown)

Abstract: Cosmic ray (CR) interactions can be a challenging source of background for neutrino oscillation and cross-section measurements in surface detectors. We present methods for CR rejection in measurements of charged-current quasielastic-like (CCQE-like) neutrino interactions, with a muon and a proton in the final state, measured using liquid argon time projection chambers (LArTPCs). Using a sample of… ▽ More Cosmic ray (CR) interactions can be a challenging source of background for neutrino oscillation and cross-section measurements in surface detectors. We present methods for CR rejection in measurements of charged-current quasielastic-like (CCQE-like) neutrino interactions, with a muon and a proton in the final state, measured using liquid argon time projection chambers (LArTPCs). Using a sample of cosmic data collected with the MicroBooNE detector, mixed with simulated neutrino scattering events, a set of event selection criteria is developed that produces an event sample with minimal contribution from CR background. Depending on the selection criteria used a purity between 50% and 80% can be achieved with a signal selection efficiency between 50% and 25%, with higher purity coming at the expense of lower efficiency. While using a specific dataset from the MicroBooNE detector and selection criteria values optimized for CCQE-like events, the concepts presented here are generic and can be adapted for various studies of exclusive νμ interactions in LArTPCs. △ Less

Submitted 2 January, 2019; v1 submitted 9 December, 2018; originally announced December 2018.

Comments: 12 pages, 10 figures, 1 table

Report number: FERMILAB-PUB-18-677-E

First Measurement of $ν_μ$ Charged-Current $π^{0}$ Production on Argon with a LArTPC

Authors: MicroBooNE collaboration, C. Adams, M. Alrashed, R. An, J. Anthony, J. Asaadi, A. Ashkenazi, M. Auger, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, M. Bass, F. Bay, A. Bhat, K. Bhattacharya, M. Bishai, A. Blake, T. Bolton, L. Camilleri, D. Caratelli, I. Caro Terrazas, R. Carr, R. Castillo Fernandez, F. Cavanna , et al. (150 additional authors not shown)

Abstract: We report the first measurement of the flux-integrated cross section of $ν_μ$ charged-current single $π^{0}$ production on argon. This measurement is performed with the MicroBooNE detector, an 85 ton active mass liquid argon time projection chamber exposed to the Booster Neutrino Beam at Fermilab. This result on argon is compared to past measurements on lighter nuclei to investigate the scaling as… ▽ More We report the first measurement of the flux-integrated cross section of $ν_μ$ charged-current single $π^{0}$ production on argon. This measurement is performed with the MicroBooNE detector, an 85 ton active mass liquid argon time projection chamber exposed to the Booster Neutrino Beam at Fermilab. This result on argon is compared to past measurements on lighter nuclei to investigate the scaling assumptions used in models of the production and transport of pions in neutrino-nucleus scattering. The techniques used are an important demonstration of the successful reconstruction and analysis of neutrino interactions producing electromagnetic final states using a liquid argon time projection chamber operating at the earth's surface. △ Less

Submitted 6 November, 2018; originally announced November 2018.

Journal ref: Phys. Rev. D 99, 091102 (2019)

Different patterns of social closeness observed in mobile phone communication

Authors: Mikaela Irene D. Fudolig, Daniel Monsivais, Kunal Bhattacharya, Hang-Hyun Jo, Kimmo Kaski

Abstract: We analyze a large-scale mobile phone call dataset containing information on the age, gender, and billing locality of users to get insight into social closeness in pairs of individuals of similar age. We show that in addition to using the demographic information, the ranking of contacts by their call frequency in egocentric networks is crucial to characterize the different communication patterns.… ▽ More We analyze a large-scale mobile phone call dataset containing information on the age, gender, and billing locality of users to get insight into social closeness in pairs of individuals of similar age. We show that in addition to using the demographic information, the ranking of contacts by their call frequency in egocentric networks is crucial to characterize the different communication patterns. We find that mutually top-ranked opposite-gender pairs show the highest levels of call frequency and daily regularity, which is consistent with the behavior of real-life romantic partners. At somewhat lower level of call frequency and daily regularity come the mutually top-ranked same-gender pairs, while the lowest call frequency and daily regularity are observed for mutually non-top-ranked pairs. We have also observed that older pairs tend to call less frequently and less regularly than younger pairs, while the average call durations exhibit a more complex dependence on age. We expect that a more detailed analysis can help us better characterize the nature of relationships between pairs of individuals and distinguish between various types of relations, such as siblings, friends, and romantic partners. △ Less

Submitted 7 August, 2019; v1 submitted 30 August, 2018; originally announced August 2018.

Comments: 17 pages, 5 figures, 2 tables. Journal of Computational Social Science (published online, 2019)

A Deep Neural Network for Pixel-Level Electromagnetic Particle Identification in the MicroBooNE Liquid Argon Time Projection Chamber

Authors: MicroBooNE collaboration, C. Adams, M. Alrashed, R. An, J. Anthony, J. Asaadi, A. Ashkenazi, M. Auger, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, M. Bass, F. Bay, A. Bhat, K. Bhattacharya, M. Bishai, A. Blake, T. Bolton, L. Camilleri, D. Caratelli, I. Caro Terrazas, R. Carr, R. Castillo Fernandez, F. Cavanna , et al. (148 additional authors not shown)

Abstract: We have developed a convolutional neural network (CNN) that can make a pixel-level prediction of objects in image data recorded by a liquid argon time projection chamber (LArTPC) for the first time. We describe the network design, training techniques, and software tools developed to train this network. The goal of this work is to develop a complete deep neural network based data reconstruction cha… ▽ More We have developed a convolutional neural network (CNN) that can make a pixel-level prediction of objects in image data recorded by a liquid argon time projection chamber (LArTPC) for the first time. We describe the network design, training techniques, and software tools developed to train this network. The goal of this work is to develop a complete deep neural network based data reconstruction chain for the MicroBooNE detector. We show the first demonstration of a network's validity on real LArTPC data using MicroBooNE collection plane images. The demonstration is performed for stopping muon and a $ν_μ$ charged current neutral pion data samples. △ Less

Submitted 22 August, 2018; originally announced August 2018.

Journal ref: Phys. Rev. D 99, 092001 (2019)

Surface modified mesoporous g-C3N4@FeNi3 as prompt and proficient magnetic adsorbent for crude oil recovery

Authors: Meenakshi Talukdar, Sushant Kumar Behera, Kakoli Bhattacharya, Pritam Deb

Abstract: Efficient oil adsorption and recovery is a generous universal importance for future energy demand and environmental protection. Adsorbents based on 2D flatland with engineered surfaces can overcome the limitations of conventional methods for selective oil adsorption. Here, we report magnetic hydrophobic/oleophilic graphitic C3N4 nanosheets that exhibit excellent oil sorption performance and rapid… ▽ More Efficient oil adsorption and recovery is a generous universal importance for future energy demand and environmental protection. Adsorbents based on 2D flatland with engineered surfaces can overcome the limitations of conventional methods for selective oil adsorption. Here, we report magnetic hydrophobic/oleophilic graphitic C3N4 nanosheets that exhibit excellent oil sorption performance and rapid removal of adsorbed oil using an external magnet. Combining porous and nanosheets structure along with magnetic FeNi3 and fatty acid surface functionalization make the system an efficient adsorbent for adsorbing and separating crude oil from water. The graphitic sheets selectively adsorb crude oil with enhancement of thickness up to 9 folds and mass by 4.5 times than the pristine nanocomposite system. The smart adsorption property of g-C3N4@FeNi3 has been realized through comprehensive adsorption kinetics and inclusive isotherm studies. The nanocomposite can be further recycled and reused in an eco-friendly manner for oil adsorption and recovery. △ Less

Submitted 19 December, 2018; v1 submitted 8 June, 2018; originally announced June 2018.

Comments: 4 pages, 5 figures

Comparison of νμ-Ar multiplicity distributions observed by MicroBooNE to GENIE model predictions

Authors: C. Adams, R. An, J. Anthony, J. Asaadi, M. Auger, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, M. Bass, F. Bay, A. Bhat, K. Bhattacharya, M. Bishai, A. Blake, T. Bolton, L. Camilleri, D. Caratelli, R. Castillo Fernandez, F. Cavanna, G. Cerati, H. Chen, Y. Chen, E. Church, D. Cianci , et al. (140 additional authors not shown)

Abstract: We measure a large set of observables in inclusive charged current muon neutrino scattering on argon with the MicroBooNE liquid argon time projection chamber operating at Fermilab. We evaluate three neutrino interaction models based on the widely used GENIE event generator using these observables. The measurement uses a data set consisting of neutrino interactions with a final state muon candidate… ▽ More We measure a large set of observables in inclusive charged current muon neutrino scattering on argon with the MicroBooNE liquid argon time projection chamber operating at Fermilab. We evaluate three neutrino interaction models based on the widely used GENIE event generator using these observables. The measurement uses a data set consisting of neutrino interactions with a final state muon candidate fully contained within the MicroBooNE detector. These data were collected in 2016 with the Fermilab Booster Neutrino Beam, which has an average neutrino energy of 800 MeV, using an exposure corresponding to 5E19 protons-on-target. The analysis employs fully automatic event selection and charged particle track reconstruction and uses a data-driven technique to separate neutrino interactions from cosmic ray background events. We find that GENIE models consistently describe the shapes of a large number of kinematic distributions for fixed observed multiplicity. △ Less

Submitted 26 March, 2019; v1 submitted 17 May, 2018; originally announced May 2018.

Comments: 31 pages, 39 figures, 10 tables

Journal ref: Eur. Phys. J. C (2019) 79:248

Social Physics: Uncovering Human Behaviour from Communication

Authors: Kunal Bhattacharya, Kimmo Kaski

Abstract: In the post year 2000 era the technologies that facilitate human communication have rapidly multiplied. While the adoption of these technologies has hugely impacted the behaviour and sociality of people, specifically in urban but also in rural environments, their "digital footprints" on different data bases have become an active area of research. The existence and accessibility of such large popul… ▽ More In the post year 2000 era the technologies that facilitate human communication have rapidly multiplied. While the adoption of these technologies has hugely impacted the behaviour and sociality of people, specifically in urban but also in rural environments, their "digital footprints" on different data bases have become an active area of research. The existence and accessibility of such large population-level datasets, has allowed scientists to study and model innate human tendencies and social patterns in an unprecedented way that complements traditional research approaches like questionnaire studies. In this review we focus on data analytics and modelling research - we call Social Physics - as it has been carried out using the mobile phone data sets to get insight into the various aspects of human sociality, burstiness in communication, mobility patterns, and daily rhythms. △ Less

Submitted 13 April, 2018; originally announced April 2018.

Comments: Review, 38 pages, 8 figures

Ionization Electron Signal Processing in Single Phase LArTPCs II. Data/Simulation Comparison and Performance in MicroBooNE

Authors: MicroBooNE collaboration, C. Adams, R. An, J. Anthony, J. Asaadi, M. Auger, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, M. Bass, F. Bay, A. Bhat, K. Bhattacharya, M. Bishai, A. Blake, T. Bolton, L. Camilleri, D. Caratelli, R. Carr, I. Caro Terrazas, R. Castillo Fernandez, F. Cavanna, G. Cerati, H. Chen , et al. (146 additional authors not shown)

Abstract: The single-phase liquid argon time projection chamber (LArTPC) provides a large amount of detailed information in the form of fine-grained drifted ionization charge from particle traces. To fully utilize this information, the deposited charge must be accurately extracted from the raw digitized waveforms via a robust signal processing chain. Enabled by the ultra-low noise levels associated with cry… ▽ More The single-phase liquid argon time projection chamber (LArTPC) provides a large amount of detailed information in the form of fine-grained drifted ionization charge from particle traces. To fully utilize this information, the deposited charge must be accurately extracted from the raw digitized waveforms via a robust signal processing chain. Enabled by the ultra-low noise levels associated with cryogenic electronics in the MicroBooNE detector, the precise extraction of ionization charge from the induction wire planes in a single-phase LArTPC is qualitatively demonstrated on MicroBooNE data with event display images, and quantitatively demonstrated via waveform-level and track-level metrics. Improved performance of induction plane calorimetry is demonstrated through the agreement of extracted ionization charge measurements across different wire planes for various event topologies. In addition to the comprehensive waveform-level comparison of data and simulation, a calibration of the cryogenic electronics response is presented and solutions to various MicroBooNE-specific TPC issues are discussed. This work presents an important improvement in LArTPC signal processing, the foundation of reconstruction and therefore physics analyses in MicroBooNE. △ Less

Submitted 11 June, 2018; v1 submitted 7 April, 2018; originally announced April 2018.

Comments: 54 pages, 36 figures; the first part of this work can be found at arXiv:1802.08709

Journal ref: JINST 13 P07007 (2018)

Group formation on a small-world: experiment and modelling

Authors: Kunal Bhattacharya, Tuomas Takko, Daniel Monsivais, Kimmo Kaski

Abstract: As a step towards studying human-agent collectives we conduct an online game with human participants cooperating on a network. The game is presented in the context of achieving group formation through local coordination. The players set initially to a small world network with limited information on the location of other players, coordinate their movements to arrange themselves into groups. To unde… ▽ More As a step towards studying human-agent collectives we conduct an online game with human participants cooperating on a network. The game is presented in the context of achieving group formation through local coordination. The players set initially to a small world network with limited information on the location of other players, coordinate their movements to arrange themselves into groups. To understand the decision making process we construct a data-driven model of agents based on probability matching. The model allows us to gather insight into the nature and degree of rationality employed by the human players. By varying the parameters in agent based simulations we are able to benchmark the human behaviour. We observe that while the players utilize the neighbourhood information in limited capacity, the perception of risk is optimal. We also find that for certain parameter ranges the agents are able to act more efficiently when compared to the human players. This approach would allow us to simulate the collective dynamics in games with agents having varying strategies playing alongside human proxies. △ Less

Submitted 10 July, 2019; v1 submitted 2 March, 2018; originally announced March 2018.

Comments: 23 pages, 9 figures

Journal ref: J. R. Soc. Interface 16:20180814 (2019)

Ionization Electron Signal Processing in Single Phase LArTPCs I. Algorithm Description and Quantitative Evaluation with MicroBooNE Simulation

Authors: MicroBooNE collaboration, C. Adams, R. An, J. Anthony, J. Asaadi, M. Auger, L. Bagby, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, M. Bass, F. Bay, A. Bhat, K. Bhattacharya, M. Bishai, A. Blake, T. Bolton, L. Camilleri, D. Caratelli, R. Castillo Fernandez, F. Cavanna, G. Cerati, H. Chen, Y. Chen , et al. (144 additional authors not shown)

Abstract: We describe the concept and procedure of drifted-charge extraction developed in the MicroBooNE experiment, a single-phase liquid argon time projection chamber (LArTPC). This technique converts the raw digitized TPC waveform to the number of ionization electrons passing through a wire plane at a given time. A robust recovery of the number of ionization electrons from both induction and collection a… ▽ More We describe the concept and procedure of drifted-charge extraction developed in the MicroBooNE experiment, a single-phase liquid argon time projection chamber (LArTPC). This technique converts the raw digitized TPC waveform to the number of ionization electrons passing through a wire plane at a given time. A robust recovery of the number of ionization electrons from both induction and collection anode wire planes will augment the 3D reconstruction, and is particularly important for tomographic reconstruction algorithms. A number of building blocks of the overall procedure are described. The performance of the signal processing is quantitatively evaluated by comparing extracted charge with the true charge through a detailed TPC detector simulation taking into account position-dependent induced current inside a single wire region and across multiple wires. Some areas for further improvement of the performance of the charge extraction procedure are also discussed. △ Less

Submitted 9 April, 2018; v1 submitted 23 February, 2018; originally announced February 2018.

Comments: 60 pages, 36 figures. The second part of this work can be found at arXiv:1804.02583

Journal ref: JINST 13 P07006 (2018)

Peer relations with mobile phone data: Best friends and family formation

Authors: Tamas David-Barrett, Anna Rotkirch, Asim Ghosh, Kunal Bhattacharya, Daniel Monsivais, Isabel Behncke, Janos Kertesz, Kimmo Kaski

Abstract: Earlier attempts to investigate the changes of the role of friendship in different life stages have failed due to lack of data. We close this gap by using a large data set of mobile phone calls from a European country in 2007, to study how the people's call patterns to their close social contacts are associated with age and gender of the callers. We hypothesize that (i) communication with peers, d… ▽ More Earlier attempts to investigate the changes of the role of friendship in different life stages have failed due to lack of data. We close this gap by using a large data set of mobile phone calls from a European country in 2007, to study how the people's call patterns to their close social contacts are associated with age and gender of the callers. We hypothesize that (i) communication with peers, defined as callers of similar age, will be most important during the period of family formation and that (ii) the importance of best friends defined as same-sex callers of exactly the same age, will be stronger for women than for men. Results show that the frequency of phone calls with the same-sex peers in this population turns out to be relatively stable through life for both men and women. In line with the first hypothesis, there was a significant increase in the length of the phone calls for callers between ages 30 to 40 years. Partly in line with the second hypothesis, the increase in phone calls turned out to be particularly pronounced among females, although there were only minor gender differences in call frequencies. Furthermore, women tended to have long phone conversations with their same-age female friend, and also with somewhat older peers. In sum, we provide evidence from big data for the adult life stages at which peers are most important, and suggest that best friends appear to have a niche of their own in human sociality. △ Less

Submitted 25 August, 2017; originally announced August 2017.

Comments: 13 pages, 2 figures

Network of families in a contemporary population: regional and cultural assortativity

Authors: Kunal Bhattacharya, Venla Berg, Asim Ghosh, Daniel Monsivais, Janos Kertesz, Anna Rotkirch, Kimmo Kaski

Abstract: Using a large dataset with individual-level demographic information of 60,000 families in contemporary Finland, we analyse the variation and cultural assortativity in a network of families. Families are considered as vertices and unions between males and females who have a common child and belong to different families are considered as edges in such a network of families. The sampled network is a… ▽ More Using a large dataset with individual-level demographic information of 60,000 families in contemporary Finland, we analyse the variation and cultural assortativity in a network of families. Families are considered as vertices and unions between males and females who have a common child and belong to different families are considered as edges in such a network of families. The sampled network is a collection of many disjoint components with the largest connected component being dominated by families rooted in one specific region. We characterize the network in terms of the basic structural properties and then explore the network transitivity and assortativity with regards to regions of origin and linguistic identity. Transitivity is seen to result from linguistic homophily in the network. Overall, our results demonstrate that geographic proximity and language strongly influence the structuring of network. △ Less

Submitted 22 August, 2017; v1 submitted 21 August, 2017; originally announced August 2017.

Comments: 25 pages, 12 figures

Migration patterns across the life course of families: Gender differences and proximity with parents and siblings in Finland

Authors: Asim Ghosh, Venla Berg, Kunal Bhattacharya, Daniel Monsivais, Janos Kertesz, Kimmo Kaski, Anna Rotkirch

Abstract: Family members' life course tendencies to remain geographically close to each other or to migrate due to education or job opportunities have been studied relatively little. Here we investigate migration patterns of parents and their children between 19 administrative regions of Finland from 1970 to 2012. Using the FinnFamily register dataset of 60 000 index individuals and their family members, we… ▽ More Family members' life course tendencies to remain geographically close to each other or to migrate due to education or job opportunities have been studied relatively little. Here we investigate migration patterns of parents and their children between 19 administrative regions of Finland from 1970 to 2012. Using the FinnFamily register dataset of 60 000 index individuals and their family members, we investigate the patterns of regional migration and regional co-residence of parents and their children. Specifically, we analyse how likely it is for children to reside in the same region as their parents at any specific age, whether parents and children who live in different regions are likely to reunite, and whether siblings function as regional attractors to each other. Results show an intense regional migration of people to the capital area. The migration propensity of individuals is high in early childhood and peaks in early adulthood. About two thirds of Finnish children live in the same region as their parents throughout their adult lives. Females show higher propensity to migrate than males, since daughters move away from their parents earlier and with a higher rate than sons do. The propensity for two full sibling brothers to be in the same region is higher than that for other types of sibling dyads. We conclude that family members serve as important geographical attractors to each other through the life course and that family attraction is stronger for sons and brothers than for daughters and sisters in contemporary Finland. △ Less

Submitted 8 August, 2017; originally announced August 2017.

Comments: 17 pages (14 from the main text and 3 from the SI), 10 figures in the main text, 2 in the SI. 1 table in the main text, 3 in the SI

Tracking Urban Human Activity from Mobile Phone Calling Patterns

Authors: Daniel Monsivais, Asim Ghosh, Kunal Bhattacharya, Robin I. M Dunbar, Kimmo Kaski

Abstract: Timings of human activities are marked by circadian clocks which in turn are entrained to different environmental signals. In an urban environment the presence of artificial lighting and various social cues tend to disrupt the natural entrainment with the sunlight. However, it is not completely understood to what extent this is the case. Here we exploit the large-scale data analysis techniques to… ▽ More Timings of human activities are marked by circadian clocks which in turn are entrained to different environmental signals. In an urban environment the presence of artificial lighting and various social cues tend to disrupt the natural entrainment with the sunlight. However, it is not completely understood to what extent this is the case. Here we exploit the large-scale data analysis techniques to study the mobile phone calling activity of people in large cities to infer the dynamics of urban daily rhythms. From the calling patterns of about 1,000,000 users spread over different cities but lying inside the same time-zone, we show that the onset and termination of the calling activity synchronizes with the east-west progression of the sun. We also find that the onset and termination of the calling activity of users follows a yearly dynamics, varying across seasons, and that its timings are entrained to solar midnight. Furthermore, we show that the average mid-sleep time of people living in urban areas depends on the age and gender of each cohort as a result of biological and social factors. △ Less

Submitted 17 October, 2017; v1 submitted 20 April, 2017; originally announced April 2017.

Comments: This manuscript includes and extends content from our old manuscript "Seasonal and geographical influence on sleeping patterns inferred from mobile phone data", which was split into two parts due to its length and new investigations added

Quantifying gender preferences across humans lifespan

Authors: Asim Ghosh, Daniel Monsivais, Kunal Bhattacharya, Robin I. M. Dunbar, Kimmo Kaski

Abstract: In human relations individuals' gender and age play a key role in the structures and dynamics of their social arrangements. In order to analyze the gender preferences of individuals in interaction with others at different stages of their lives we study a large mobile phone dataset. To do this we consider four fundamental gender-related caller and callee combinations of human interactions, namely m… ▽ More In human relations individuals' gender and age play a key role in the structures and dynamics of their social arrangements. In order to analyze the gender preferences of individuals in interaction with others at different stages of their lives we study a large mobile phone dataset. To do this we consider four fundamental gender-related caller and callee combinations of human interactions, namely male to male, male to female, female to male, and female to female, which together with age, kinship, and different levels of friendship give rise to a wide scope of human sociality. Here we analyse the relative strength of these four types of interaction using a large dataset of mobile phone communication records. Our analysis suggests strong age dependence for an ego of one gender choosing to call an individual of either gender. We observe a strong opposite sex bonding across most of their reproductive age. However, older women show a strong tendency to connect to another female that is one generation younger in a way that is suggestive of the \emph{grandmothering effect}. We also find that the relative strength among the four possible interactions depends on phone call duration. For calls of medium and long duration, opposite gender interactions are significantly more probable than same gender interactions during the reproductive years, suggesting potential emotional exchange between spouses. By measuring the fraction of calls to other generations we find that mothers tend to make calls more to their daughters than to their sons, whereas fathers make calls more to their sons than to their daughters. For younger people, most of their calls go to same generation alters, while older people call the younger people more frequently, which supports the suggestion that \emph{affection flows downward}. △ Less

Submitted 18 November, 2016; originally announced November 2016.

Comments: 12 pages, 8 figures

Spectrum-splitting approach for Fermi-operator expansion in all-electron Kohn-Sham DFT calculations

Authors: Phani Motamarri, Vikram Gavini, Kaushik Bhattacharya, Michael Ortiz

Abstract: We present a spectrum-splitting approach to conduct all-electron Kohn-Sham density functional theory (DFT) calculations by employing Fermi-operator expansion of the Kohn-Sham Hamiltonian. The proposed approach splits the subspace containing the occupied eigenspace into a core-subspace, spanned by the core eigenfunctions, and its complement, the valence-subspace, and thereby enables an efficient co… ▽ More We present a spectrum-splitting approach to conduct all-electron Kohn-Sham density functional theory (DFT) calculations by employing Fermi-operator expansion of the Kohn-Sham Hamiltonian. The proposed approach splits the subspace containing the occupied eigenspace into a core-subspace, spanned by the core eigenfunctions, and its complement, the valence-subspace, and thereby enables an efficient computation of the Fermi-operator expansion by reducing the expansion to the valence-subspace projected Kohn-Sham Hamiltonian. The key ideas used in our approach are: (i) employ Chebyshev filtering to compute a subspace containing the occupied states followed by a localization procedure to generate non-orthogonal localized functions spanning the Chebyshev-filtered subspace; (ii) compute the Kohn-Sham Hamiltonian projected onto the valence-subspace; (iii) employ Fermi-operator expansion in terms of the valence-subspace projected Hamiltonian to compute the density matrix, electron-density and band energy. We demonstrate the accuracy and performance of the method on benchmark materials systems involving silicon nano-clusters up to 1330 electrons, a single gold atom and a six-atom gold nano-cluster. The benchmark studies on silicon nano-clusters revealed a staggering five-fold reduction in the Fermi-operator expansion polynomial degree by using the spectrum-splitting approach for accuracies in the ground-state energies of $\sim 10^{-4} Ha/atom$ with respect to reference calculations. Further, numerical investigations on gold suggest that spectrum splitting is indispensable to achieve meaningful accuracies, while employing Fermi-operator expansion. △ Less

Submitted 28 August, 2016; originally announced August 2016.

Comments: 12 pages, 5 figures