Showing 1–50 of 97 results for author: Park, J W

CAT: Contrastive Adapter Training for Personalized Image Generation

Authors: Jae Wan Park, Sang Hyun Park, Jun Young Koh, Junha Lee, Min Song

Abstract: The emergence of various adapters, including Low-Rank Adaptation (LoRA) applied from the field of natural language processing, has allowed diffusion models to personalize image generation at a low cost. However, due to the various challenges including limited datasets and shortage of regularization and computation resources, adapter training often results in unsatisfactory outcomes, leading to the… ▽ More The emergence of various adapters, including Low-Rank Adaptation (LoRA) applied from the field of natural language processing, has allowed diffusion models to personalize image generation at a low cost. However, due to the various challenges including limited datasets and shortage of regularization and computation resources, adapter training often results in unsatisfactory outcomes, leading to the corruption of the backbone model's prior knowledge. One of the well known phenomena is the loss of diversity in object generation, especially within the same class which leads to generating almost identical objects with minor variations. This poses challenges in generation capabilities. To solve this issue, we present Contrastive Adapter Training (CAT), a simple yet effective strategy to enhance adapter training through the application of CAT loss. Our approach facilitates the preservation of the base model's original knowledge when the model initiates adapters. Furthermore, we introduce the Knowledge Preservation Score (KPS) to evaluate CAT's ability to keep the former information. We qualitatively and quantitatively compare CAT's improvement. Finally, we mention the possibility of CAT in the aspects of multi-concept adapter and optimization. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: CVPRW 2024

Dual-species Bose-Einstein condensates of $^{23}$Na and $^{41}$K with tunable interactions

Authors: Jaeryeong Chang, Sungjun Lee, Yoonsoo Kim, Younghoon Lim, Jee Woo Park

Abstract: We report the creation of dual-species Bose-Einstein condensates (BECs) of $^{23}$Na and $^{41}$K. Favorable background scattering lengths enable efficient sympathetic cooling of $^{41}$K via forced evaporative cooling of $^{23}$Na in a plugged magnetic trap and an optical dipole trap. The $1/e$ lifetime of the thermal mixture in the stretched hyperfine state exceeds 5 s in the presence of backgro… ▽ More We report the creation of dual-species Bose-Einstein condensates (BECs) of $^{23}$Na and $^{41}$K. Favorable background scattering lengths enable efficient sympathetic cooling of $^{41}$K via forced evaporative cooling of $^{23}$Na in a plugged magnetic trap and an optical dipole trap. The $1/e$ lifetime of the thermal mixture in the stretched hyperfine state exceeds 5 s in the presence of background scattering. At the end of evaporation, we create dual BECs in the immiscible phase, with about $3\times10^5$ $^{23}$Na atoms surrounding $5\times10^4$ $^{41}$K atoms. To further enable the tuning of the interspecies interaction strength, we locate multiple Feshbach resonances at magnetic fields up to 100 G. The broadest $s$-wave resonance located at 73.4(3) G features a favorable width of 1.8(2) G. This work sets the stage for the creation of ultracold gases of strongly dipolar bosonic $^{23}$Na$^{41}$K molecules as well as the exploration of many-body physics in bosonic $^{23}$Na-$^{41}$K mixtures. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 6 pages, 4 figures

Integration of Graph Neural Network and Neural-ODEs for Tumor Dynamic Prediction

Authors: Omid Bazgir, Zichen Wang, Ji Won Park, Marc Hafner, James Lu

Abstract: In anti-cancer drug development, a major scientific challenge is disentangling the complex relationships between high-dimensional genomics data from patient tumor samples, the corresponding tumor's organ of origin, the drug targets associated with given treatments and the resulting treatment response. Furthermore, to realize the aspirations of precision medicine in identifying and adjusting treatm… ▽ More In anti-cancer drug development, a major scientific challenge is disentangling the complex relationships between high-dimensional genomics data from patient tumor samples, the corresponding tumor's organ of origin, the drug targets associated with given treatments and the resulting treatment response. Furthermore, to realize the aspirations of precision medicine in identifying and adjusting treatments for patients depending on the therapeutic response, there is a need for building tumor dynamic models that can integrate both longitudinal tumor size as well as multimodal, high-content data. In this work, we take a step towards enhancing personalized tumor dynamic predictions by proposing a heterogeneous graph encoder that utilizes a bipartite Graph Convolutional Neural network (GCN) combined with Neural Ordinary Differential Equations (Neural-ODEs). We applied the methodology to a large collection of patient-derived xenograft (PDX) data, spanning a wide variety of treatments (as well as their combinations) on tumors that originated from a number of different organs. We first show that the methodology is able to discover a tumor dynamic model that significantly improves upon an empirical model which is in current use. Additionally, we show that the graph encoder is able to effectively utilize multimodal data to enhance tumor predictions. Our findings indicate that the methodology holds significant promise and offers potential applications in pre-clinical settings. △ Less

Submitted 27 March, 2024; v1 submitted 2 October, 2023; originally announced October 2023.

Emergent Quantum Phenomena of Noncentrosymmetric Charge-Density Wave in 1T-Transition Metal Dichalcogenides

Authors: Cheong-Eung Ahn, Kyung-Hwan Jin, Young-Jae Choi, Jae Whan Park, Han Woong Yeom, Ara Go, Yong Baek Kim, Gil Young Cho

Abstract: 1T-transition metal dichalcogenides (TMD) have been an exciting platform for exploring the intertwinement of charge density waves and strong correlation phenomena. While the David star structure has been conventionally considered as the underlying charge order in the literature, recent scanning tunneling probe experiments on several monolayer 1T-TMD materials have motivated a new, alternative stru… ▽ More 1T-transition metal dichalcogenides (TMD) have been an exciting platform for exploring the intertwinement of charge density waves and strong correlation phenomena. While the David star structure has been conventionally considered as the underlying charge order in the literature, recent scanning tunneling probe experiments on several monolayer 1T-TMD materials have motivated a new, alternative structure, namely the anion-centered David star structure. In this Letter, we show that this novel anion-centered David star structure manifestly breaks inversion symmetry, resulting in flat bands with pronounced Rashba spin-orbit couplings. These distinctive features unlock novel possibilities and functionalities for 1T-TMDs, including the giant spin Hall effect, the emergence of Chern bands, and spin liquid that spontaneously breaks crystalline rotational symmetry. Our findings establish promising avenues for exploring emerging quantum phenomena of monolayer 1T-TMDs with this novel noncentrosymmetric structure. △ Less

Submitted 14 June, 2024; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 4.5+12 pages, 4+14 figures

Journal ref: Phys. Rev. Lett. 132, 226401 (2024)

Blind Biological Sequence Denoising with Self-Supervised Set Learning

Authors: Nathan Ng, Ji Won Park, Jae Hyeon Lee, Ryan Lewis Kelly, Stephen Ra, Kyunghyun Cho

Abstract: Biological sequence analysis relies on the ability to denoise the imprecise output of sequencing platforms. We consider a common setting where a short sequence is read out repeatedly using a high-throughput long-read platform to generate multiple subreads, or noisy observations of the same sequence. Denoising these subreads with alignment-based approaches often fails when too few subreads are avai… ▽ More Biological sequence analysis relies on the ability to denoise the imprecise output of sequencing platforms. We consider a common setting where a short sequence is read out repeatedly using a high-throughput long-read platform to generate multiple subreads, or noisy observations of the same sequence. Denoising these subreads with alignment-based approaches often fails when too few subreads are available or error rates are too high. In this paper, we propose a novel method for blindly denoising sets of sequences without directly observing clean source sequence labels. Our method, Self-Supervised Set Learning (SSSL), gathers subreads together in an embedding space and estimates a single set embedding as the midpoint of the subreads in both the latent and sequence spaces. This set embedding represents the "average" of the subreads and can be decoded into a prediction of the clean sequence. In experiments on simulated long-read DNA data, SSSL methods denoise small reads of $\leq 6$ subreads with 17% fewer errors and large reads of $>6$ subreads with 8% fewer errors compared to the best baseline. On a real dataset of antibody sequences, SSSL improves over baselines on two self-supervised metrics, with a significant improvement on difficult small reads that comprise over 60% of the test set. By accurately denoising these reads, SSSL promises to better realize the potential of high-throughput DNA sequencing data for downstream scientific applications. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Kinkless electronic junction along one dimensional electronic channel

Authors: Qirong Yao, Jae Whan Park, Choongjae Won, Sang-Wook Cheong, Han Woong Yeom

Abstract: Here we report the formation of type-A and type-B electronic junctions without any structural discontinuity along a well-defined 1-nm-wide one-dimensional electronic channel within a van der Waals layer. We employ scanning tunneling microscopy and spectroscopy techniques to investigate the atomic and electronic structure along peculiar domain walls formed on the charge-density-wave phase of 1T-TaS… ▽ More Here we report the formation of type-A and type-B electronic junctions without any structural discontinuity along a well-defined 1-nm-wide one-dimensional electronic channel within a van der Waals layer. We employ scanning tunneling microscopy and spectroscopy techniques to investigate the atomic and electronic structure along peculiar domain walls formed on the charge-density-wave phase of 1T-TaS2. We find distinct kinds of abrupt electronic junctions with discontinuities of the band gap along the domain walls, which do not have any structural kinks and defects. Our density-functional calculations reveal a novel mechanism of the electronic junction formation; they are formed by a kinked domain wall in the layer underneath through substantial electronic interlayer coupling. This work demonstrates that the interlayer electronic coupling can be an effective control knob over several-nanometer-scale electronic property of two-dimensional atomic monolayers. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Topological soliton molecule in quasi 1D charge density wave

Authors: Taehwan Im, Sun Kyu Song, Jae Whan Park, Han Woong Yeom

Abstract: Soliton molecules, bound states of two solitons, can be important for the informatics using solitons and the quest for exotic particles in a wide range of physical systems from unconventional superconductors to nuclear matter and Higgs field, but have been observed only in temporal dimension for classical wave optical systems. Here, we identify a topological soliton molecule formed spatially in an… ▽ More Soliton molecules, bound states of two solitons, can be important for the informatics using solitons and the quest for exotic particles in a wide range of physical systems from unconventional superconductors to nuclear matter and Higgs field, but have been observed only in temporal dimension for classical wave optical systems. Here, we identify a topological soliton molecule formed spatially in an electronic system, a quasi 1D charge density wave of indium atomic wires. This system is composed of two coupled Peierls chains, which are endowed with a Z$_4$ topology and three distinct, right-chiral, left-chiral, and non-chiral, solitons. Our scanning tunneling microscopy measurements identify a bound state of right- and left-chiral solitons with distinct in-gap states and net zero phase shift. Our density functional theory calculations reveal the attractive interaction of these solitons and the hybridization of their electronic states. This result initiates the study of the interaction between solitons in electronic systems, which can provide novel manybody electronic states and extra data-handling capacity beyond the given soliton topology. △ Less

Submitted 10 August, 2023; originally announced August 2023.

BOtied: Multi-objective Bayesian optimization with tied multivariate ranks

Authors: Ji Won Park, Nataša Tagasovska, Michael Maser, Stephen Ra, Kyunghyun Cho

Abstract: Many scientific and industrial applications require the joint optimization of multiple, potentially competing objectives. Multi-objective Bayesian optimization (MOBO) is a sample-efficient framework for identifying Pareto-optimal solutions. At the heart of MOBO is the acquisition function, which determines the next candidate to evaluate by navigating the best compromises among the objectives. In t… ▽ More Many scientific and industrial applications require the joint optimization of multiple, potentially competing objectives. Multi-objective Bayesian optimization (MOBO) is a sample-efficient framework for identifying Pareto-optimal solutions. At the heart of MOBO is the acquisition function, which determines the next candidate to evaluate by navigating the best compromises among the objectives. In this paper, we show a natural connection between non-dominated solutions and the extreme quantile of the joint cumulative distribution function (CDF). Motivated by this link, we propose the Pareto-compliant CDF indicator and the associated acquisition function, BOtied. BOtied inherits desirable invariance properties of the CDF, and an efficient implementation with copulas allows it to scale to many objectives. Our experiments on a variety of synthetic and real-world problems demonstrate that BOtied outperforms state-of-the-art MOBO acquisition functions while being computationally efficient for many objectives. △ Less

Submitted 7 June, 2024; v1 submitted 1 June, 2023; originally announced June 2023.

Comments: 12 pages (+9 appendix), 13 figures. Accepted at ICML 2024

Chain of Log-Concave Markov Chains

Authors: Saeed Saremi, Ji Won Park, Francis Bach

Abstract: We introduce a theoretical framework for sampling from unnormalized densities based on a smoothing scheme that uses an isotropic Gaussian kernel with a single fixed noise scale. We prove one can decompose sampling from a density (minimal assumptions made on the density) into a sequence of sampling from log-concave conditional densities via accumulation of noisy measurements with equal noise levels… ▽ More We introduce a theoretical framework for sampling from unnormalized densities based on a smoothing scheme that uses an isotropic Gaussian kernel with a single fixed noise scale. We prove one can decompose sampling from a density (minimal assumptions made on the density) into a sequence of sampling from log-concave conditional densities via accumulation of noisy measurements with equal noise levels. Our construction is unique in that it keeps track of a history of samples, making it non-Markovian as a whole, but it is lightweight algorithmically as the history only shows up in the form of a running empirical mean of samples. Our sampling algorithm generalizes walk-jump sampling (Saremi & Hyvärinen, 2019). The "walk" phase becomes a (non-Markovian) chain of (log-concave) Markov chains. The "jump" from the accumulated measurements is obtained by empirical Bayes. We study our sampling algorithm quantitatively using the 2-Wasserstein metric and compare it with various Langevin MCMC algorithms. We also report a remarkable capacity of our algorithm to "tunnel" between modes of a distribution. △ Less

Submitted 28 September, 2023; v1 submitted 30 May, 2023; originally announced May 2023.

Metal halide thermoelectrics: prediction of high-performance CsCu2I3

Authors: Jong Woong Park, Young-Kwang Jung, Aron Walsh

Abstract: Thermoelectric devices can directly convert waste heat into electricity, which makes them an important clean energy technology. The underlying materials performance can be evaluated by the dimensionless figure of merit ZT. Metal halides are attractive candidates due to their chemical flexibility and ease of processing; however, the maximum ZT realized (ZT = 0.15) falls far below the level needed f… ▽ More Thermoelectric devices can directly convert waste heat into electricity, which makes them an important clean energy technology. The underlying materials performance can be evaluated by the dimensionless figure of merit ZT. Metal halides are attractive candidates due to their chemical flexibility and ease of processing; however, the maximum ZT realized (ZT = 0.15) falls far below the level needed for commercialization (ZT > 1). Using a first-principles procedure we assess the thermoelectric potential of copper halide CsCu2I3, which features 1D Cu-I connectivity. The n-type crystal is predicted to exhibit a maximum ZT of 2.2 at 600 K along the b-axis. The strong phonon anharmonicity of this system is shown by locally stable non-centrosymmetric Amm2 structures that are averaged to form the observed centrosymmetric Cmcm space group. Our work provides insights into the structure-property relations in metal halide thermoelectrics and suggests a path forward to engineer higher-performance heat-to-electricity conversion. △ Less

Submitted 19 May, 2023; originally announced May 2023.

Comments: 12 pages, 10 figures

Journal ref: PRX Energy 2, 043004 (2023)

Latent Stochastic Differential Equations for Modeling Quasar Variability and Inferring Black Hole Properties

Authors: Joshua Fagin, Ji Won Park, Henry Best, James Hung-Hsu Chan, K. E Saavik Ford, Matthew J. Graham, V. Ashley Villar, Shirley Ho, Matthew O'Dowd

Abstract: Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar's brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described… ▽ More Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar's brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described by a stochastic differential equation (SDE). Upcoming wide-field telescopes such as the Rubin Observatory Legacy Survey of Space and Time (LSST) are expected to observe tens of millions of AGN in multiple filters over a ten year period, so there is a need for efficient and automated modeling techniques that can handle the large volume of data. Latent SDEs are machine learning models well suited for modeling quasar variability, as they can explicitly capture the underlying stochastic dynamics. In this work, we adapt latent SDEs to jointly reconstruct multivariate quasar light curves and infer their physical properties such as the black hole mass, inclination angle, and temperature slope. Our model is trained on realistic simulations of LSST ten year quasar light curves, and we demonstrate its ability to reconstruct quasar light curves even in the presence of long seasonal gaps and irregular sampling across different bands, outperforming a multioutput Gaussian process regression baseline. Our method has the potential to provide a deeper understanding of the physical properties of quasars and is applicable to a wide range of other multivariate time series with missing data and irregular sampling. △ Less

Submitted 16 April, 2024; v1 submitted 9 April, 2023; originally announced April 2023.

Comments: 22 pages, 12 figures

Journal ref: The Astrophysical Journal, Volume 965, Number 2, April 2024

Mobile Kink Solitons in a Van der Waals Charge-Density-Wave Layer

Authors: Jinwon Lee, Jae Whan Park, Gil-Young Cho, Han Woong Yeom

Abstract: Kinks, point-like geometrical defects along dislocations, domain walls, and DNA, are stable and mobile, as solutions of a sine-Gordon wave equation. While they are widely investigated for crystal deformations and domain wall motions, electronic properties of individual kinks have received little attention. In this work, electronically and topologically distinct kinks are discovered along electroni… ▽ More Kinks, point-like geometrical defects along dislocations, domain walls, and DNA, are stable and mobile, as solutions of a sine-Gordon wave equation. While they are widely investigated for crystal deformations and domain wall motions, electronic properties of individual kinks have received little attention. In this work, electronically and topologically distinct kinks are discovered along electronic domain walls in a correlated van der Waals insulator of 1$T$-TaS$_2$. Mobile kinks and antikinks are identified as trapped by pinning defects and imaged in scanning tunneling microscopy. Their atomic structures and in-gap electronic states are unveiled, which are mapped approximately into Su-Schrieffer-Heeger solitons. The twelve-fold degeneracy of the domain walls in the present system guarantees an extraordinarily large number of distinct kinks and antikinks to emerge. Such large degeneracy together with the robust geometrical nature may be useful for handling multilevel information in van der Waals materials architectures. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Comments: 12 pages, 4 figures

Journal ref: Advanced Materials (2023)

SupSiam: Non-contrastive Auxiliary Loss for Learning from Molecular Conformers

Authors: Michael Maser, Ji Won Park, Joshua Yao-Yu Lin, Jae Hyeon Lee, Nathan C. Frey, Andrew Watkins

Abstract: We investigate Siamese networks for learning related embeddings for augmented samples of molecular conformers. We find that a non-contrastive (positive-pair only) auxiliary task aids in supervised training of Euclidean neural networks (E3NNs) and increases manifold smoothness (MS) around point-cloud geometries. We demonstrate this property for multiple drug-activity prediction tasks while maintain… ▽ More We investigate Siamese networks for learning related embeddings for augmented samples of molecular conformers. We find that a non-contrastive (positive-pair only) auxiliary task aids in supervised training of Euclidean neural networks (E3NNs) and increases manifold smoothness (MS) around point-cloud geometries. We demonstrate this property for multiple drug-activity prediction tasks while maintaining relevant performance metrics, and propose an extension of MS to probabilistic and regression settings. We provide an analysis of representation collapse, finding substantial effects of task-weighting, latent dimension, and regularization. We expect the presented protocol to aid in the development of reliable E3NNs from molecular conformers, even for small-data drug discovery programs. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: Submitted to the MLDD workshop, ICLR 2023

Hierarchical Inference of the Lensing Convergence from Photometric Catalogs with Bayesian Graph Neural Networks

Authors: Ji Won Park, Simon Birrer, Madison Ueland, Miles Cranmer, Adriano Agnello, Sebastian Wagner-Carena, Philip J. Marshall, Aaron Roodman, the LSST Dark Energy Science Collaboration

Abstract: We present a Bayesian graph neural network (BGNN) that can estimate the weak lensing convergence ($κ$) from photometric measurements of galaxies along a given line of sight. The method is of particular interest in strong gravitational time delay cosmography (TDC), where characterizing the "external convergence" ($κ_{\rm ext}$) from the lens environment and line of sight is necessary for precise in… ▽ More We present a Bayesian graph neural network (BGNN) that can estimate the weak lensing convergence ($κ$) from photometric measurements of galaxies along a given line of sight. The method is of particular interest in strong gravitational time delay cosmography (TDC), where characterizing the "external convergence" ($κ_{\rm ext}$) from the lens environment and line of sight is necessary for precise inference of the Hubble constant ($H_0$). Starting from a large-scale simulation with a $κ$ resolution of $\sim$1$'$, we introduce fluctuations on galaxy-galaxy lensing scales of $\sim$1$''$ and extract random sightlines to train our BGNN. We then evaluate the model on test sets with varying degrees of overlap with the training distribution. For each test set of 1,000 sightlines, the BGNN infers the individual $κ$ posteriors, which we combine in a hierarchical Bayesian model to yield constraints on the hyperparameters governing the population. For a test field well sampled by the training set, the BGNN recovers the population mean of $κ$ precisely and without bias, resulting in a contribution to the $H_0$ error budget well under 1\%. In the tails of the training set with sparse samples, the BGNN, which can ingest all available information about each sightline, extracts more $κ$ signal compared to a simplified version of the traditional method based on matching galaxy number counts, which is limited by sample variance. Our hierarchical inference pipeline using BGNNs promises to improve the $κ_{\rm ext}$ characterization for precision TDC. The implementation of our pipeline is available as a public Python package, Node to Joy. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 15 pages, 8 figures (+ 6 pages, 2 figures in Appendix). Submitted to ApJ. Code at https://github.com/jiwoncpark/node-to-joy

PropertyDAG: Multi-objective Bayesian optimization of partially ordered, mixed-variable properties for biological sequence design

Authors: Ji Won Park, Samuel Stanton, Saeed Saremi, Andrew Watkins, Henri Dwyer, Vladimir Gligorijevic, Richard Bonneau, Stephen Ra, Kyunghyun Cho

Abstract: Bayesian optimization offers a sample-efficient framework for navigating the exploration-exploitation trade-off in the vast design space of biological sequences. Whereas it is possible to optimize the various properties of interest jointly using a multi-objective acquisition function, such as the expected hypervolume improvement (EHVI), this approach does not account for objectives with a hierarch… ▽ More Bayesian optimization offers a sample-efficient framework for navigating the exploration-exploitation trade-off in the vast design space of biological sequences. Whereas it is possible to optimize the various properties of interest jointly using a multi-objective acquisition function, such as the expected hypervolume improvement (EHVI), this approach does not account for objectives with a hierarchical dependency structure. We consider a common use case where some regions of the Pareto frontier are prioritized over others according to a specified $\textit{partial ordering}$ in the objectives. For instance, when designing antibodies, we would like to maximize the binding affinity to a target antigen only if it can be expressed in live cell culture -- modeling the experimental dependency in which affinity can only be measured for antibodies that can be expressed and thus produced in viable quantities. In general, we may want to confer a partial ordering to the properties such that each property is optimized conditioned on its parent properties satisfying some feasibility condition. To this end, we present PropertyDAG, a framework that operates on top of the traditional multi-objective BO to impose this desired ordering on the objectives, e.g. expression $\rightarrow$ affinity. We demonstrate its performance over multiple simulated active learning iterations on a penicillin production task, toy numerical problem, and a real-world antibody design task. △ Less

Submitted 8 October, 2022; originally announced October 2022.

Comments: 9 pages, 7 figures. Submitted to NeurIPS 2022 AI4Science Workshop

Multi-segment preserving sampling for deep manifold sampler

Authors: Daniel Berenberg, Jae Hyeon Lee, Simon Kelow, Ji Won Park, Andrew Watkins, Vladimir Gligorijević, Richard Bonneau, Stephen Ra, Kyunghyun Cho

Abstract: Deep generative modeling for biological sequences presents a unique challenge in reconciling the bias-variance trade-off between explicit biological insight and model flexibility. The deep manifold sampler was recently proposed as a means to iteratively sample variable-length protein sequences by exploiting the gradients from a function predictor. We introduce an alternative approach to this guide… ▽ More Deep generative modeling for biological sequences presents a unique challenge in reconciling the bias-variance trade-off between explicit biological insight and model flexibility. The deep manifold sampler was recently proposed as a means to iteratively sample variable-length protein sequences by exploiting the gradients from a function predictor. We introduce an alternative approach to this guided sampling procedure, multi-segment preserving sampling, that enables the direct inclusion of domain-specific knowledge by designating preserved and non-preserved segments along the input sequence, thereby restricting variation to only select regions. We present its effectiveness in the context of antibody design by training two models: a deep manifold sampler and a GPT-2 language model on nearly six million heavy chain sequences annotated with the IGHV1-18 gene. During sampling, we restrict variation to only the complementarity-determining region 3 (CDR3) of the input. We obtain log probability scores from a GPT-2 model for each sampled CDR3 and demonstrate that multi-segment preserving sampling generates reasonable designs while maintaining the desired, preserved regions. △ Less

Submitted 9 May, 2022; originally announced May 2022.

Z3 Charge Density Wave of Silicon Atomic Chains on a Vicinal Silicon Surface

Authors: Euihwan Do, Jae Whan Park, Oleksandr Stetsovych, Pavel Jelinek, Han Woong Yeom

Abstract: An ideal one-dimensional electronic system is formed along atomic chains on Au-decorated vicinal silicon surfaces but the nature of its low temperature phases has been puzzled for last two decades. Here, we unambiguously identify the low temperature structural distortion of this surface using high resolution atomic force microscopy and scanning tunneling microscopy. The most important structural i… ▽ More An ideal one-dimensional electronic system is formed along atomic chains on Au-decorated vicinal silicon surfaces but the nature of its low temperature phases has been puzzled for last two decades. Here, we unambiguously identify the low temperature structural distortion of this surface using high resolution atomic force microscopy and scanning tunneling microscopy. The most important structural ingredient of this surface, the step-edge Si chains are found to be strongly buckled, every third atoms down, forming trimer unitcells. This observation is consistent with the recent model of rehybridized dangling bonds and rules out the antiferromagnetic spin ordering proposed earlier. The spectroscopy and electronic structure calculation indicate a charge density wave insulator with a Z3 topology making it possible to exploit topological phases and excitations. Tunneling current was found to substantially lower the energy barrier between three degenerate CDW states, which induces a dynamically fluctuating CDW at very low temperature. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Journal ref: ACS Nano 16, 6598 (2022)

DeepZipper: A Novel Deep Learning Architecture for Lensed Supernovae Identification

Authors: Robert Morgan, B. Nord, K. Bechtol, S. J. González, E. Buckley-Geer, A. Möller, J. W. Park, A. G. Kim, S. Birrer, M. Aguena, J. Annis, S. Bocquet, D. Brooks, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, R. Cawthon, L. N. da Costa, T. M. Davis, J. De Vicente, P. Doel, I. Ferrero, D. Friedel, J. Frieman, J. García-Bellido , et al. (26 additional authors not shown)

Abstract: Large-scale astronomical surveys have the potential to capture data on large numbers of strongly gravitationally lensed supernovae (LSNe). To facilitate timely analysis and spectroscopic follow-up before the supernova fades, an LSN needs to be identified soon after it begins. To quickly identify LSNe in optical survey datasets, we designed ZipperNet, a multi-branch deep neural network that combine… ▽ More Large-scale astronomical surveys have the potential to capture data on large numbers of strongly gravitationally lensed supernovae (LSNe). To facilitate timely analysis and spectroscopic follow-up before the supernova fades, an LSN needs to be identified soon after it begins. To quickly identify LSNe in optical survey datasets, we designed ZipperNet, a multi-branch deep neural network that combines convolutional layers (traditionally used for images) with long short-term memory (LSTM) layers (traditionally used for time series). We tested ZipperNet on the task of classifying objects from four categories -- no lens, galaxy-galaxy lens, lensed type Ia supernova, lensed core-collapse supernova -- within high-fidelity simulations of three cosmic survey data sets -- the Dark Energy Survey (DES), Rubin Observatory's Legacy Survey of Space and Time (LSST), and a Dark Energy Spectroscopic Instrument (DESI) imaging survey. Among our results, we find that for the LSST-like dataset, ZipperNet classifies LSNe with a receiver operating characteristic area under the curve of 0.97, predicts the spectroscopic type of the lensed supernovae with 79\% accuracy, and demonstrates similarly high performance for LSNe 1-2 epochs after first detection. We anticipate that a model like ZipperNet, which simultaneously incorporates spatial and temporal information, can play a significant role in the rapid identification of lensed transient systems in cosmic survey experiments. △ Less

Submitted 19 May, 2022; v1 submitted 2 December, 2021; originally announced December 2021.

Comments: Published in ApJ

Report number: FERMILAB-PUB-21-392-E-SCD

Analytical Gradient Theory for Spin-Free State-Averaged Second-Order Driven Similarity Renormalization Group Perturbation Theory (SA-DSRG-MRPT2) and Its Applications for Conical Intersection Optimizations

Authors: Jae Woo Park

Abstract: The second-order multireference driven similarity renormalization group perturbation theory (DSRG-MRPT2) theory provides an efficient means of correcting the dynamical correlation with the multiconfiguration reference function. The state-averaged DSRG-MRPT2 (SA-DSRG-MRPT2) method is the simplest means of treating the excited states with DSRG-MRPT2. In this method, the Hamiltonian dressed with dyna… ▽ More The second-order multireference driven similarity renormalization group perturbation theory (DSRG-MRPT2) theory provides an efficient means of correcting the dynamical correlation with the multiconfiguration reference function. The state-averaged DSRG-MRPT2 (SA-DSRG-MRPT2) method is the simplest means of treating the excited states with DSRG-MRPT2. In this method, the Hamiltonian dressed with dynamical correlation is diagonalized in the CASCI state subspace (SA-DSRG-MRPT2c) or the configuration subspace (SA-DSRG-MRPT2). This work develops the analytical gradient theory for spin-free SA-DSRG-MRPT2(c) with the density-fitting (DF) approximation. We check the accuracy of the analytical gradients against the numerical gradients. We present applications for optimizing minimum energy conical intersections (MECI) of ethylene and retinal model chromophores (PSB3 and RPSB6). We investigate the dependence of the optimized geometries and energies on the flow parameter and reference relaxations. The smoothness of the SA-DSRG-MRPT2(c) potential energy surfaces near the reference (CASSCF) MECI is comparable to the XMCQDPT2 one. These results render the SA-DSRG-MRPT2(c) theory a promising approach for studies of conical intersections. △ Less

Submitted 28 January, 2022; v1 submitted 14 November, 2021; originally announced November 2021.

Comments: 40 pages, 8 figures, 4 tables

Creation and annihilation of mobile fractional solitons in atomic chains

Authors: Jae Whan Park, Eui Hwan Do, Jin Sung Shin, Sun Kyu Song, Oleksandr Stetsovych, Pavel Jelinek, Han Woong Yeom

Abstract: Localized modes in one dimensional topological systems, such as Majonara modes in topological superconductors, are promising platforms for robust information processing. In one dimensional topological insulators, mobile topological solitons are expected but have not been fully realized yet. We discover fractionalized phase defects moving along trimer silicon atomic chains formed along step edges o… ▽ More Localized modes in one dimensional topological systems, such as Majonara modes in topological superconductors, are promising platforms for robust information processing. In one dimensional topological insulators, mobile topological solitons are expected but have not been fully realized yet. We discover fractionalized phase defects moving along trimer silicon atomic chains formed along step edges of a vicinal silicon surface. Tunneling microscopy identifies local defects with phase shifts of 2π/3 and 4π/3 with their electronic states within the band gap and with their motions activated above 100 K. Theoretical calculations reveal the topological soliton origin of the phase defects with fractional charges of {\pm}2e/3 and {\pm}4e/3. An individual soliton can be created and annihilated at a desired location by current pulse from the probe tip. Mobile and manipulatable topological solitons discovered here provide a new platform of robustly-protected informatics with extraordinary functionalities. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: 18 pages, 4 figures

Analytical Gradient Theory for Resolvent-Fitted Second-Order Extended Multiconfiguration Perturbation Theory (XMCQDPT2)

Authors: Jae Woo Park

Abstract: We present the formulation and implementation of an analytical gradient algorithm for extended multiconfiguration quasidegenerate perturbation theory (XMCQDPT2) with the resolvent-fitting approximation by Granovsky. This algorithm is powerful when optimizing molecular configurations with a moderate-sized active space and many electronic states. First, we present the powerfulness and accuracy of re… ▽ More We present the formulation and implementation of an analytical gradient algorithm for extended multiconfiguration quasidegenerate perturbation theory (XMCQDPT2) with the resolvent-fitting approximation by Granovsky. This algorithm is powerful when optimizing molecular configurations with a moderate-sized active space and many electronic states. First, we present the powerfulness and accuracy of resolvent-fitting approximations compared to the canonical XMCQDPT2 theory. Then, we demonstrate the utility of the current algorithm in frequency analyses, optimizing the minimum energy conical intersection (MECI) geometries of the retinal chromophore model RPSB6, and evaluating nuclear gradients when there are many electronic states. Furthermore, we parallelize the algorithm using the OpenMP/MPI hybrid approach. Additionally, we report the computational cost and parallel efficiency of the program. △ Less

Submitted 27 August, 2021; v1 submitted 20 June, 2021; originally announced June 2021.

Comments: 37 pages, 4 figures, 1 scheme, 4 tables

lenstronomy II: A gravitational lensing software ecosystem

Authors: Simon Birrer, Anowar J. Shajib, Daniel Gilman, Aymeric Galan, Jelle Aalbers, Martin Millon, Robert Morgan, Giulia Pagano, Ji Won Park, Luca Teodori, Nicolas Tessore, Madison Ueland, Lyne Van de Vyvere, Sebastian Wagner-Carena, Ewoud Wempe, Lilan Yang, Xuheng Ding, Thomas Schmidt, Dominique Sluse, Ming Zhang, Adam Amara

Abstract: lenstronomy is an Astropy-affiliated Python package for gravitational lensing simulations and analyses. lenstronomy was introduced by Birrer and Amara (2018) and is based on the linear basis set approach by Birrer et a. (2015). The user and developer base of lenstronomy has substantially grown since then, and the software has become an integral part of a wide range of recent analyses, such as meas… ▽ More lenstronomy is an Astropy-affiliated Python package for gravitational lensing simulations and analyses. lenstronomy was introduced by Birrer and Amara (2018) and is based on the linear basis set approach by Birrer et a. (2015). The user and developer base of lenstronomy has substantially grown since then, and the software has become an integral part of a wide range of recent analyses, such as measuring the Hubble constant with time-delay strong lensing or constraining the nature of dark matter from resolved and unresolved small scale lensing distortion statistics. The modular design has allowed the community to incorporate innovative new methods, as well as to develop enhanced software and wrappers with more specific aims on top of the lenstronomy API. Through community engagement and involvement, lenstronomy has become a foundation of an ecosystem of affiliated packages extending the original scope of the software and proving its robustness and applicability at the forefront of the strong gravitational lensing community in an open source and reproducible manner. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Comments: published by JOSS. Software available at https://github.com/sibirrer/lenstronomy. Comments, issues and pull requests welcome!

Journal ref: Journal of Open Source Software, 6(62), 3283 (2021)

Inferring Black Hole Properties from Astronomical Multivariate Time Series with Bayesian Attentive Neural Processes

Authors: Ji Won Park, Ashley Villar, Yin Li, Yan-Fei Jiang, Shirley Ho, Joshua Yao-Yu Lin, Philip J. Marshall, Aaron Roodman

Abstract: Among the most extreme objects in the Universe, active galactic nuclei (AGN) are luminous centers of galaxies where a black hole feeds on surrounding matter. The variability patterns of the light emitted by an AGN contain information about the physical properties of the underlying black hole. Upcoming telescopes will observe over 100 million AGN in multiple broadband wavelengths, yielding a large… ▽ More Among the most extreme objects in the Universe, active galactic nuclei (AGN) are luminous centers of galaxies where a black hole feeds on surrounding matter. The variability patterns of the light emitted by an AGN contain information about the physical properties of the underlying black hole. Upcoming telescopes will observe over 100 million AGN in multiple broadband wavelengths, yielding a large sample of multivariate time series with long gaps and irregular sampling. We present a method that reconstructs the AGN time series and simultaneously infers the posterior probability density distribution (PDF) over the physical quantities of the black hole, including its mass and luminosity. We apply this method to a simulated dataset of 11,000 AGN and report precision and accuracy of 0.4 dex and 0.3 dex in the inferred black hole mass. This work is the first to address probabilistic time series reconstruction and parameter inference for AGN in an end-to-end fashion. △ Less

Submitted 18 June, 2021; v1 submitted 2 June, 2021; originally announced June 2021.

Comments: 6 pages, 4 figures, 1 table, written for non-astronomers, submitted to the ICML 2021 Time Series and Uncertainty and Robustness in Deep Learning Workshops. Comments welcome! Added affiliations and references for Fig 1

Near-Exact CASSCF-Level Geometry Optimization with a Large Active Space using Adaptive Sampling Configuration Interaction Self-Consistent Field Corrected with Second-Order Perturbation Theory (ASCI-SCF-PT2)

Authors: Jae Woo Park

Abstract: An accurate description of electron correlation is one of the most challenging problems in quantum chemistry. The exact electron correlation can be obtained by means of full configuration interaction (FCI). A simple strategy for approximating FCI at a reduced computational cost is selected CI (SCI), which diagonalizes the Hamiltonian within only the chosen configuration space. Recovery of the cont… ▽ More An accurate description of electron correlation is one of the most challenging problems in quantum chemistry. The exact electron correlation can be obtained by means of full configuration interaction (FCI). A simple strategy for approximating FCI at a reduced computational cost is selected CI (SCI), which diagonalizes the Hamiltonian within only the chosen configuration space. Recovery of the contributions of the remaining configurations is possible with second-order perturbation theory. Here, we apply adaptive sampling configuration interaction (ASCI) combined with molecular orbital optimizations (ASCI-SCF) corrected with second-order perturbation theory (ASCI-SCF-PT2) for geometry optimization by implementing the analytical nuclear gradient algorithm for ASCI-PT2 with the Z-vector (Lagrangian) formalism. We demonstrate that for phenalenyl radicals and anthracene, optimized geometries and the number of unpaired electrons can be obtained at nearly the CASSCF accuracy by incorporating PT2 corrections and extrapolating them. We demonstrate the current algorithm's utility for optimizing the equilibrium geometries and electronic structures of 6-ring-fused polycyclic aromatic hydrocarbons and 4-periacene. △ Less

Submitted 12 May, 2021; v1 submitted 17 March, 2021; originally announced March 2021.

Comments: 9 figures, 5 tables

DESC DC2 Data Release Note

Authors: LSST Dark Energy Science Collaboration, Bela Abolfathi, Robert Armstrong, Humna Awan, Yadu N. Babuji, Franz Erik Bauer, George Beckett, Rahul Biswas, Joanne R. Bogart, Dominique Boutigny, Kyle Chard, James Chiang, Johann Cohen-Tanugi, Andrew J. Connolly, Scott F. Daniel, Seth W. Digel, Alex Drlica-Wagner, Richard Dubois, Eric Gawiser, Thomas Glanzman, Salman Habib, Andrew P. Hearin, Katrin Heitmann, Fabio Hernandez, Renée Hložek , et al. (32 additional authors not shown)

Abstract: In preparation for cosmological analyses of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), the LSST Dark Energy Science Collaboration (LSST DESC) has created a 300 deg$^2$ simulated survey as part of an effort called Data Challenge 2 (DC2). The DC2 simulated sky survey, in six optical bands with observations following a reference LSST observing cadence, was processed with th… ▽ More In preparation for cosmological analyses of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), the LSST Dark Energy Science Collaboration (LSST DESC) has created a 300 deg$^2$ simulated survey as part of an effort called Data Challenge 2 (DC2). The DC2 simulated sky survey, in six optical bands with observations following a reference LSST observing cadence, was processed with the LSST Science Pipelines (19.0.0). In this Note, we describe the public data release of the resulting object catalogs for the coadded images of five years of simulated observations along with associated truth catalogs. We include a brief description of the major features of the available data sets. To enable convenient access to the data products, we have developed a web portal connected to Globus data services. We describe how to access the data and provide example Jupyter Notebooks in Python to aid first interactions with the data. We welcome feedback and questions about the data release via a GitHub repository. △ Less

Submitted 13 June, 2022; v1 submitted 12 January, 2021; originally announced January 2021.

Comments: 25 pages, 3 figures; 9 tables. A detailed changelog can be found in Appendix A. To obtain data, visit the DESC Data Portal at https://data.lsstdesc.org/

Large-Scale Gravitational Lens Modeling with Bayesian Neural Networks for Accurate and Precise Inference of the Hubble Constant

Authors: Ji Won Park, Sebastian Wagner-Carena, Simon Birrer, Philip J. Marshall, Joshua Yao-Yu Lin, Aaron Roodman

Abstract: We investigate the use of approximate Bayesian neural networks (BNNs) in modeling hundreds of time-delay gravitational lenses for Hubble constant ($H_0$) determination. Our BNN was trained on synthetic HST-quality images of strongly lensed active galactic nuclei (AGN) with lens galaxy light included. The BNN can accurately characterize the posterior PDFs of model parameters governing the elliptica… ▽ More We investigate the use of approximate Bayesian neural networks (BNNs) in modeling hundreds of time-delay gravitational lenses for Hubble constant ($H_0$) determination. Our BNN was trained on synthetic HST-quality images of strongly lensed active galactic nuclei (AGN) with lens galaxy light included. The BNN can accurately characterize the posterior PDFs of model parameters governing the elliptical power-law mass profile in an external shear field. We then propagate the BNN-inferred posterior PDFs into ensemble $H_0$ inference, using simulated time delay measurements from a plausible dedicated monitoring campaign. Assuming well-measured time delays and a reasonable set of priors on the environment of the lens, we achieve a median precision of $9.3$\% per lens in the inferred $H_0$. A simple combination of 200 test-set lenses results in a precision of 0.5 $\textrm{km s}^{-1} \textrm{ Mpc}^{-1}$ ($0.7\%$), with no detectable bias in this $H_0$ recovery test. The computation time for the entire pipeline -- including the training set generation, BNN training, and $H_0$ inference -- translates to 9 minutes per lens on average for 200 lenses and converges to 6 minutes per lens as the sample size is increased. Being fully automated and efficient, our pipeline is a promising tool for exploring ensemble-level systematics in lens modeling for $H_0$ inference. △ Less

Submitted 11 April, 2021; v1 submitted 30 November, 2020; originally announced December 2020.

Comments: 21 pages (+2 appendix), 17 figures. Published in ApJ. Code at https://github.com/jiwoncpark/h0rton. Datasets, trained models, and inference results at https://zenodo.org/record/4300382

Journal ref: ApJ 910 39 (2021)

Hierarchical Inference With Bayesian Neural Networks: An Application to Strong Gravitational Lensing

Authors: Sebastian Wagner-Carena, Ji Won Park, Simon Birrer, Philip J. Marshall, Aaron Roodman, Risa H. Wechsler

Abstract: In the past few years, approximate Bayesian Neural Networks (BNNs) have demonstrated the ability to produce statistically consistent posteriors on a wide range of inference problems at unprecedented speed and scale. However, any disconnect between training sets and the distribution of real-world objects can introduce bias when BNNs are applied to data. This is a common challenge in astrophysics an… ▽ More In the past few years, approximate Bayesian Neural Networks (BNNs) have demonstrated the ability to produce statistically consistent posteriors on a wide range of inference problems at unprecedented speed and scale. However, any disconnect between training sets and the distribution of real-world objects can introduce bias when BNNs are applied to data. This is a common challenge in astrophysics and cosmology, where the unknown distribution of objects in our Universe is often the science goal. In this work, we incorporate BNNs with flexible posterior parameterizations into a hierarchical inference framework that allows for the reconstruction of population hyperparameters and removes the bias introduced by the training distribution. We focus on the challenge of producing posterior PDFs for strong gravitational lens mass model parameters given Hubble Space Telescope (HST) quality single-filter, lens-subtracted, synthetic imaging data. We show that the posterior PDFs are sufficiently accurate (i.e., statistically consistent with the truth) across a wide variety of power-law elliptical lens mass distributions. We then apply our approach to test data sets whose lens parameters are drawn from distributions that are drastically different from the training set. We show that our hierarchical inference framework mitigates the bias introduced by an unrepresentative training set's interim prior. Simultaneously, given a sufficiently broad training set, we can precisely reconstruct the population hyperparameters governing our test distributions. Our full pipeline, from training to hierarchical inference on thousands of lenses, can be run in a day. The framework presented here will allow us to efficiently exploit the full constraining power of future ground- and space-based surveys. △ Less

Submitted 22 March, 2021; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: Accepted by ApJ. Code available at https://github.com/swagnercarena/ovejero

Journal ref: ApJ 909 187 (2021)

The LSST DESC DC2 Simulated Sky Survey

Authors: LSST Dark Energy Science Collaboration, Bela Abolfathi, David Alonso, Robert Armstrong, Éric Aubourg, Humna Awan, Yadu N. Babuji, Franz Erik Bauer, Rachel Bean, George Beckett, Rahul Biswas, Joanne R. Bogart, Dominique Boutigny, Kyle Chard, James Chiang, Chuck F. Claver, Johann Cohen-Tanugi, Céline Combet, Andrew J. Connolly, Scott F. Daniel, Seth W. Digel, Alex Drlica-Wagner, Richard Dubois, Emmanuel Gangler, Eric Gawiser , et al. (55 additional authors not shown)

Abstract: We describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses… ▽ More We describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep (WFD) area of approximately 300 deg^2 as well as a deep drilling field (DDF) of approximately 1 deg^2. We simulate 5 years of the planned 10-year survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the dataset to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic testbed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time-domain cosmology. △ Less

Submitted 26 January, 2021; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: 39 pages, 19 figures, version accepted for publication in ApJS

Atomic structures and electronic correlation of monolayer 1T-TaSe2

Authors: Jae Whan Park, Han Woong Yeom

Abstract: We investigate atomic and electronic structures of monolayer 1T-TaSe2 using density functional theory calculations. Monolayers of 1T-TaSe2 were recently grown on graphene substrates and suggested as an intriguing Mott insulator [Nat. Phys. 16, 218 (2020)]. However, the prevailing structural model for the model system of 1T-TaS2, the cation-centered cluster of a David-star shape with strong electro… ▽ More We investigate atomic and electronic structures of monolayer 1T-TaSe2 using density functional theory calculations. Monolayers of 1T-TaSe2 were recently grown on graphene substrates and suggested as an intriguing Mott insulator [Nat. Phys. 16, 218 (2020)]. However, the prevailing structural model for the model system of 1T-TaS2, the cation-centered cluster of a David-star shape with strong electron correlation, could not explain the characteristic and unusual orbital splitting observed in scanning tunneling spectroscopy experiments. We suggest an alternative structure model, an anion-centered cluster structure, which can reproduce most of the unusual spectroscopic characteristics with electron doping from the substrate without electron correlation. The unusual spectroscopic features observed, thus, seems to indicate a simple and usual band insulating state. This work indicates the importance of a large structural degree of freedom given for a cluster Mott insulator. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: 4 page, 3 figures

Network reinforcement driven drug repurposing for COVID-19 by exploiting disease-gene-drug associations

Authors: Yonghyun Nam, Jae-Seung Yun, Seung Mi Lee, Ji Won Park, Ziqi Chen, Brian Lee, Anurag Verma, Xia Ning, Li Shen, Dokyoon Kim

Abstract: Currently, the number of patients with COVID-19 has significantly increased. Thus, there is an urgent need for developing treatments for COVID-19. Drug repurposing, which is the process of reusing already-approved drugs for new medical conditions, can be a good way to solve this problem quickly and broadly. Many clinical trials for COVID-19 patients using treatments for other diseases have already… ▽ More Currently, the number of patients with COVID-19 has significantly increased. Thus, there is an urgent need for developing treatments for COVID-19. Drug repurposing, which is the process of reusing already-approved drugs for new medical conditions, can be a good way to solve this problem quickly and broadly. Many clinical trials for COVID-19 patients using treatments for other diseases have already been in place or will be performed at clinical sites in the near future. Additionally, patients with comorbidities such as diabetes mellitus, obesity, liver cirrhosis, kidney diseases, hypertension, and asthma are at higher risk for severe illness from COVID-19. Thus, the relationship of comorbidity disease with COVID-19 may help to find repurposable drugs. To reduce trial and error in finding treatments for COVID-19, we propose building a network-based drug repurposing framework to prioritize repurposable drugs. First, we utilized knowledge of COVID-19 to construct a disease-gene-drug network (DGDr-Net) representing a COVID-19-centric interactome with components for diseases, genes, and drugs. DGDr-Net consisted of 592 diseases, 26,681 human genes and 2,173 drugs, and medical information for 18 common comorbidities. The DGDr-Net recommended candidate repurposable drugs for COVID-19 through network reinforcement driven scoring algorithms. The scoring algorithms determined the priority of recommendations by utilizing graph-based semi-supervised learning. From the predicted scores, we recommended 30 drugs, including dexamethasone, resveratrol, methotrexate, indomethacin, quercetin, etc., as repurposable drugs for COVID-19, and the results were verified with drugs that have been under clinical trials. The list of drugs via a data-driven computational approach could help reduce trial-and-error in finding treatment for COVID-19. △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: 4 figures

Continual BERT: Continual Learning for Adaptive Extractive Summarization of COVID-19 Literature

Authors: Jong Won Park

Abstract: The scientific community continues to publish an overwhelming amount of new research related to COVID-19 on a daily basis, leading to much literature without little to no attention. To aid the community in understanding the rapidly flowing array of COVID-19 literature, we propose a novel BERT architecture that provides a brief yet original summarization of lengthy papers. The model continually lea… ▽ More The scientific community continues to publish an overwhelming amount of new research related to COVID-19 on a daily basis, leading to much literature without little to no attention. To aid the community in understanding the rapidly flowing array of COVID-19 literature, we propose a novel BERT architecture that provides a brief yet original summarization of lengthy papers. The model continually learns on new data in online fashion while minimizing catastrophic forgetting, thus fitting to the need of the community. Benchmark and manual examination of its performance show that the model provide a sound summary of new scientific literature. △ Less

Submitted 27 September, 2020; v1 submitted 7 July, 2020; originally announced July 2020.

Comments: 6 pages, 3 figures

Time Delay Lens Modelling Challenge

Authors: X. Ding, T. Treu, S. Birrer, G. C. -F. Chen, J. Coles, P. Denzel, M. Frigo A. Galan, P. J. Marshall, M. Millon, A. More, A. J. Shajib, D. Sluse, H. Tak, D. Xu, M. W. Auger, V. Bonvin, H. Chand, F. Courbin, G. Despali, C. D. Fassnacht, D. Gilman, S. Hilbert, S. R. Kumar, Y. -Y. Lin, J. W. Park , et al. (4 additional authors not shown)

Abstract: In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic… ▽ More In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated datasets. The results in Rung 1 and Rung 2 show that methods that use only the point source positions tend to have lower precision ($10 - 20\%$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic datasets can recover $H_0$ within the target accuracy ($ |A| < 2\%$) and precision ($< 6\%$ per system), even in the presence of poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix, and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges. △ Less

Submitted 17 February, 2021; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: 23 pages, 12 figures, 6 tables, MNRAS accepted

Journal ref: MNRAS, (2021), 503, 1096

Resonant dipolar collisions of ultracold molecules induced by microwave dressing

Authors: Zoe Z. Yan, Jee Woo Park, Yiqi Ni, Huanqian Loh, Sebastian Will, Tijs Karman, Martin Zwierlein

Abstract: We demonstrate microwave dressing on ultracold, fermionic ${}^{23}$Na${}^{40}$K ground-state molecules and observe resonant dipolar collisions with cross sections exceeding three times the $s$-wave unitarity limit. The origin of these collisions is the resonant alignment of the approaching molecules' dipoles along the intermolecular axis, which leads to strong attraction. We explain our observatio… ▽ More We demonstrate microwave dressing on ultracold, fermionic ${}^{23}$Na${}^{40}$K ground-state molecules and observe resonant dipolar collisions with cross sections exceeding three times the $s$-wave unitarity limit. The origin of these collisions is the resonant alignment of the approaching molecules' dipoles along the intermolecular axis, which leads to strong attraction. We explain our observations with a conceptually simple two-state picture based on the Condon approximation. Furthermore, we perform coupled-channels calculations that agree well with the experimentally observed collision rates. While collisions are observed here as laser-induced loss, microwave dressing on chemically stable molecules trapped in box potentials may enable the creation of strongly interacting dipolar gases of molecules. △ Less

Submitted 5 March, 2020; originally announced March 2020.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. Lett. 125, 063401 (2020)

Artificial Relativistic Molecules

Authors: Jae Whan Park, Hyo Sung Kim, Thomas Brumme, Thomas Heine, Han Woong Yeom

Abstract: We fabricate artificial molecules composed of heavy atom lead on a van der Waals crystal. Pb atoms templated on a honeycomb charge-order superstructure of IrTe2 form clusters ranging from dimers to heptamers including benzene-shaped ring hexamers. Tunneling spectroscopy and electronic structure calculations reveal the formation of unusual relativistic molecular orbitals within the clusters. The sp… ▽ More We fabricate artificial molecules composed of heavy atom lead on a van der Waals crystal. Pb atoms templated on a honeycomb charge-order superstructure of IrTe2 form clusters ranging from dimers to heptamers including benzene-shaped ring hexamers. Tunneling spectroscopy and electronic structure calculations reveal the formation of unusual relativistic molecular orbitals within the clusters. The spin-orbit coupling is essential both in forming such Dirac electronic states and stabilizing the artificial molecules by reducing the adatom-substrate interaction. Lead atoms are found to be ideally suited for a maximized relativistic effect. This work initiates the use of novel two dimensional orderings to guide the fabrication of artificial molecules of unprecedented properties. △ Less

Submitted 3 February, 2020; originally announced February 2020.

Comments: 7 page, 3 figures

Multireference electron correlation methods: Journeys along potential energy surfaces

Authors: Jae Woo Park, Rachael Al-Saadon, Matthew K. MacLeod, Toru Shiozaki, Bess Vlaisavljevich

Abstract: Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way, and therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field (MCSCF) theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces (PES). This incl… ▽ More Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way, and therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field (MCSCF) theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces (PES). This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections, and on-the-fly photodynamics simulations; both depend heavily on the ability of the method to properly explore the PES. Since such applications require the nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly improves the utility of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. To motivate the readers, we first summarize the notable applications of multireference electron correlation methods to mainstream chemistry, including geometry optimizations and on-the-fly dynamics. Subsequently, we review the analytical nuclear gradient and derivative coupling theories for these methods, and the software infrastructure that allows one to make use of these quantities in applications. The future prospects are discussed at the end of this review. △ Less

Submitted 8 August, 2019; originally announced November 2019.

Comments: Submitted to Chem. Rev. Comments welcome

Analytical Gradient Theory for Strongly Contracted (SC-) and Partially Contracted (PC-) N-Electron Valence State Perturbation Theory (NEVPT2)

Authors: Jae Woo Park

Abstract: An analytical gradient theory for single-state N-electron valence state perturbation theory (NEVPT2), using both strongly contracted (SC) and partially contracted (PC) internal contraction schemes, is developed. We demonstrate the utility of the developed algorithm in the optimization of the single-state molecular geometry of acrolein, benzyne, benzene, the retinal chromophore PSB3, the GFP chromo… ▽ More An analytical gradient theory for single-state N-electron valence state perturbation theory (NEVPT2), using both strongly contracted (SC) and partially contracted (PC) internal contraction schemes, is developed. We demonstrate the utility of the developed algorithm in the optimization of the single-state molecular geometry of acrolein, benzyne, benzene, the retinal chromophore PSB3, the GFP chromophore pHBI, and porphine, with the cc-pVTZ basis sets. The SC-NEVPT2 analytical gradients exhibit numerical instability due to the lack of invariance with respect to the rotations among the inactive orbitals. On the other hand, PC-NEVPT2 gives molecular geometries comparable to CASPT2 in any tested cases. We discuss possible future developments that will make the NEVPT2 gradient algorithm a powerful tool for optimizing the molecular geometries and conducting molecular dynamics simulations of correlated systems. △ Less

Submitted 1 September, 2019; v1 submitted 23 July, 2019; originally announced July 2019.

Comments: 32 pages, 2 figures, 3 tables

Journal ref: J. Chem. Theory Comput. 15, 5417 (2019)

Emergent Honeycomb Network of Topological Excitations in Correlated Charge Density Wave

Authors: Jae Whan Park, Gil Young Cho, Jinwon Lee, Han Woong Yeom

Abstract: When two periodic potentials compete in materials, one may adopt the other, which straightforwardly generates topological defects. Of particular interest are domain walls in charge-, dipole-, and spin-ordered systems, which govern macroscopic properties and important functionality. However, detailed atomic and electronic structures of domain walls have often been uncertain and the microscopic mech… ▽ More When two periodic potentials compete in materials, one may adopt the other, which straightforwardly generates topological defects. Of particular interest are domain walls in charge-, dipole-, and spin-ordered systems, which govern macroscopic properties and important functionality. However, detailed atomic and electronic structures of domain walls have often been uncertain and the microscopic mechanism of their functionality has been elusive. Here, we clarify the complete atomic and electronic structures of the domain wall network, a honeycomb network connected by Z$_{3}$ vortices, in the nearly commensurate Mott charge-density wave (CDW) phase of 1T-TaS$_{2}$. Scanning tunneling microscopy resolves characteristic charge orders within domain walls and their vortices. Density functional theory calculations disclose their unique atomic relaxations and the metallic in-gap states confined tightly therein. A generic theory is constructed, which connects this emergent honeycomb network of conducting electrons to the enhanced superconductivity. △ Less

Submitted 21 July, 2019; originally announced July 2019.

Comments: 5 figures

Stable Flatbands, Topology, and Superconductivity of Magic Honeycomb Networks

Authors: Jongjun M. Lee, Chenhua Geng, Jae Whan Park, Masaki Oshikawa, Sung-Sik Lee, Han Woong Yeom, Gil Young Cho

Abstract: We propose a new principle to realize flatbands which are robust in real materials, based on a network superstructure of one-dimensional segments. This mechanism is naturally realized in the nearly commensurate charge-density wave of 1T-TaS${}_2$ with the honeycomb network of conducting domain walls, and the resulting flatband can naturally explain the enhanced superconductivity. We also show that… ▽ More We propose a new principle to realize flatbands which are robust in real materials, based on a network superstructure of one-dimensional segments. This mechanism is naturally realized in the nearly commensurate charge-density wave of 1T-TaS${}_2$ with the honeycomb network of conducting domain walls, and the resulting flatband can naturally explain the enhanced superconductivity. We also show that corner states, which are a hallmark of the higher-order topological insulators, appear in the network superstructure. △ Less

Submitted 27 March, 2020; v1 submitted 28 June, 2019; originally announced July 2019.

Comments: 7+16 pages, 2+14 figures, Accepted version to PRL

Journal ref: Phys. Rev. Lett. 124, 137002 (2020)

Imaginary shift in CASPT2 nuclear gradient and derivative coupling theory

Authors: Jae Woo Park, Rachael Al-Saadon, Nils E. Strand, Toru Shiozaki

Abstract: We report the analytical nuclear gradient theory for complete active space second-order perturbation theory (CASPT2) with imaginary shift, which is commonly used to avoid divergence of the perturbation expression. Our formulation is based on the Lagrangian approach and is an extension of the algorithm for CASPT2 nuclear gradients with real shift. The working equations are derived and implemented i… ▽ More We report the analytical nuclear gradient theory for complete active space second-order perturbation theory (CASPT2) with imaginary shift, which is commonly used to avoid divergence of the perturbation expression. Our formulation is based on the Lagrangian approach and is an extension of the algorithm for CASPT2 nuclear gradients with real shift. The working equations are derived and implemented into an efficient parallel program. Numerical examples are presented for the ground- and excited-state geometries and conical intersections of a green fluorescent protein model chromophore, $p$-HBDI$^-$. We also report timing benchmarks with adenine, $p$-HBDI$^-$, and iron porphyrin. It is demonstrated that the energies and geometries obtained with the imaginary shift improve accuracy at a minor additional cost which is mainly associated with evaluating the effective density matrix elements for the imaginary shift term. △ Less

Submitted 14 April, 2019; originally announced April 2019.

Pseudogap and weak multifractality in disordered Mott charge-density-wave insulator

Authors: Jianhua Gao, Jae Whan Park, Kiseok Kim, Sun Kyu Song, Fangchu Chen, Xuan Luo, Yuping Sun, Han Woong Yeom

Abstract: The competition, coexistence and cooperation of various orders in low-dimensional materials like spin, charge, topological orders and charge-density-wave has been one of the most intriguing issues in condensed matter physics. In particular, layered transition metal dichalcogenides provide an ideal platform for studying such an interplay with a notable case of 1${T}$-TaS$_{2}$ featuring Mott-insula… ▽ More The competition, coexistence and cooperation of various orders in low-dimensional materials like spin, charge, topological orders and charge-density-wave has been one of the most intriguing issues in condensed matter physics. In particular, layered transition metal dichalcogenides provide an ideal platform for studying such an interplay with a notable case of 1${T}$-TaS$_{2}$ featuring Mott-insulating ground state, charge-density-wave, spin frustration and emerging superconductivity together. We investigated local electronic states of Se-substituted 1${T}$-TaS$_{2}$ by scanning tunneling microscopy/spectroscopy (STM/STS), where superconductivity emerges from the unique Mott-CDW state. Spatially resolved STS measurements reveal that an apparent V-shape pseudogap forms at the Fermi Level (E$_{F}$), with the origin of the electronic states splitting and transformation from the Mott states, and the CDW gaps are largely preserved. The formation of the pseudogap has little correlation to the variation of local Se concentration, but appears to be a global characteristics. Furthermore, the correlation length of local density of states (LDOS) diverges at the Fermi energy and decays rapidly at high energies. The spatial correlation shows a power-law decay close to the Fermi energy. Our statistics analysis of the LDOS indicates that our system exhibits weak multifractal behavior of the wave functions. These findings strongly support a correlated metallic state induced by disorder in our system, which provides an new insight into the novel mechanism of emerging superconductivity in the two-dimensional correlated electronic systems. △ Less

Submitted 9 April, 2019; originally announced April 2019.

Discovering Options for Exploration by Minimizing Cover Time

Authors: Yuu Jinnai, Jee Won Park, David Abel, George Konidaris

Abstract: One of the main challenges in reinforcement learning is solving tasks with sparse reward. We show that the difficulty of discovering a distant rewarding state in an MDP is bounded by the expected cover time of a random walk over the graph induced by the MDP's transition dynamics. We therefore propose to accelerate exploration by constructing options that minimize cover time. The proposed algorithm… ▽ More One of the main challenges in reinforcement learning is solving tasks with sparse reward. We show that the difficulty of discovering a distant rewarding state in an MDP is bounded by the expected cover time of a random walk over the graph induced by the MDP's transition dynamics. We therefore propose to accelerate exploration by constructing options that minimize cover time. The proposed algorithm finds an option which provably diminishes the expected number of steps to visit every state in the state space by a uniform random walk. We show empirically that the proposed algorithm improves the learning time in several domains with sparse rewards. △ Less

Submitted 16 March, 2019; v1 submitted 1 March, 2019; originally announced March 2019.

Kibble-Zurek universality in a strongly interacting Fermi superfluid

Authors: Bumsuk Ko, Jee Woo Park, Y. Shin

Abstract: Near a continuous phase transition, systems with different microscopic origins display universal dynamics if their underlying symmetries are compatible. In a thermally quenched system, the Kibble-Zurek mechanism for the creation of topological defects unveils this universality through a characteristic power-law exponent, which captures the dependence of the defect density on the quench rate. Here,… ▽ More Near a continuous phase transition, systems with different microscopic origins display universal dynamics if their underlying symmetries are compatible. In a thermally quenched system, the Kibble-Zurek mechanism for the creation of topological defects unveils this universality through a characteristic power-law exponent, which captures the dependence of the defect density on the quench rate. Here, we report the observation of the Kibble-Zurek universality in a strongly interacting Fermi superfluid. As the system's microscopic description is tuned from bosonic to fermionic, the quench formation of vortices reveals a constant scaling exponent arising from the $U(1)$ gauge symmetry of the system. For rapid quenches, destructive vortex collisions lead to the saturation of their densities, whose values can be universally scaled by the interaction-dependent area of the vortex cores. △ Less

Submitted 19 February, 2019; originally announced February 2019.

Journal ref: Nat. Phys. (2019)

Recent results on Choi's orthogonal Latin squares

Authors: Jon-Lark Kim, Dong Eun Ohk, Doo Young Park, Jae Woo Park

Abstract: Choi Seok-Jeong studied Latin squares at least 60 years earlier than Euler although this was less known. He introduced a pair of orthogonal Latin squares of order 9 in his book. Interestingly, his two orthogonal non-double-diagonal Latin squares produce a magic square of order 9, whose theoretical reason was not studied. There have been a few studies on Choi's Latin squares of order 9. The most re… ▽ More Choi Seok-Jeong studied Latin squares at least 60 years earlier than Euler although this was less known. He introduced a pair of orthogonal Latin squares of order 9 in his book. Interestingly, his two orthogonal non-double-diagonal Latin squares produce a magic square of order 9, whose theoretical reason was not studied. There have been a few studies on Choi's Latin squares of order 9. The most recent one is Ko-Wei Lih's construction of Choi's Latin squares of order 9 based on the two $3 \times 3$ orthogonal Latin squares. In this paper, we give a new generalization of Choi's orthogonal Latin squares of order 9 to orthogonal Latin squares of size $n^2$ using the Kronecker product including Lih's construction. We find a geometric description of Choi's orthogonal Latin squares of order 9 using the dihedral group $D_8$. We also give a new way to construct magic squares from two orthogonal non-double-diagonal Latin squares, which explains why Choi's Latin squares produce a magic square of order 9. △ Less

Submitted 28 September, 2021; v1 submitted 5 December, 2018; originally announced December 2018.

Comments: 1. 16 pages revised on Sep. 8, 2021. 2. Journal of Algebra Combinatorics Discrete Structures and Applications, Accepted on Sep. 15, 2021

MSC Class: 05B15; 05B20

Critical Vortex Shedding in a Strongly Interacting Fermionic Superfluid

Authors: Jee Woo Park, Bumsuk Ko, Yong-il Shin

Abstract: We study the critical vortex shedding in a strongly interacting fermionic superfluid of $^{6}$Li across the BEC-BCS crossover. By moving an optical obstacle in the sample and directly imaging the vortices after time of flight, the critical velocity $v_{\rm c}$ for vortex shedding is measured as a function of the obstacle travel distance $L$. The observed $v_{\rm c}$ increases with decreasing $L$,… ▽ More We study the critical vortex shedding in a strongly interacting fermionic superfluid of $^{6}$Li across the BEC-BCS crossover. By moving an optical obstacle in the sample and directly imaging the vortices after time of flight, the critical velocity $v_{\rm c}$ for vortex shedding is measured as a function of the obstacle travel distance $L$. The observed $v_{\rm c}$ increases with decreasing $L$, where the rate of increase is the highest in the unitary regime. In the deep BEC regime, an empirical dissipation model well captures the dependence of $v_{\rm c}$ on $L$, characterized by a constant value of $η= -\frac{{\rm d}(1/v_{\rm c})}{{\rm d}(1/L)}$. However, as the system is tuned across the resonance, a step increase of $η$ develops about a characteristic distance $L_{\rm c}$ as $L$ is increased, where $L_{\rm c}$ is comparable to the obstacle size. This bimodal behavior is strengthened as the system is tuned towards the BCS regime. We attribute this evolution of $v_{\rm c}$ to the participation of pair breaking in the vortex shedding dynamics of a fermionic superfluid. △ Less

Submitted 22 May, 2018; originally announced May 2018.

Comments: 8 pages, 8 figures

Journal ref: Phys. Rev. Lett. 121, 225301 (2018)

On the accuracy of retinal protonated Schiff base models

Authors: Jae Woo Park, Toru Shiozaki

Abstract: We investigate the molecular geometries of the ground state and the minimal energy conical intersections (MECIs) between the ground and first excited states of the models for the retinal protonated Schiff base in the gas phase using the extended multistate complete active space second-order perturbation theory (XMS-CASPT2). The biggest model in this work is the rhodopsin chromophore truncated betw… ▽ More We investigate the molecular geometries of the ground state and the minimal energy conical intersections (MECIs) between the ground and first excited states of the models for the retinal protonated Schiff base in the gas phase using the extended multistate complete active space second-order perturbation theory (XMS-CASPT2). The biggest model in this work is the rhodopsin chromophore truncated between the ε and δ carbon atoms, which consists of 54 atoms and 12-orbital π conjugation. The results are compared with those obtained by the state-averaged complete active space self-consistent field (SA-CASSCF). The XMS-CASPT2 results suggest that the minimum energy conical intersection associated with the so-called 13-14 isomerization is thermally inaccessible, which is in contrast to the SA-CASSCF results. The differences between the geometries of the conical intersections computed by SA-CASSCF and XMS-CASPT2 are ascribed to the fact that the charge transfer states are more stabilized by dynamical electron correlation than the diradicaloid states. The impact of the various choices of active spaces, basis sets, and state averaging schemes is also examined. △ Less

Submitted 31 January, 2018; originally announced February 2018.

Comments: Contribution to the special issue in honor of the 80th birthday of Professor Michael Baer

On-the-fly CASPT2 surface hopping dynamics

Authors: Jae Woo Park, Toru Shiozaki

Abstract: We report the development of programs for on-the-fly surface hopping dynamics simulations in the gas and condensed phases on the potential energy surfaces computed by multistate multireference perturbation theory (XMS-CASPT2) with full internal contraction. On-the-fly nonadiabatic dynamics simulations are made possible by improving the algorithm for XMS-CASPT2 nuclear energy gradient and derivativ… ▽ More We report the development of programs for on-the-fly surface hopping dynamics simulations in the gas and condensed phases on the potential energy surfaces computed by multistate multireference perturbation theory (XMS-CASPT2) with full internal contraction. On-the-fly nonadiabatic dynamics simulations are made possible by improving the algorithm for XMS-CASPT2 nuclear energy gradient and derivative coupling evaluation. The program is interfaced to a surface hopping dynamics program, Newton-X, and a classical molecular dynamics package, tinker, to realize such simulations. On-the-fly XMS-CASPT2 surface-hopping dynamics simulations of 9H-adenine and an anionic GFP model chromophore (para-hydroxybenzilideneimidazolin-5-one) in water are presented to demonstrate the applicability of our program to sizable systems. Our program is implemented in the bagel package, which is publicly available under the GNU General Public License. △ Less

Submitted 29 June, 2017; v1 submitted 31 May, 2017; originally announced June 2017.

Analytical derivative coupling for multistate CASPT2 theory

Authors: Jae Woo Park, Toru Shiozaki

Abstract: The probability of non-radiative transitions in photochemical dynamics is determined by the derivative couplings, the couplings between different electronic states through the nuclear degrees of freedom. Efficient and accurate evaluation of the derivative couplings is, therefore, of central importance to realize reliable computer simulations of photochemical reactions. In this work, the derivative… ▽ More The probability of non-radiative transitions in photochemical dynamics is determined by the derivative couplings, the couplings between different electronic states through the nuclear degrees of freedom. Efficient and accurate evaluation of the derivative couplings is, therefore, of central importance to realize reliable computer simulations of photochemical reactions. In this work, the derivative couplings for multistate multireference second-order perturbation theory (MS-CASPT2) and its 'extended' variant (XMS-CASPT2) are studied, in which we present an algorithm for their analytical evaluation. The computational costs for evaluating the derivative couplings are essentially the same as those for calculating the nuclear energy gradients. The geometries and energies calculated with XMS-CASPT2 for small molecules at minimum energy conical intersections (MECIs) are in good agreement with those computed by multireference configuration interaction. As numerical examples, MECIs are optimized using XMS-CASPT2 for stilbene and a GFP model chromophore (the 4-para-hydroxybenzylidene-1,2-dimethyl-imidazolin-5-one anion). △ Less

Submitted 3 May, 2017; v1 submitted 9 January, 2017; originally announced January 2017.

Second-Scale Nuclear Spin Coherence Time of Trapped Ultracold $^{23}$Na$^{40}$K Molecules

Authors: Jee Woo Park, Zoe Z. Yan, Huanqian Loh, Sebastian A. Will, Martin W. Zwierlein

Abstract: Coherence, the stability of the relative phase between quantum states, lies at the heart of quantum mechanics. Applications such as precision measurement, interferometry, and quantum computation are enabled by physical systems that have quantum states with robust coherence. With the creation of molecular ensembles at sub-$μ$K temperatures, diatomic molecules have become a novel system under full q… ▽ More Coherence, the stability of the relative phase between quantum states, lies at the heart of quantum mechanics. Applications such as precision measurement, interferometry, and quantum computation are enabled by physical systems that have quantum states with robust coherence. With the creation of molecular ensembles at sub-$μ$K temperatures, diatomic molecules have become a novel system under full quantum control. Here, we report on the observation of stable coherence between a pair of nuclear spin states of ultracold fermionic NaK molecules in the singlet rovibrational ground state. Employing microwave fields, we perform Ramsey spectroscopy and observe coherence times on the scale of one second. This work opens the door for the exploration of single molecules as a versatile quantum memory. Switchable long-range interactions between dipolar molecules can further enable two-qubit gates, allowing quantum storage and processing in the same physical system. Within the observed coherence time, $10^4$ one- and two-qubit gate operations will be feasible. △ Less

Submitted 13 June, 2016; originally announced June 2016.

Comments: 8 pages, 5 figures

Journal ref: Science 357, 372-375 (2017)

Metallic Indium Monolayers on Si(111)

Authors: Jae Whan Park, Myung Ho Kang

Abstract: Density-functional calculations are used to identify one-atom-thick metallic In overlayers on the Si(111) surface, which have long been sought in quest of the ultimate two-dimensional (2D) limit of free-electron-like metallic properties. We predict two metastable single-layer In phases, one $\sqrt{7}\times\sqrt{3}$ phase with a coverage of 1.4 monolayer (ML; here 1 ML refers to one In atom per top… ▽ More Density-functional calculations are used to identify one-atom-thick metallic In overlayers on the Si(111) surface, which have long been sought in quest of the ultimate two-dimensional (2D) limit of free-electron-like metallic properties. We predict two metastable single-layer In phases, one $\sqrt{7}\times\sqrt{3}$ phase with a coverage of 1.4 monolayer (ML; here 1 ML refers to one In atom per top Si atom) and the other $\sqrt{7}\times\sqrt{7}$ phase with 1.43 ML, which indeed match well with experimental evidences. Both phases reveal quasi-1D arrangements of protruded In atoms, leading to 2D-metallic but anisotropic band structures and Fermi surfaces. This directional feature contrasts with the free-electron-like In-overlayer properties that are known to persist up to the double-layer thickness, implying that we may have achieved the 2D limit of free-electron-like In overlayers in previous studies of double-layer In phases. △ Less

Submitted 27 May, 2016; originally announced May 2016.

Comments: 5 pages, 6 figures

Journal ref: Phys. Rev. Lett. 117, 116102 (2016)

Coherent Microwave Control of Ultracold $^{23}$Na$^{40}$K Molecules

Authors: Sebastian A. Will, Jee Woo Park, Zoe Z. Yan, Huanqian Loh, Martin W. Zwierlein

Abstract: We demonstrate coherent microwave control of rotational and hyperfine states of trapped, ultracold, and chemically stable $^{23}$Na$^{40}$K molecules. Starting with all molecules in the absolute rovibrational and hyperfine ground state, we study rotational transitions in combined magnetic and electric fields and explain the rich hyperfine structure. Following the transfer of the entire molecular e… ▽ More We demonstrate coherent microwave control of rotational and hyperfine states of trapped, ultracold, and chemically stable $^{23}$Na$^{40}$K molecules. Starting with all molecules in the absolute rovibrational and hyperfine ground state, we study rotational transitions in combined magnetic and electric fields and explain the rich hyperfine structure. Following the transfer of the entire molecular ensemble into a single hyperfine level of the first rotationally excited state, $J{=}1$, we observe collisional lifetimes of more than $3\, \rm s$, comparable to those in the rovibrational ground state, $J{=}0$. Long-lived ensembles and full quantum state control are prerequisites for the use of ultracold molecules in quantum simulation, precision measurements and quantum information processing. △ Less

Submitted 31 March, 2016; originally announced April 2016.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. Lett. 116, 225306 (2016)