Showing 1–50 of 213 results for author: Chung, E

Correlation between dust continuum and CN line emissions in high-mass star-forming regions

Authors: Jihye Hwang, Chang Won Lee, Jongsoo Kim, Eun Jung Chung, Kee-Tae Kim

Abstract: Measuring the strength of three dimensional (3D) magnetic field vector is challenging as it is not easy to recognize whether its line-of-sight (LOS) and plane-of-sky (POS) components are obtained from the same region. CN ($N = 1 - 0$) emission has been used to get the LOS component of a magnetic field (B$_\mathrm{LOS}$) from its Zeeman splitting lines, while dust continuum emission has been used t… ▽ More Measuring the strength of three dimensional (3D) magnetic field vector is challenging as it is not easy to recognize whether its line-of-sight (LOS) and plane-of-sky (POS) components are obtained from the same region. CN ($N = 1 - 0$) emission has been used to get the LOS component of a magnetic field (B$_\mathrm{LOS}$) from its Zeeman splitting lines, while dust continuum emission has been used to get the POS component of a magnetic field (B$_\mathrm{POS}$). We use the CN ($N = 1 - 0$) data observed with the Taeduk Radio Astronomy Observatory (TRAO) 14-m telescope and the dust continuum data from $Herschel$ archive toward six high-mass star-forming regions in order to test whether CN line and dust continuum emission can trace a similar region and thus can be used for inferring 3D magnetic field strength. Our comparison between CN and H$_2$ column densities for all targets indicates that CN line emission tends to be strong toward bright continuum regions. The positions of peak CN column densities are particularly well correlated with those of peak H$_2$ column densities at least over the H$_2$ column density of 8.0 $\times$ 10$^{22}$ cm$^{-2}$ within one or two telescope beam size in all targets, implying that CN line and dust continuum emitting regions are likely spatially coincident. This enabled us to make the reliable measurement of 3D magnetic field strengths of five targets by taking a vector sum of their B$_\mathrm{LOS}$ and B$_\mathrm{POS}$, helping to decide the magnetical criticality of the targets as supercritical or transcritical. △ Less

Submitted 19 August, 2024; originally announced August 2024.

Comments: Accepted in ApJ, 14th Aug 2024

Magnetic Fields in Massive Star-forming Regions (MagMaR) IV: Tracing the Magnetic Fields in the O-type protostellar system IRAS 16547$-$4247

Authors: Luis A. Zapata, Manuel Fernández-López, Patricio Sanhueza, Josep M. Girart, Luis F. Rodríguez, Paulo Cortes, Koch Patrick, María T. Beltrán, Kate Pattle, Henrik Beuther, Piyali Saha, Wenyu Jiao, Fengwei Xu, Xing Walker Lu, Fernando Olguin, Shanghuo Li, Ian W. Stephens, Ji-hyun Kang, Yu Cheng, Spandan Choudhury, Kaho Morii, Eun Jung Chung, Jia-Wei Wang, Jihye Hwang, A-Ran Lyo , et al. (2 additional authors not shown)

Abstract: The formation of the massive stars, and in particular, the role that the magnetic fields play in their early evolutionary phase is still far from being completely understood. Here, we present Atacama Large Millimeter/Submillimeter Array (ALMA) 1.2 mm full polarized continuum, and H$^{13}$CO$^+$(3$-$2), CS(5$-$4), and HN$^{13}$C(3$-$2) line observations with a high angular resolution ($\sim$0.4… ▽ More The formation of the massive stars, and in particular, the role that the magnetic fields play in their early evolutionary phase is still far from being completely understood. Here, we present Atacama Large Millimeter/Submillimeter Array (ALMA) 1.2 mm full polarized continuum, and H$^{13}$CO$^+$(3$-$2), CS(5$-$4), and HN$^{13}$C(3$-$2) line observations with a high angular resolution ($\sim$0.4$''$ or 1100 au). In the 1.2 mm continuum emission, we reveal a dusty envelope surrounding the massive protostars, IRAS16547-E and IRAS16547-W, with dimensions of $\sim$10,000 au. This envelope has a bi-conical structure likely carved by the powerful thermal radio jet present in region. The magnetic fields vectors follow very-well the bi-conical envelope. The polarization fraction is $\sim$2.0\% in this region. Some of these vectors seem to converge to IRAS 16547-E, and IRAS 16547-W, the most massive protostars. Moreover, the velocity fields revealed from the spectral lines H$^{13}$CO$^+$(3$-$2), and HN$^{13}$C(3$-$2) show velocity gradients with a good correspondence with the magnetic fields, that maybe are tracing the cavities of molecular outflows or maybe in some parts infall. We derived a magnetic field strength in some filamentary regions that goes from 2 to 6.1\,mG. We also find that the CS(5$-$4) molecular line emission reveals multiple outflow cavities or bow-shocks with different orientations, some of which seem to follow the NW-SE radio thermal jet. △ Less

Submitted 19 August, 2024; originally announced August 2024.

Comments: Accepted by the Astrophysical Journal, 13 pages

Constraint Energy Minimizing Generalized Multiscale Finite Element Method for Convection Diffusion Equations with Inhomogeneous Boundary Conditions

Authors: Po Chai Wong, Eric T. Chung, Changqing Ye, Lina Zhao

Abstract: In this paper, we develop the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) for convection-diffusion equations with inhomogeneous Dirichlet, Neumann and Robin boundary conditions, along with high-contrast coefficients. For time independent problems, boundary correctors $\mathcal{D}^m$ and $\mathcal{N}^{m}$ for Dirichlet, Neumann, and Robin conditions are de… ▽ More In this paper, we develop the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) for convection-diffusion equations with inhomogeneous Dirichlet, Neumann and Robin boundary conditions, along with high-contrast coefficients. For time independent problems, boundary correctors $\mathcal{D}^m$ and $\mathcal{N}^{m}$ for Dirichlet, Neumann, and Robin conditions are designed. For time dependent problems, a scheme to update the boundary correctors is formulated. Error analysis in both cases is given to show the first-order convergence in energy norm with respect to the coarse mesh size $H$ and second-order convergence in $L^2-$norm, as verified by numerical examples, with which different finite difference schemes are compared for temporal discretization. Nonlinear problems are also demonstrated in combination with Strang splitting. △ Less

Submitted 1 August, 2024; originally announced August 2024.

Comments: 36 pages.11 figures. Submitted to Journal of Computational Mathematics

MSC Class: 65M12; 65M15; 65N30

Multiscale modeling for a class of high-contrast heterogeneous sign-changing problems

Authors: Changqing Ye, Xingguang Jin, Patrick Ciarlet Jr., Eric T. Chung

Abstract: The mathematical formulation of sign-changing problems involves a linear second-order partial differential equation in the divergence form, where the coefficient can assume positive and negative values in different subdomains. These problems find their physical background in negative-index metamaterials, either as inclusions embedded into common materials as the matrix or vice versa. In this paper… ▽ More The mathematical formulation of sign-changing problems involves a linear second-order partial differential equation in the divergence form, where the coefficient can assume positive and negative values in different subdomains. These problems find their physical background in negative-index metamaterials, either as inclusions embedded into common materials as the matrix or vice versa. In this paper, we propose a numerical method based on the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) specifically designed for sign-changing problems. The construction of auxiliary spaces in the original CEM-GMsFEM is tailored to accommodate the sign-changing setting. The numerical results demonstrate the effectiveness of the proposed method in handling sophisticated coefficient profiles and the robustness of coefficient contrast ratios. Under several technical assumptions and by applying the \texttt{T}-coercivity theory, we establish the inf-sup stability and provide an a priori error estimate for the proposed method. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Magnetic Fields in Massive Star-forming Regions (MagMaR): Unveiling an Hourglass Magnetic Field in G333.46-0.16 using ALMA

Authors: Piyali Saha, Patricio Sanhueza, Marco Padovani, Josep M. Girart, Paulo Cortes, Kaho Morii, Junhao Liu, A. Sanchez-Monge, Daniele Galli, Shantanu Basu, Patrick M. Koch, Maria T. Beltran, Shanghuo Li, Henrik Beuther, Ian W. Stephens, Fumitaka Nakamura, Qizhou Zhang, Wenyu Jiao, M. Fernandez-Lopez, Jihye Hwang, Eun Jung Chung, Kate Pattle, Luis A. Zapata, Fengwei Xu, Fernando A. Olguin , et al. (11 additional authors not shown)

Abstract: The contribution of the magnetic field to the formation of high-mass stars is poorly understood. We report the high-angular resolution ($\sim0.3^{\prime\prime}$, 870 au) map of the magnetic field projected on the plane of the sky (B$_\mathrm{POS}$) towards the high-mass star forming region G333.46$-$0.16 (G333), obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.2 mm as par… ▽ More The contribution of the magnetic field to the formation of high-mass stars is poorly understood. We report the high-angular resolution ($\sim0.3^{\prime\prime}$, 870 au) map of the magnetic field projected on the plane of the sky (B$_\mathrm{POS}$) towards the high-mass star forming region G333.46$-$0.16 (G333), obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.2 mm as part of the Magnetic Fields in Massive Star-forming Regions (MagMaR) survey. The B$_\mathrm{POS}$ morphology found in this region is consistent with a canonical ``hourglass'' which suggest a dynamically important field. This region is fragmented into two protostars separated by $\sim1740$ au. Interestingly, by analysing H$^{13}$CO$^{+}$ ($J=3-2$) line emission, we find no velocity gradient over the extend of the continuum which is consistent with a strong field. We model the B$_\mathrm{POS}$, obtaining a marginally supercritical mass-to-flux ratio of 1.43, suggesting an initially strongly magnetized environment. Based on the Davis-Chandrasekhar-Fermi method, the magnetic field strength towards G333 is estimated to be 5.7 mG. The absence of strong rotation and outflows towards the central region of G333 suggests strong magnetic braking, consistent with a highly magnetized environment. Our study shows that despite being a strong regulator, the magnetic energy fails to prevent the process of fragmentation, as revealed by the formation of the two protostars in the central region. △ Less

Submitted 23 July, 2024; originally announced July 2024.

Robust Multiscale Methods for Helmholtz equations in high contrast heterogeneous media

Authors: Xingguang Jin, Changqing Ye, Eric T. Chung

Abstract: In this paper, we provide the constraint energy minimization generalized multiscale finite element method (CEM-GMsFEM) to solve Helmholtz equations in heterogeneous medium. This novel multiscale method is specifically designed to overcome problems related to pollution effect, high-contrast coefficients, and the loss of hermiticity of operators. We establish the inf-sup stability and give an a prio… ▽ More In this paper, we provide the constraint energy minimization generalized multiscale finite element method (CEM-GMsFEM) to solve Helmholtz equations in heterogeneous medium. This novel multiscale method is specifically designed to overcome problems related to pollution effect, high-contrast coefficients, and the loss of hermiticity of operators. We establish the inf-sup stability and give an a priori error estimate for this method under a number of established assumptions and resolution conditions. The theoretical results are validated by a set of numerical tests, which further show that the multiscale technique can effectively capture pertinent physical phenomena. △ Less

Submitted 8 July, 2024; v1 submitted 5 July, 2024; originally announced July 2024.

Thinking Inside The Box: Privacy Against Stronger Adversaries

Authors: Eldon Chung

Abstract: In this thesis, we study extensions of statistical cryptographic primitives. In particular we study leakage-resilient secret sharing, non-malleable extractors, and immunized ideal one-way functions. The thesis is divided into three main chapters. In the first chapter, we show that 2-out-of-2 leakage resilient (and also non-malleable) secret sharing requires randomness sources that are also extract… ▽ More In this thesis, we study extensions of statistical cryptographic primitives. In particular we study leakage-resilient secret sharing, non-malleable extractors, and immunized ideal one-way functions. The thesis is divided into three main chapters. In the first chapter, we show that 2-out-of-2 leakage resilient (and also non-malleable) secret sharing requires randomness sources that are also extractable. This rules out the possibility of using min-entropic sources. In the second, we introduce collision-resistant seeded extractors and show that any seeded extractor can be made collision resistant at a small overhead in seed length. We then use it to give a two-source non-malleable extractor with entropy rate 0.81 in one source and polylogarithmic in the other. The non-malleable extractor lead to the first statistical privacy amplification protocol against memory tampering adversaries. In the final chapter, we study the hardness of the data structure variant of the 3SUM problem which is motivated by a recent construction to immunise random oracles against pre-processing adversaries. We give worst-case data structure hardness for the 3SUM problem matching known barriers in data structures for adaptive adversaries. We also give a slightly stronger lower bound in the case of non-adaptivity. Lastly, we give a novel result in the bit-probe setting. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: PhD thesis

MagMar III -- Resisting the Pressure, Is the Magnetic Field Overwhelmed in NGC6334I?

Authors: Paulo C. Cortes, Josep M. Girart, Patricio Sanhueza, Junhao Liu, Sergio Martin, Ian W. Stephens, Henrik Beuther, Patrick M. Koch, M. Fernandez-Lopez, Alvaro Sanchez-Monge, Jia-Wei Wang, Kaho Morii, Shanghuo Li, Piyali Saha, Qizhou Zhang, David Rebolledo, Luis A. Zapata, Ji-hyun Kang, Wenyu Jiao, Jongsoo Kim, Yu Cheng, Jihye Hwang, Eun Jung Chung, Spandan Choudhury, A-Ran Lyo , et al. (1 additional authors not shown)

Abstract: We report on ALMA observations of polarized dust emission at 1.2 mm from NGC6334I, a source known for its significant flux outbursts. Between five months, our data show no substantial change in total intensity and a modest 8\% variation in linear polarization, suggesting a phase of stability or the conclusion of the outburst. The magnetic field, inferred from this polarized emission, displays a pr… ▽ More We report on ALMA observations of polarized dust emission at 1.2 mm from NGC6334I, a source known for its significant flux outbursts. Between five months, our data show no substantial change in total intensity and a modest 8\% variation in linear polarization, suggesting a phase of stability or the conclusion of the outburst. The magnetic field, inferred from this polarized emission, displays a predominantly radial pattern from North-West to South-East with intricate disturbances across major cores, hinting at spiral structures. Energy analysis of CS$(J=5 \rightarrow 4)$ emission yields an outflow energy of approximately $3.5\times10^{45}$ ergs, aligning with previous interferometric studies. Utilizing the Davis-Chandrasekhar-Fermi method, we determined magnetic field strengths ranging from 1 to 11 mG, averaging at 1.9 mG. This average increases to 4 $\pm 1$ mG when incorporating Zeeman measurements. Comparative analyses using gravitational, thermal, and kinetic energy maps reveal that magnetic energy is significantly weaker, possibly explaining the observed field morphology. We also find that the energy in the outflows and the expanding cometary {\HII} region is also larger than the magnetic energy, suggesting that protostellar feedback maybe the dominant driver behind the injection of turbulence in NGC6334I at the scales sampled by our data. The gas in NGC6334I predominantly exhibits supersonic and trans-Alfvenic conditions, transitioning towards a super-Alfvenic regime, underscoring a diminished influence of the magnetic field with increasing gas density. These observations are in agreement with prior polarization studies at 220 GHz, enriching our understanding of the dynamic processes in high-mass star-forming regions. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: Accepted for Publication at the Astrophysical Journal

An iterative constraint energy minimizing generalized multiscale finite element method for contact problem

Authors: Zishang Li, Changqing Ye, Eric T. Chung

Abstract: This work presents an Iterative Constraint Energy Minimizing Generalized Multiscale Finite Element Method (ICEM-GMsFEM) for solving the contact problem with high contrast coefficients. The model problem can be characterized by a variational inequality, where we add a penalty term to convert this problem into a non-smooth and non-linear unconstrained minimizing problem. The characterization of the… ▽ More This work presents an Iterative Constraint Energy Minimizing Generalized Multiscale Finite Element Method (ICEM-GMsFEM) for solving the contact problem with high contrast coefficients. The model problem can be characterized by a variational inequality, where we add a penalty term to convert this problem into a non-smooth and non-linear unconstrained minimizing problem. The characterization of the minimizer satisfies the variational form of a mixed Dirilect-Neumann-Robin boundary value problem. So we apply CEM-GMsFEM iteratively and introduce special boundary correctors along with multiscale spaces to achieve an optimal convergence rate. Numerical results are conducted for different highly heterogeneous permeability fields, validating the fast convergence of the CEM-GMsFEM iteration in handling the contact boundary and illustrating the stability of the proposed method with different sets of parameters. We also prove the fast convergence of the proposed iterative CEM-GMsFEM method and provide an error estimate of the multiscale solution under a mild assumption. △ Less

Submitted 5 June, 2024; originally announced June 2024.

The Formation of Filaments and Dense cores in the Cocoon Nebula (IC~5146)

Authors: Eun Jung Chung, Chang Won Lee, Shinyoung Kim, Mario Tafalla, Hyunju Yoo, Jungyeon Cho, Woojin Kwon

Abstract: We present 850~$μ$m linear polarization and C$^{18}$O~(3-2) and $^{13}$CO~(3-2) molecular line observations toward the filaments (F13 and F13S) in the Cocoon Nebula (IC~5146) using the JCMT POL-2 and HARP instruments. F13 and F13S are found to be thermally supercritical with identified dense cores along their crests. Our findings include that the polarization fraction decreases in denser regions,… ▽ More We present 850~$μ$m linear polarization and C$^{18}$O~(3-2) and $^{13}$CO~(3-2) molecular line observations toward the filaments (F13 and F13S) in the Cocoon Nebula (IC~5146) using the JCMT POL-2 and HARP instruments. F13 and F13S are found to be thermally supercritical with identified dense cores along their crests. Our findings include that the polarization fraction decreases in denser regions, indicating reduced dust grain alignment efficiency. The magnetic field vectors at core scales tend to be parallel to the filaments, but disturbed at the high density regions. Magnetic field strengths measured using the Davis-Chandrasekhar-Fermi method are 58$\pm$31 and 40$\pm$9~$μ$G for F13 and F13S, respectively, and it reveals subcritical and sub-Alfvénic filaments, emphasizing the importance of magnetic fields in the Cocoon region. Sinusoidal C$^{18}$O~(3-2) velocity and density distributions are observed along the filaments' skeletons, and their variations are mostly displaced by $\sim1/4 \times$wavelength of the sinusoid, indicating core formation occurred through the fragmentation of a gravitationally unstable filament, but with shorter core spacings than predicted. Large scale velocity fields of F13 and F13S, studied using $^{13}$CO~(3-2) data, present V-shape transverse velocity structure. We propose a scenario for the formation and evolution of F13 and F13S, along with the dense cores within them. A radiation shock front generated by a B-type star collided with a sheet-like cloud about 1.4~Myr ago, the filaments became thermally critical due to mass infall through self-gravity $\sim$1~Myr ago, and subsequently dense cores formed through gravitational fragmentation, accompanied by the disturbance of the magnetic field. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: 26 pages, 19 figures, 3 tables, Accepted for publication in ApJ

CEM-GMsFEM for Poisson equations in heterogeneous perforated domains

Authors: Wei Xie, Yin Yang, Eric Chung, Yunqing Huang

Abstract: In this paper, we propose a novel multiscale model reduction strategy tailored to address the Poisson equation within heterogeneous perforated domains. The numerical simulation of this intricate problem is impeded by its multiscale characteristics, necessitating an exceptionally fine mesh to adequately capture all relevant details. To overcome the challenges inherent in the multiscale nature of th… ▽ More In this paper, we propose a novel multiscale model reduction strategy tailored to address the Poisson equation within heterogeneous perforated domains. The numerical simulation of this intricate problem is impeded by its multiscale characteristics, necessitating an exceptionally fine mesh to adequately capture all relevant details. To overcome the challenges inherent in the multiscale nature of the perforations, we introduce a coarse space constructed using the Constraint Energy Minimizing Generalized Multiscale Finite Element Method (CEM-GMsFEM). This involves constructing basis functions through a sequence of local energy minimization problems over eigenspaces containing localized information pertaining to the heterogeneities. Through our analysis, we demonstrate that the oversampling layers depend on the local eigenvalues, thereby implicating the local geometry as well. Additionally, we provide numerical examples to illustrate the efficacy of the proposed scheme. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Improving Retrieval for RAG based Question Answering Models on Financial Documents

Authors: Spurthi Setty, Harsh Thakkar, Alyssa Lee, Eden Chung, Natan Vidra

Abstract: The effectiveness of Large Language Models (LLMs) in generating accurate responses relies heavily on the quality of input provided, particularly when employing Retrieval Augmented Generation (RAG) techniques. RAG enhances LLMs by sourcing the most relevant text chunk(s) to base queries upon. Despite the significant advancements in LLMs' response quality in recent years, users may still encounter i… ▽ More The effectiveness of Large Language Models (LLMs) in generating accurate responses relies heavily on the quality of input provided, particularly when employing Retrieval Augmented Generation (RAG) techniques. RAG enhances LLMs by sourcing the most relevant text chunk(s) to base queries upon. Despite the significant advancements in LLMs' response quality in recent years, users may still encounter inaccuracies or irrelevant answers; these issues often stem from suboptimal text chunk retrieval by RAG rather than the inherent capabilities of LLMs. To augment the efficacy of LLMs, it is crucial to refine the RAG process. This paper explores the existing constraints of RAG pipelines and introduces methodologies for enhancing text retrieval. It delves into strategies such as sophisticated chunking techniques, query expansion, the incorporation of metadata annotations, the application of re-ranking algorithms, and the fine-tuning of embedding algorithms. Implementing these approaches can substantially improve the retrieval quality, thereby elevating the overall performance and reliability of LLMs in processing and responding to queries. △ Less

Submitted 31 July, 2024; v1 submitted 22 March, 2024; originally announced April 2024.

A fast cosine transformation accelerated method for predicting effective thermal conductivity

Authors: Changqing Ye, Shubin Fu, Eric T. Chung

Abstract: Predicting effective thermal conductivity by solving a Partial Differential Equation (PDE) defined on a high-resolution Representative Volume Element (RVE) is a computationally intensive task. In this paper, we tackle the task by proposing an efficient and implementation-friendly computational method that can fully leverage the computing power offered by hardware accelerators, namely, graphical pr… ▽ More Predicting effective thermal conductivity by solving a Partial Differential Equation (PDE) defined on a high-resolution Representative Volume Element (RVE) is a computationally intensive task. In this paper, we tackle the task by proposing an efficient and implementation-friendly computational method that can fully leverage the computing power offered by hardware accelerators, namely, graphical processing units (GPUs). We first employ the Two-Point Flux-Approximation scheme to discretize the PDE and then utilize the preconditioned conjugate gradient method to solve the resulting algebraic linear system. The construction of the preconditioner originates from FFT-based homogenization methods, and an engineered linear programming technique is utilized to determine the homogeneous reference parameters. The fundamental observation presented in this paper is that the preconditioner system can be effectively solved using multiple Fast Cosine Transformations (FCT) and parallel tridiagonal matrix solvers. Regarding the fact that default multiple FCTs are unavailable on the CUDA platform, we detail how to derive FCTs from FFTs with nearly optimal memory usage. Numerical experiments including the stability comparison with standard preconditioners are conducted for 3D RVEs. Our performance reports indicate that the proposed method can achieve a $5$-fold acceleration on the GPU platform over the pure CPU platform and solve the problems with $512^3$ degrees of freedom and reasonable contrast ratios in less than $30$ seconds. △ Less

Submitted 2 April, 2024; originally announced April 2024.

A robust two-level overlapping preconditioner for Darcy flow in high-contrast media

Authors: Changqing Ye, Shubin Fu, Eric T. Chung, Jizu Huang

Abstract: In this article, a two-level overlapping domain decomposition preconditioner is developed for solving linear algebraic systems obtained from simulating Darcy flow in high-contrast media. Our preconditioner starts at a mixed finite element method for discretizing the partial differential equation by Darcy's law with the no-flux boundary condition and is then followed by a velocity elimination techn… ▽ More In this article, a two-level overlapping domain decomposition preconditioner is developed for solving linear algebraic systems obtained from simulating Darcy flow in high-contrast media. Our preconditioner starts at a mixed finite element method for discretizing the partial differential equation by Darcy's law with the no-flux boundary condition and is then followed by a velocity elimination technique to yield a linear algebraic system with only unknowns of pressure. Then, our main objective is to design a robust and efficient domain decomposition preconditioner for this system, which is accomplished by engineering a multiscale coarse space that is capable of characterizing high-contrast features of the permeability field. A generalized eigenvalue problem is solved in each non-overlapping coarse element in a communication-free manner to form the global solver, which is accompanied by local solvers originated from additive Schwarz methods but with a non-Galerkin discretization to derive the two-level preconditioner. We provide a rigorous analysis that indicates that the condition number of the preconditioned system could be bounded above with several assumptions. Extensive numerical experiments with various types of three-dimensional high-contrast models are exhibited. In particular, we study the robustness against the contrast of the media as well as the influences of numbers of eigenfunctions, oversampling sizes, and subdomain partitions on the efficiency of the proposed preconditioner. Besides, strong and weak scalability performances are also examined. △ Less

Submitted 28 March, 2024; originally announced March 2024.

An efficient multiscale multigrid preconditioner for Darcy flow in high-contrast media

Authors: Changqing Ye, Shubin Fu, Eric T. Chung, Jizu Huang

Abstract: In this paper, we develop a multigrid preconditioner to solve Darcy flow in highly heterogeneous porous media. The key component of the preconditioner is to construct a sequence of nested subspaces $W_{\mathcal{L}}\subset W_{\mathcal{L}-1}\subset\cdots\subset W_1=W_h$. An appropriate spectral problem is defined in the space of $W_{i-1}$, then the eigenfunctions of the spectral problems are utilize… ▽ More In this paper, we develop a multigrid preconditioner to solve Darcy flow in highly heterogeneous porous media. The key component of the preconditioner is to construct a sequence of nested subspaces $W_{\mathcal{L}}\subset W_{\mathcal{L}-1}\subset\cdots\subset W_1=W_h$. An appropriate spectral problem is defined in the space of $W_{i-1}$, then the eigenfunctions of the spectral problems are utilized to form $W_i$. The preconditioner is applied to solve a positive semidefinite linear system which results from discretizing the Darcy flow equation with the lowest order Raviart-Thomas spaces and adopting a trapezoidal quadrature rule. Theoretical analysis and numerical investigations of this preconditioner will be presented. In particular, we will consider several typical highly heterogeneous permeability fields whose resolutions are up to $1024^3$ and examine the computational performance of the preconditioner in several aspects, such as strong scalability, weak scalability, and robustness against the contrast of the media. We also demonstrate an application of this preconditioner for solving a two-phase flow benchmark problem. △ Less

Submitted 28 March, 2024; originally announced March 2024.

Enhancing Large Language Model Performance To Answer Questions and Extract Information More Accurately

Authors: Liang Zhang, Katherine Jijo, Spurthi Setty, Eden Chung, Fatima Javid, Natan Vidra, Tommy Clifford

Abstract: Large Language Models (LLMs) generate responses to questions; however, their effectiveness is often hindered by sub-optimal quality of answers and occasional failures to provide accurate responses to questions. To address these challenges, a fine-tuning process is employed, involving feedback and examples to refine models. The objective is to enhance AI models through continuous feedback loops, ut… ▽ More Large Language Models (LLMs) generate responses to questions; however, their effectiveness is often hindered by sub-optimal quality of answers and occasional failures to provide accurate responses to questions. To address these challenges, a fine-tuning process is employed, involving feedback and examples to refine models. The objective is to enhance AI models through continuous feedback loops, utilizing metrics such as cosine similarity, LLM evaluation and Rouge-L scores to evaluate the models. Leveraging LLMs like GPT-3.5, GPT4ALL, and LLaMA2, and Claude, this approach is benchmarked on financial datasets, including the FinanceBench and RAG Instruct Benchmark Tester Dataset, illustrating the necessity of fine-tuning. The results showcase the capability of fine-tuned models to surpass the accuracy of zero-shot LLMs, providing superior question and answering capabilities. Notably, the combination of fine-tuning the LLM with a process known as Retrieval Augmented Generation (RAG) proves to generate responses with improved accuracy. △ Less

Submitted 26 January, 2024; originally announced February 2024.

Filamentary Network and Magnetic Field Structures Revealed with BISTRO in the High-Mass Star-Forming Region NGC2264 : Global Properties and Local Magnetogravitational Configurations

Authors: Jia-Wei Wang, Patrick M. Koch, Seamus D. Clarke, Gary Fuller, Nicolas Peretto, Ya-Wen Tang, Hsi-Wei Yen, Shih-Ping Lai, Nagayoshi Ohashi, Doris Arzoumanian, Doug Johnstone, Ray Furuya, Shu-ichiro Inutsuka, Chang Won Lee, Derek Ward-Thompson, Valentin J. M. Le Gouellec, Hong-Li Liu, Lapo Fanciullo, Jihye Hwang, Kate Pattle, Frédérick Poidevin, Mehrnoosh Tahani, Takashi Onaka, Mark G. Rawlings, Eun Jung Chung , et al. (132 additional authors not shown)

Abstract: We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from… ▽ More We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astrophysical Journal. 43 pages, 32 figures, and 4 tables (including Appendix)

Improving Classification Performance With Human Feedback: Label a few, we label the rest

Authors: Natan Vidra, Thomas Clifford, Katherine Jijo, Eden Chung, Liang Zhang

Abstract: In the realm of artificial intelligence, where a vast majority of data is unstructured, obtaining substantial amounts of labeled data to train supervised machine learning models poses a significant challenge. To address this, we delve into few-shot and active learning, where are goal is to improve AI models with human feedback on a few labeled examples. This paper focuses on understanding how a co… ▽ More In the realm of artificial intelligence, where a vast majority of data is unstructured, obtaining substantial amounts of labeled data to train supervised machine learning models poses a significant challenge. To address this, we delve into few-shot and active learning, where are goal is to improve AI models with human feedback on a few labeled examples. This paper focuses on understanding how a continuous feedback loop can refine models, thereby enhancing their accuracy, recall, and precision through incremental human input. By employing Large Language Models (LLMs) such as GPT-3.5, BERT, and SetFit, we aim to analyze the efficacy of using a limited number of labeled examples to substantially improve model accuracy. We benchmark this approach on the Financial Phrasebank, Banking, Craigslist, Trec, Amazon Reviews datasets to prove that with just a few labeled examples, we are able to surpass the accuracy of zero shot large language models to provide enhanced text classification performance. We demonstrate that rather than needing to manually label millions of rows of data, we just need to label a few and the model can effectively predict the rest. △ Less

Submitted 17 January, 2024; originally announced January 2024.

Multiple Beads-on-a-string: Dark Matter-Deficient Galaxy Formation in a Mini-bullet Satellite-satellite Galaxy Collision

Authors: Joohyun Lee, Eun-jin Shin, Ji-hoon Kim, Paul R. Shapiro, Eunwoo Chung

Abstract: Dark matter-deficient galaxies (DMDGs) discovered in the survey of ultra-diffuse galaxies (UDGs), in apparent conflict with standard CDM, may be produced by high-velocity galaxy-galaxy collisions, the $\textit{Mini-bullet}$ scenario. Recent observations of an aligned trail of $7-11$ UDGs near NGC1052, including DMDGs DF2 and DF4, suggesting a common formation event, $\sim8.9\pm1.5$ Gyr ago, provid… ▽ More Dark matter-deficient galaxies (DMDGs) discovered in the survey of ultra-diffuse galaxies (UDGs), in apparent conflict with standard CDM, may be produced by high-velocity galaxy-galaxy collisions, the $\textit{Mini-bullet}$ scenario. Recent observations of an aligned trail of $7-11$ UDGs near NGC1052, including DMDGs DF2 and DF4, suggesting a common formation event, $\sim8.9\pm1.5$ Gyr ago, provide a test. Hydro/N-body simulations, supplemented by galaxy orbit integrations, demonstrate that satellite-satellite collisions outside the host-galaxy virial radius can reproduce the observed UDGs in the NGC1052 group. A trail of $\sim10$ DMDGs is shown to form, including two massive ones that replicate the observed motions of DF2 and DF4. The linear relation, $v=Ax+v_{0}$, conjectured previously to relate positions ($x$) and velocities ($v$) of the aligned DMDGs as a signature of the collision event, is approximately obeyed, but individual DMDGs can deviate significantly from it. The progenitors whose collision spawned the trail of DMDGs survive the collision without, themselves, becoming DMDGs. We predict one progenitor is located at the end of the trail, testable by observing the difference between its stars, formed pre-collision, from those of the DMDGs, formed post-collision. By contrast, stellar ages and metallicities of the DMDGs are nearly identical. We further offer a hint that the tidal field of host NGC1052 may contribute to making DMDGs diffuse. $Λ$CDM simulation in a 100 cMpc box finds our required initial conditions $\sim10$ times at $z<3$. These results indicate current observations are consistent with the $\textit{Mini-bullet}$ scenario. △ Less

Submitted 9 February, 2024; v1 submitted 18 December, 2023; originally announced December 2023.

Comments: 27 pages, 14 figures, minor revisions for clarification only in response to referee report, + 2 new figures. For associated movies and information, see https://joohyun-lee.github.io/research/DMDG3/

Report number: UT-WI-11-2024

Journal ref: ApJ, 966, 72 (2024)

A second-order exponential integration constraint energy minimizing generalized multiscale method for parabolic problems

Authors: Leonardo A. Poveda, Juan Galvis, Eric Chung

Abstract: This paper investigates an efficient exponential integrator generalized multiscale finite element method for solving a class of time-evolving partial differential equations in bounded domains. The proposed method first performs the spatial discretization of the model problem using constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). This approach consists of two… ▽ More This paper investigates an efficient exponential integrator generalized multiscale finite element method for solving a class of time-evolving partial differential equations in bounded domains. The proposed method first performs the spatial discretization of the model problem using constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). This approach consists of two stages. First, the auxiliary space is constructed by solving local spectral problems, where the basis functions corresponding to small eigenvalues are captured. The multiscale basis functions are obtained in the second stage using the auxiliary space by solving local energy minimization problems over the oversampling domains. The basis functions have exponential decay outside the corresponding local oversampling regions. We shall consider the first and second-order explicit exponential Runge-Kutta approach for temporal discretization and to build a fully discrete numerical solution. The exponential integration strategy for the time variable allows us to take full advantage of the CEM-GMsFEM as it enables larger time steps due to its stability properties. We derive the error estimates in the energy norm under the regularity assumption. Finally, we will provide some numerical experiments to sustain the efficiency of the proposed method. △ Less

Submitted 5 July, 2024; v1 submitted 17 October, 2023; originally announced October 2023.

MSC Class: 65M15; 65M60; 65M12; 65M22

Microscaling Data Formats for Deep Learning

Authors: Bita Darvish Rouhani, Ritchie Zhao, Ankit More, Mathew Hall, Alireza Khodamoradi, Summer Deng, Dhruv Choudhary, Marius Cornea, Eric Dellinger, Kristof Denolf, Stosic Dusan, Venmugil Elango, Maximilian Golub, Alexander Heinecke, Phil James-Roxby, Dharmesh Jani, Gaurav Kolhe, Martin Langhammer, Ada Li, Levi Melnick, Maral Mesmakhosroshahi, Andres Rodriguez, Michael Schulte, Rasoul Shafipour, Lei Shao , et al. (8 additional authors not shown)

Abstract: Narrow bit-width data formats are key to reducing the computational and storage costs of modern deep learning applications. This paper evaluates Microscaling (MX) data formats that combine a per-block scaling factor with narrow floating-point and integer types for individual elements. MX formats balance the competing needs of hardware efficiency, model accuracy, and user friction. Empirical result… ▽ More Narrow bit-width data formats are key to reducing the computational and storage costs of modern deep learning applications. This paper evaluates Microscaling (MX) data formats that combine a per-block scaling factor with narrow floating-point and integer types for individual elements. MX formats balance the competing needs of hardware efficiency, model accuracy, and user friction. Empirical results on over two dozen benchmarks demonstrate practicality of MX data formats as a drop-in replacement for baseline FP32 for AI inference and training with low user friction. We also show the first instance of training generative language models at sub-8-bit weights, activations, and gradients with minimal accuracy loss and no modifications to the training recipe. △ Less

Submitted 19 October, 2023; v1 submitted 16 October, 2023; originally announced October 2023.

Wavelet-based Edge Multiscale Finite Element Methods for Singularly Perturbed Convection-Diffusion Equations

Authors: Shubin Fu, Eric Chung, Guanglian Li

Abstract: We propose a novel efficient and robust Wavelet-based Edge Multiscale Finite Element Method (WEMsFEM) motivated by \cite{MR3980476,GL18} to solve the singularly perturbed convection-diffusion equations. The main idea is to first establish a local splitting of the solution over a local region by a local bubble part and local Harmonic extension part, and then derive a global splitting by means of Pa… ▽ More We propose a novel efficient and robust Wavelet-based Edge Multiscale Finite Element Method (WEMsFEM) motivated by \cite{MR3980476,GL18} to solve the singularly perturbed convection-diffusion equations. The main idea is to first establish a local splitting of the solution over a local region by a local bubble part and local Harmonic extension part, and then derive a global splitting by means of Partition of Unity. This facilitates a representation of the solution as a summation of a global bubble part and a global Harmonic extension part, where the first part can be computed locally in parallel. To approximate the second part, we construct an edge multiscale ansatz space locally with hierarchical bases as the local boundary data that has a guaranteed approximation rate \noteLg{both inside and outside of the layers}. The key innovation of this proposed WEMsFEM lies in a provable convergence rate with little restriction on the mesh size. Its convergence rate with respect to the computational degree of freedom is rigorously analyzed, which is verified by extensive 2-d and 3-d numerical tests. △ Less

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

TRAO Survey of Nearby Filamentary Molecular clouds, the Universal Nursery of Stars (TRAO-FUNS). III. Filaments and dense cores in the NGC 2068 and NGC 2071 regions of Orion B

Authors: Hyunju Yoo, Chang Won Lee, Eun Jung Chung, Shinyoung Kim, Mario Tafalla, Paola Caselli, Philip C. Myers, Kyoung Hee Kim, Tie Liu, Woojin Kwon, Archana Soam, Jongsoo Kim

Abstract: We present the results of molecular line observations performed toward the NGC 2068 and NGC 2071 regions of the Orion B cloud as the TRAO-FUNS project to study the roles of the filamentary structure in the formation of dense cores and stars in the clouds. Gaussian decomposition for the C$^{18}$O spectra with multiple velocity components and application of a Friends-of-Friends algorithm for the dec… ▽ More We present the results of molecular line observations performed toward the NGC 2068 and NGC 2071 regions of the Orion B cloud as the TRAO-FUNS project to study the roles of the filamentary structure in the formation of dense cores and stars in the clouds. Gaussian decomposition for the C$^{18}$O spectra with multiple velocity components and application of a Friends-of-Friends algorithm for the decomposed components allowed us to identify a few tens of velocity coherent filaments. We also identified 48 dense cores from the observations of N$_2$H$^{+}$ using a core finding tool, FellWalker. We made the virial analysis for these filaments and dense cores, finding that the filaments with N$_2$H$^{+}$ dense core are thermally supercritical, and the filaments with larger ratio between the line mass and the thermal critical line mass tend to have more dense cores. We investigated the contribution of the nonthermal motions in dense cores and filaments, showing the dense cores are mostly in transonic/subsonic motions while their natal filaments are mostly in supersonic motions. This may indicates that gas turbulent motions in the filaments have been dissipated at the core scale to form the dense cores there. The filaments with (dynamically evolved) dense cores in infalling motions or with NH$_2$D bright (or chemically evolved) dense cores are all found to be gravitationally critical. Therefore, the criticality of the filament is thought to provide a key condition for its fragmentation, the formation of dense cores, and their kinematical and chemical evolution. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 43 pages, 34 figures (more figures are available in the online journal), 4 tables, Accepted for publication in Astrophysical Journal

Multicontinuum homogenization. General theory and applications

Authors: E. Chung, Y. Efendiev, J. Galvis, W. T. Leung

Abstract: In this paper, we discuss a general framework for multicontinuum homogenization. Multicontinuum models are widely used in many applications and some derivations for these models are established. In these models, several macroscopic variables at each macroscale point are defined and the resulting multicontinuum equations are formulated. In this paper, we propose a general formulation and associated… ▽ More In this paper, we discuss a general framework for multicontinuum homogenization. Multicontinuum models are widely used in many applications and some derivations for these models are established. In these models, several macroscopic variables at each macroscale point are defined and the resulting multicontinuum equations are formulated. In this paper, we propose a general formulation and associated ingredients that allow performing multicontinuum homogenization. Our derivation consists of several main parts. In the first part, we propose a general expansion, where the solution is expressed via the product of multiple macro variables and associated cell problems. The second part consists of formulating the cell problems. The cell problems are formulated as saddle point problems with constraints for each continua. Defining the continua via test functions, we set the constraints as an integral representation. Finally, substituting the expansion to the original system, we obtain multicontinuum systems. We present an application to the mixed formulation of elliptic equations. This is a challenging system as the system does not have symmetry. We discuss the local problems and various macroscale representations for the solution and its gradient. Using various order approximations, one can obtain different systems of equations. We discuss the applicability of multicontinuum homogenization and relate this to high contrast in the cell problem. Numerical results are presented. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Photometry of Type II Supernova SN 2023ixf with a Worldwide Citizen Science Network

Authors: Lauren A. Sgro, Thomas M. Esposito, Guillaume Blaclard, Sebastian Gomez, Franck Marchis, Alexei V. Filippenko, Daniel O'Conner Peluso, Stephen S. Lawrence, Aad Verveen, Andreas Wagner, Anouchka Nardi, Barbara Wiart, Benjamin Mirwald, Bill Christensen, Bob Eramia, Bruce Parker, Bruno Guillet, Byungki Kim, Chelsey A. Logan, Christopher C. M. Kyba, Christopher Toulmin, Claudio G. Vantaggiato, Dana Adhis, Dave Gary, Dave Goodey , et al. (66 additional authors not shown)

Abstract: We present highly sampled photometry of the supernova (SN) 2023ixf, a Type II SN in M101, beginning 2 days before its first known detection. To gather these data, we enlisted the global Unistellar Network of citizen scientists. These 252 observations from 115 telescopes show the SN's rising brightness associated with shock emergence followed by gradual decay. We measure a peak $M_{V}$ = -18.18… ▽ More We present highly sampled photometry of the supernova (SN) 2023ixf, a Type II SN in M101, beginning 2 days before its first known detection. To gather these data, we enlisted the global Unistellar Network of citizen scientists. These 252 observations from 115 telescopes show the SN's rising brightness associated with shock emergence followed by gradual decay. We measure a peak $M_{V}$ = -18.18 $\pm$ 0.09 mag at 2023-05-25 21:37 UTC in agreement with previously published analyses. △ Less

Submitted 7 July, 2023; originally announced July 2023.

Comments: 4 pages, 1 figure

Journal ref: Res. Notes AAS 7 141 (2023)

The JCMT BISTRO Survey: Studying the Complex Magnetic Field of L43

Authors: Janik Karoly, Derek Ward-Thompson, Kate Pattle, David Berry, Anthony Whitworth, Jason Kirk, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar… ▽ More We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $μ$G in the main starless core and up to $\sim$90$\pm$40 $μ$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfvénic. We also present a new method of data reduction for these denser but fainter objects like starless cores. △ Less

Submitted 22 May, 2023; v1 submitted 18 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ. 23 pages, 9 figures (7 main text, 2 appendix)

Magnetic Fields and Fragmentation of Filaments in the Hub of California-X

Authors: Eun Jung Chung, Chang Won Lee, Woojin Kwon, Mario Tafalla, Shinyoung Kim, Archana Soam, Jungyeon Cho

Abstract: We present 850 $μ$m polarization and $\rm C^{18}O (3-2)$ molecular line observations toward the X-shaped nebula in the California molecular cloud using the JCMT SCUBA-2/POL-2 and HARP instruments. The 850 $μ$m emission shows that the observed region includes two elongated filamentary structures (Fil1 and Fil2) having chains of regularly spaced cores. We measured the mass per unit length of the fil… ▽ More We present 850 $μ$m polarization and $\rm C^{18}O (3-2)$ molecular line observations toward the X-shaped nebula in the California molecular cloud using the JCMT SCUBA-2/POL-2 and HARP instruments. The 850 $μ$m emission shows that the observed region includes two elongated filamentary structures (Fil1 and Fil2) having chains of regularly spaced cores. We measured the mass per unit length of the filament and found that Fil1 and Fil2 are thermally super- and subcritical, respectively, but both are subcritical if nonthermal turbulence is considered. The mean projected spacings ($Δ\bar S$) of cores in Fil1 and Fil2 are 0.13 and 0.16 pc, respectively. $Δ\bar S$ are smaller than $4\times$filament width expected in the classical cylinder fragmentation model. The large-scale magnetic field orientations shown by Planck are perpendicular to the long axes of Fil1 and Fil2, while those in the filaments obtained from the high-resolution polarization data of JCMT are disturbed, but those in Fil1 tend to have longitudinal orientations. Using the modified Davis-Chandrasekhar-Fermi (DCF) method, we estimated the magnetic field strengths ($B_{\rm pos}$) of filaments which are 110$\pm$80 and 90$\pm$60 $μ$G. We calculated the gravitational, kinematic, and magnetic energies of the filaments, and found that the fraction of magnetic energy is larger than 60 % in both filaments. We propose that a dominant magnetic energy may lead the filament to be fragmented into aligned cores as suggested by Tang et al., and a shorter core spacing can be due to a projection effect via the inclined geometry of filaments or due to a non-negligible, longitudinal magnetic fields in case of Fil1. △ Less

Submitted 17 May, 2023; originally announced May 2023.

Comments: Accepted for publication in the ApJ. 18 pages, 11 figures, 2 tables

Separating Key Agreement and Computational Differential Privacy

Authors: Vipul Arora, Eldon Chung, Zeyong Li, Thomas Tan

Abstract: Two party differential privacy allows two parties who do not trust each other, to come together and perform a joint analysis on their data whilst maintaining individual-level privacy. We show that any efficient, computationally differentially private protocol that has black-box access to key agreement (and nothing stronger), is also an efficient, information-theoretically differentially private pr… ▽ More Two party differential privacy allows two parties who do not trust each other, to come together and perform a joint analysis on their data whilst maintaining individual-level privacy. We show that any efficient, computationally differentially private protocol that has black-box access to key agreement (and nothing stronger), is also an efficient, information-theoretically differentially private protocol. In other words, the existence of efficient key agreement protocols is insufficient for efficient, computationally differentially private protocols. In doing so, we make progress in answering an open question posed by Vadhan about the minimal computational assumption needed for computational differential privacy. Combined with the information-theoretic lower bound due to McGregor, Mironov, Pitassi, Reingold, Talwar, and Vadhan in [FOCS'10], we show that there is no fully black-box reduction from efficient, computationally differentially private protocols for computing the Hamming distance (or equivalently inner product over the integers) on $n$ bits, with additive error lower than $O\left(\frac{\sqrt{n}}{e^ε\log(n)}\right)$, to key agreement. This complements the result by Haitner, Mazor, Silbak, and Tsfadia in [STOC'22], which showed that computing the Hamming distance implies key agreement. We conclude that key agreement is \emph{strictly} weaker than computational differential privacy for computing the inner product, thereby answering their open question on whether key agreement is sufficient. △ Less

Submitted 28 August, 2023; v1 submitted 14 April, 2023; originally announced April 2023.

Comments: A key step in relating the probability that can be computed by the PSPACE algorithm to the statistical distinguishing probability is missing and not yet shown. Our arguments in this work so far have not yet been able to show this step. Thus the final conclusion that key agreement is black-box insufficient for CDP is not yet proven

Convergence of the CEM-GMsFEM for compressible flow in highly heterogeneous media

Authors: Leonardo A. Poveda, Shubin Fu, Eric T. Chung, Lina Zhao

Abstract: This paper presents and analyses a Constraint Energy Minimization Generalized Multiscale Finite Element Method (CEM-GMsFEM) for solving single-phase non-linear compressible flows in highly heterogeneous media. The construction of CEM-GMsFEM hinges on two crucial steps: First, the auxiliary space is constructed by solving local spectral problems, where the basis functions corresponding to small eig… ▽ More This paper presents and analyses a Constraint Energy Minimization Generalized Multiscale Finite Element Method (CEM-GMsFEM) for solving single-phase non-linear compressible flows in highly heterogeneous media. The construction of CEM-GMsFEM hinges on two crucial steps: First, the auxiliary space is constructed by solving local spectral problems, where the basis functions corresponding to small eigenvalues are captured. Then the basis functions are obtained by solving local energy minimization problems over the oversampling domains using the auxiliary space. The basis functions have exponential decay outside the corresponding local oversampling regions. The convergence of the proposed method is provided, and we show that this convergence only depends on the coarse grid size and is independent of the heterogeneities. An online enrichment guided by \emph{a posteriori} error estimator is developed to enhance computational efficiency. Several numerical experiments on a three-dimensional case to confirm the theoretical findings are presented, illustrating the performance of the method and giving efficient and accurate numerical. △ Less

Submitted 30 March, 2023; originally announced March 2023.

MSC Class: 65M12; 65M60; 65M22

Automated deep learning segmentation of high-resolution 7 T postmortem MRI for quantitative analysis of structure-pathology correlations in neurodegenerative diseases

Authors: Pulkit Khandelwal, Michael Tran Duong, Shokufeh Sadaghiani, Sydney Lim, Amanda Denning, Eunice Chung, Sadhana Ravikumar, Sanaz Arezoumandan, Claire Peterson, Madigan Bedard, Noah Capp, Ranjit Ittyerah, Elyse Migdal, Grace Choi, Emily Kopp, Bridget Loja, Eusha Hasan, Jiacheng Li, Alejandra Bahena, Karthik Prabhakaran, Gabor Mizsei, Marianna Gabrielyan, Theresa Schuck, Winifred Trotman, John Robinson , et al. (12 additional authors not shown)

Abstract: Postmortem MRI allows brain anatomy to be examined at high resolution and to link pathology measures with morphometric measurements. However, automated segmentation methods for brain mapping in postmortem MRI are not well developed, primarily due to limited availability of labeled datasets, and heterogeneity in scanner hardware and acquisition protocols. In this work, we present a high resolution… ▽ More Postmortem MRI allows brain anatomy to be examined at high resolution and to link pathology measures with morphometric measurements. However, automated segmentation methods for brain mapping in postmortem MRI are not well developed, primarily due to limited availability of labeled datasets, and heterogeneity in scanner hardware and acquisition protocols. In this work, we present a high resolution of 135 postmortem human brain tissue specimens imaged at 0.3 mm$^{3}$ isotropic using a T2w sequence on a 7T whole-body MRI scanner. We developed a deep learning pipeline to segment the cortical mantle by benchmarking the performance of nine deep neural architectures, followed by post-hoc topological correction. We then segment four subcortical structures (caudate, putamen, globus pallidus, and thalamus), white matter hyperintensities, and the normal appearing white matter. We show generalizing capabilities across whole brain hemispheres in different specimens, and also on unseen images acquired at 0.28 mm^3 and 0.16 mm^3 isotropic T2*w FLASH sequence at 7T. We then compute localized cortical thickness and volumetric measurements across key regions, and link them with semi-quantitative neuropathological ratings. Our code, Jupyter notebooks, and the containerized executables are publicly available at: https://pulkit-khandelwal.github.io/exvivo-brain-upenn △ Less

Submitted 17 October, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: Preprint submitted to NeuroImage Project website: https://pulkit-khandelwal.github.io/exvivo-brain-upenn

First BISTRO observations of the dark cloud Taurus L1495A-B10: the role of the magnetic field in the earliest stages of low-mass star formation

Authors: Derek Ward-Thompson, Janik Karoly, Kate Pattle, Anthony Whitworth, Jason Kirk, David Berry, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to… ▽ More We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965). △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages, 5 figures. ApJ accepted

With Shared Microexponents, A Little Shifting Goes a Long Way

Authors: Bita Rouhani, Ritchie Zhao, Venmugil Elango, Rasoul Shafipour, Mathew Hall, Maral Mesmakhosroshahi, Ankit More, Levi Melnick, Maximilian Golub, Girish Varatkar, Lei Shao, Gaurav Kolhe, Dimitry Melts, Jasmine Klar, Renee L'Heureux, Matt Perry, Doug Burger, Eric Chung, Zhaoxia Deng, Sam Naghshineh, Jongsoo Park, Maxim Naumov

Abstract: This paper introduces Block Data Representations (BDR), a framework for exploring and evaluating a wide spectrum of narrow-precision formats for deep learning. It enables comparison of popular quantization standards, and through BDR, new formats based on shared microexponents (MX) are identified, which outperform other state-of-the-art quantization approaches, including narrow-precision floating-p… ▽ More This paper introduces Block Data Representations (BDR), a framework for exploring and evaluating a wide spectrum of narrow-precision formats for deep learning. It enables comparison of popular quantization standards, and through BDR, new formats based on shared microexponents (MX) are identified, which outperform other state-of-the-art quantization approaches, including narrow-precision floating-point and block floating-point. MX utilizes multiple levels of quantization scaling with ultra-fine scaling factors based on shared microexponents in the hardware. The effectiveness of MX is demonstrated on real-world models including large-scale generative pretraining and inferencing, and production-scale recommendation systems. △ Less

Submitted 12 April, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

IC 5146 dark Streamer: is a first reliable candidate of edge collapse, hub-filament systems, and intertwined sub-filaments?

Authors: L. K. Dewangan, N. K. Bhadari, A. Men'shchikov, E. J. Chung, R. Devaraj, C. W. Lee, A. K. Maity, T. Baug

Abstract: The paper presents an analysis of multi-wavelength data of a nearby star-forming site IC 5146 dark Streamer (d $\sim$600 pc), which has been treated as a single and long filament, fl. Two hub-filament systems (HFSs) are known toward the eastern and the western ends of fl. Earlier published results favor the simultaneous evidence of HFSs and the end-dominated collapse (EDC) in fl. Herschel column d… ▽ More The paper presents an analysis of multi-wavelength data of a nearby star-forming site IC 5146 dark Streamer (d $\sim$600 pc), which has been treated as a single and long filament, fl. Two hub-filament systems (HFSs) are known toward the eastern and the western ends of fl. Earlier published results favor the simultaneous evidence of HFSs and the end-dominated collapse (EDC) in fl. Herschel column density map (resolution $\sim$13$''$.5) reveals two intertwined sub-filaments (i.e., fl-A and fl-B) toward fl, displaying a nearly double helix-like structure. This picture is also supported by the C$^{18}$O(3-2) emission. The scenario "fray and fragment" may explain the origin of intertwined sub-filaments. In the direction of fl, two cloud components around 2 and 4 km s$^{-1}$ are depicted using the $^{13}$CO(1-0) and C$^{18}$O(1-0) emission, and are connected in velocity space. The HFSs are spatially found at the overlapping areas of these cloud components and can be explained by the cloud-cloud collision scenario. Non-thermal gas motion in fl with larger Mach number is found. The magnetic field position angle measured from the filament's long axis shows a linear trend along the filament. This signature is confirmed in the other nearby EDC filaments, presenting a more quantitative confirmation of the EDC scenario. Based on our observational outcomes, we witness multiple processes operational in IC 5146 Streamer. Overall, the Streamer can be recognized as the first reliable candidate of edge collapse, HFSs, and intertwined sub-filaments together. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 22 pages, 11 figures, Accepted for publication in The Astrophysical Journal

JCMT BISTRO Observations: Magnetic Field Morphology of Bubbles Associated with NGC 6334

Authors: Mehrnoosh Tahani, Pierre Bastien, Ray S. Furuya, Kate Pattle, Doug Johnstone, Doris Arzoumanian, Yasuo Doi, Tetsuo Hasegawa, Shu-ichiro Inutsuka, Simon Coudé, Laura Fissel, Michael Chun-Yuan Chen, Frédérick Poidevin, Sarah Sadavoy, Rachel Friesen, Patrick M. Koch, James Di Francesco, Gerald H. Moriarty-Schieven, Zhiwei Chen, Eun Jung Chung, Chakali Eswaraiah, Lapo Fanciullo, Tim Gledhill, Valentin J. M. Le Gouellec, Thiem Hoang , et al. (120 additional authors not shown)

Abstract: We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles… ▽ More We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles indicate that the gas and magnetic field lines have been pushed away from the bubble, toward an almost tangential (to the bubble) magnetic field morphology. In the densest part of NGC 6334, where the magnetic field morphology is similar to an hourglass, the polarization observations do not exhibit observable impact from HII regions. We detect two nested radial polarization patterns in a bubble to the south of NGC 6334 that correspond to the previously observed bipolar structure in this bubble. Finally, using the results of this study, we present steps (incorporating computer vision; circular Hough Transform) that can be used in future studies to identify bubbles that have physically impacted magnetic field lines. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astrophysical Journal (ApJ)

Physics-preserving IMPES based multiscale methods for immiscible two-phase flow in highly heterogeneous porous media

Authors: Yiran Wang, Eric Chung, Shuyu Sun

Abstract: In this paper, we propose a physics-preserving multiscale method to solve an immiscible two-phase flow problem, which is modeled as a coupling system consisting of Darcy's law and mass conservation equations. We use a new Physics-preserving IMplicit Pressure Explicit Saturation (P-IMPES) scheme in order to maintain the local conservation of mass for both phases. Besides, this scheme is unbiased an… ▽ More In this paper, we propose a physics-preserving multiscale method to solve an immiscible two-phase flow problem, which is modeled as a coupling system consisting of Darcy's law and mass conservation equations. We use a new Physics-preserving IMplicit Pressure Explicit Saturation (P-IMPES) scheme in order to maintain the local conservation of mass for both phases. Besides, this scheme is unbiased and if the time step is smaller than a certain value, the saturation of both phases are bounds-preserving. When updating velocity, MGMsFEM serves as an efficient solver by computing the unknowns on a coarse grid. We follow the operation splitting techinque to deal with the two-phase flow. In particular, we use an upwind strategy to iterate the saturation explicitly and the MGMsFEM is utilized to compute velocity with a decoupled system on a coarse mesh. To show the efficiency and robustness of the proposed method, we design a set of interesting experiments. A rigorous analysis is also included to serve as a theoretical base of the method, which is well verified by the numerical results. Both simulations and analysis indicate that the method attains a good balance between accuracy and computation cost. △ Less

Submitted 12 December, 2022; originally announced December 2022.

Comments: 38 pages

The JCMT BISTRO-2 Survey: Magnetic Fields of the Massive DR21 Filament

Authors: Tao-Chung Ching, Keping Qiu, Di Li, Zhiyuan Ren, Shih-Ping Lai, David Berry, Kate Pattle, Ray Furuya, Derek Ward-Thompson, Doug Johnstone, Patrick M. Koch, Chang Won Lee, Thiem Hoang, Tetsuo Hasegawa, Woojin Kwon, Pierre Bastien, Chakali Eswaraiah, Jia-Wei Wang, Kyoung Hee Kim, Jihye Hwang, Archana Soam, A-Ran Lyo, Junhao Liu, Valentin J. M. Le Gouellec, Doris Arzoumanian , et al. (132 additional authors not shown)

Abstract: We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parall… ▽ More We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1--10 pc scales. The magnetic fields revealed in the Planck data are well aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6--1.0 mG in the DR21 filament and $\sim$ 0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by JCMT. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and sub-filaments. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 26 pages, 13 figures, ApJ accepted

A multiscale method for inhomogeneous elastic problems with high contrast coefficients

Authors: Zhongqian Wang, Changqing Ye, Eric T. Chung

Abstract: In this paper, we develop the constrained energy minimizing generalized multiscale finite element method (CEM-GMsFEM) with mixed boundary conditions (Dirichlet and Neumann) for the elasticity equations in high contrast media. By a special treatment of mixed boundary conditions separately, and combining the construction of the relaxed and constraint version of the CEM-GMsFEM, we discover that the m… ▽ More In this paper, we develop the constrained energy minimizing generalized multiscale finite element method (CEM-GMsFEM) with mixed boundary conditions (Dirichlet and Neumann) for the elasticity equations in high contrast media. By a special treatment of mixed boundary conditions separately, and combining the construction of the relaxed and constraint version of the CEM-GMsFEM, we discover that the method offers some advantages such as the independence of the target region's contrast from precision, while the sizes of oversampling domains have a significant impact on numerical accuracy. Moreover, to our best knowledge, this is the first proof of the convergence of the CEM-GMsFEM with mixed boundary conditions for the elasticity equations given. Some numerical experiments are provided to demonstrate the method's performance. △ Less

Submitted 20 October, 2022; originally announced October 2022.

Comments: CEM-GMsFEM; mixed boundary conditions; high contrast media

The JCMT BISTRO Survey: A Spiral Magnetic Field in a Hub-filament Structure, Monoceros R2

Authors: Jihye Hwang, Jongsoo Kim, Kate Pattle, Chang Won Lee, Patrick M. Koch, Doug Johnstone, Kohji Tomisaka, Anthony Whitworth, Ray S. Furuya, Ji-hyun Kang, A-Ran Lyo, Eun Jung Chung, Doris Arzoumanian, Geumsook Park, Woojin Kwon, Shinyoung Kim, Motohide Tamura, Jungmi Kwon, Archana Soam, Ilseung Han, Thiem Hoang, Kyoung Hee Kim, Takashi Onaka, Eswaraiah Chakali, Derek Ward-Thompson , et al. (135 additional authors not shown)

Abstract: We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R… ▽ More We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well-described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis-Chandrasekhar-Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from $Herschel$ data and the C$^{18}$O ($J$ = 3-2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 $\pm$ 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 $\pm$ 0.02. Additionally, the mean Alfvén Mach number is 0.35 $\pm$ 0.01. This suggests that in Mon R2, magnetic fields provide resistance against large-scale gravitational collapse, and magnetic pressure exceeds turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically sub-critical. △ Less

Submitted 13 December, 2022; v1 submitted 12 October, 2022; originally announced October 2022.

Comments: This paper is accepted to the ApJ

Role of Filamentary Structures in the Formation of Two Dense Cores, L1544 and L694-2

Authors: Shinyoung Kim, Chang Won Lee, Mario Tafalla, Maheswar Gophinathan, Paola Caselli, Philip C. Myers, Eun Jung Chung, Shanghuo Li

Abstract: We present mapping results of two prestellar cores, L1544 and L694-2, embedded in filamentary clouds in C$^{18}$O (3-2), $^{13}$CO (3-2), $^{12}$CO (3-2), HCO$^+$ (4-3), and H$^{13}$CO$^+$ (4-3) with the JCMT telescope to examine the role of the filamentary structures in the formation of dense cores in the clouds, with new distance estimates for L1544 ($175_{-3}^{+4}$ pc) and L694-2 (… ▽ More We present mapping results of two prestellar cores, L1544 and L694-2, embedded in filamentary clouds in C$^{18}$O (3-2), $^{13}$CO (3-2), $^{12}$CO (3-2), HCO$^+$ (4-3), and H$^{13}$CO$^+$ (4-3) with the JCMT telescope to examine the role of the filamentary structures in the formation of dense cores in the clouds, with new distance estimates for L1544 ($175_{-3}^{+4}$ pc) and L694-2 ($203_{-7}^{+6}$ pc). From these observations, we found that the non-thermal velocity dispersion of two prestellar cores and their surrounding clouds are smaller than or comparable to the sound speed. This may indicate that the turbulence has already been dissipated for both filaments and cores during their formation time. We also found a $λ/4$ shift between the periodic oscillations in the velocity and the column density distributions implying the possible presence of gravitational core-forming flow motion along the axis of the filament. The mass accretion rates due to these flow motions are estimated to be 2-3 M$_\odot$ Myr$^{-1}$, being comparable to that for Serpens cloud but much smaller than those for the Hub filaments, cluster, or high mass forming filaments by 1 or 2 order of magnitudes. From this study, we suggest that the filaments in our targets might be formed from the shock-compression of colliding clouds, and then the cores are formed by gravitational fragmentation of the filaments to evolve to the prestellar stage. We conclude that the filamentary structures in the clouds play an important role in the entire process of formation of dense cores and their evolution. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 30 pages, 22 figures, 6 tables, accepted for publication in ApJ

The JCMT BISTRO Survey: Multi-wavelength polarimetry of bright regions in NGC 2071 in the far-infrared/submillimetre range, with POL-2 and HAWC+

Authors: L. Fanciullo, F. Kemper, K. Pattle, P. M. Koch, S. Sadavoy, S. Coudé, A. Soam, T. Hoang, T. Onaka, V. J. M. Le Gouellec, D. Arzoumanian, D. Berry, C. Eswaraiah, E. J. Chung, R. Furuya, C. L. H. Hull, J. Hwang, D. Johnstone, J. -h. Kang, K. H. Kim, F. Kirchschlager, V. Könyves, J. Kwon, W. Kwon, S. -P. Lai , et al. (9 additional authors not shown)

Abstract: Polarized dust emission is a key tracer in the study of interstellar medium and of star formation. The observed polarization, however, is a product of magnetic field structure, dust grain properties and grain alignment efficiency, as well as their variations in the line of sight, making it difficult to interpret polarization unambiguously. The comparison of polarimetry at multiple wavelengths is a… ▽ More Polarized dust emission is a key tracer in the study of interstellar medium and of star formation. The observed polarization, however, is a product of magnetic field structure, dust grain properties and grain alignment efficiency, as well as their variations in the line of sight, making it difficult to interpret polarization unambiguously. The comparison of polarimetry at multiple wavelengths is a possible way of mitigating this problem. We use data from HAWC+/SOFIA and from SCUBA-2/POL-2 (from the BISTRO survey) to analyse the NGC 2071 molecular cloud at 154, 214 and 850 $μ$m. The polarization angle changes significantly with wavelength over part of NGC 2071, suggesting a change in magnetic field morphology on the line of sight as each wavelength best traces different dust populations. Other possible explanations are the existence of more than one polarization mechanism in the cloud or scattering from very large grains. The observed change of polarization fraction with wavelength, and the 214-to-154 $μ$m polarization ratio in particular, are difficult to reproduce with current dust models under the assumption of uniform alignment efficiency. We also show that the standard procedure of using monochromatic intensity as a proxy for column density may produce spurious results at HAWC+ wavelengths. Using both long-wavelength (POL-2, 850 $μ$m) and short-wavelength (HAWC+, $\lesssim 200\, μ$m) polarimetry is key in obtaining these results. This study clearly shows the importance of multi-wavelength polarimetry at submillimeter bands to understand the dust properties of molecular clouds and the relationship between magnetic field and star formation. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: Main article: 18 pages, 11 figures. Online supplemental material: 2 pages, 3 figures

Journal ref: Monthly Notices of the Royal Astronomical Society vol. 512 (2022) pp. 1985-2002

Evolution of the Hub-filament Structures in IC 5146 in the Context of the Energy Balance of Gravity, Turbulence, and Magnetic Field

Authors: Eun Jung Chung, Chang Won Lee, Woojin Kwon, Hyunju Yoo, Archana Soam, Jungyeon Cho

Abstract: We present the results of 850 $μ$m polarization and C$^{18}$O (3-2) line observations toward the western hub-filament structure (W-HFS) of the dark Streamer in IC 5146 using the James Clerk Maxwell Telescope (JCMT) SCUBA-2/POL-2 and HARP instruments. We aim to investigate how the relative importance of the magnetic field, gravity, and turbulence affects core formation in HFS by comparing the energ… ▽ More We present the results of 850 $μ$m polarization and C$^{18}$O (3-2) line observations toward the western hub-filament structure (W-HFS) of the dark Streamer in IC 5146 using the James Clerk Maxwell Telescope (JCMT) SCUBA-2/POL-2 and HARP instruments. We aim to investigate how the relative importance of the magnetic field, gravity, and turbulence affects core formation in HFS by comparing the energy budget of this region. We identified four 850 $μ$m cores and estimated the magnetic field strengths ($B_{\rm pos}$) of the cores and the hub and filament using the Davis-Chandrasekhar-Fermi method. The estimated $B_{\rm pos}$ is $\sim$80 to 1200 $μ$G. From Wang et al., $B_{\rm pos}$ of E-47, a core in the eastern hub (E-hub), and E-hub were re-estimated to be 500 and 320 $μ$G, respectively, with the same method. We measured the gravitational ($E_{\rm G}$), kinematic ($E_{\rm K}$), and magnetic energies ($E_{\rm B}$) in the filament and hubs and compared the relative importance among them. We found that an $E_{\rm B}$-dominant filament has $aligned$ fragmentation type, while $E_{\rm G}$-dominant hubs show $no$ and $clustered$ fragmentation types. In the $E_{\rm G}$ dominant hubs, it seems that the portion of $E_{\rm K}$ determines whether the hub becomes to have $clustered$ (the portion of $E_{\rm K}\sim20\%$) or $no$ fragmentation type ($\sim10\%$). We propose an evolutionary scenario for the E- and W-HFSs, where the HFS forms first by the collision of turbulent flows, and then the hubs and filaments can go into various types of fragmentation depending on their energy balance of gravity, turbulence, and magnetic field. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: 19 pages, 10 figures, accepted for publication in AJ

A discontinuous Galerkin based multiscale method for heterogeneous elastic wave equations

Authors: Zhongqian Wang, Shubin Fu, Zishang Li, Eric Chung

Abstract: In this paper, we develop a local multiscale model reduction strategy for the elastic wave equation in strongly heterogeneous media, which is achieved by solving the problem in a coarse mesh with multiscale basis functions. We use the interior penalty discontinuous Galerkin (IPDG) to couple the multiscale basis functions that contain important heterogeneous media information. The construction of e… ▽ More In this paper, we develop a local multiscale model reduction strategy for the elastic wave equation in strongly heterogeneous media, which is achieved by solving the problem in a coarse mesh with multiscale basis functions. We use the interior penalty discontinuous Galerkin (IPDG) to couple the multiscale basis functions that contain important heterogeneous media information. The construction of efficient multiscale basis functions starts with extracting dominant modes of carefully defined spectral problems to represent important media feature, which is followed by solving a constraint energy minimization problems. Then a Petrov-Galerkin projection and systematization onto the coarse grid is applied. As a result, an explicit and energy conserving scheme is obtained for fast online simulation. The method exhibits both coarse-mesh and spectral convergence as long as one appropriately chose the oversampling size. We rigorously analyze the stability and convergence of the proposed method. Numerical results are provided to show the performance of the multiscale method and confirm the theoretical results. △ Less

Submitted 10 July, 2022; originally announced July 2022.

AI-based automated Meibomian gland segmentation, classification and reflection correction in infrared Meibography

Authors: Ripon Kumar Saha, A. M. Mahmud Chowdhury, Kyung-Sun Na, Gyu Deok Hwang, Youngsub Eom, Jaeyoung Kim, Hae-Gon Jeon, Ho Sik Hwang, Euiheon Chung

Abstract: Purpose: Develop a deep learning-based automated method to segment meibomian glands (MG) and eyelids, quantitatively analyze the MG area and MG ratio, estimate the meiboscore, and remove specular reflections from infrared images. Methods: A total of 1600 meibography images were captured in a clinical setting. 1000 images were precisely annotated with multiple revisions by investigators and graded… ▽ More Purpose: Develop a deep learning-based automated method to segment meibomian glands (MG) and eyelids, quantitatively analyze the MG area and MG ratio, estimate the meiboscore, and remove specular reflections from infrared images. Methods: A total of 1600 meibography images were captured in a clinical setting. 1000 images were precisely annotated with multiple revisions by investigators and graded 6 times by meibomian gland dysfunction (MGD) experts. Two deep learning (DL) models were trained separately to segment areas of the MG and eyelid. Those segmentation were used to estimate MG ratio and meiboscores using a classification-based DL model. A generative adversarial network was implemented to remove specular reflections from original images. Results: The mean ratio of MG calculated by investigator annotation and DL segmentation was consistent 26.23% vs 25.12% in the upper eyelids and 32.34% vs. 32.29% in the lower eyelids, respectively. Our DL model achieved 73.01% accuracy for meiboscore classification on validation set and 59.17% accuracy when tested on images from independent center, compared to 53.44% validation accuracy by MGD experts. The DL-based approach successfully removes reflection from the original MG images without affecting meiboscore grading. Conclusions: DL with infrared meibography provides a fully automated, fast quantitative evaluation of MG morphology (MG Segmentation, MG area, MG ratio, and meiboscore) which are sufficiently accurate for diagnosing dry eye disease. Also, the DL removes specular reflection from images to be used by ophthalmologists for distraction-free assessment. △ Less

Submitted 31 May, 2022; originally announced May 2022.

Comments: 11 pages, 13 Figures, 5 Supplementary Figures

Homogenization with quasistatic Tresca's friction law: qualitative and quantitative results

Authors: Changqing Ye, Eric T. Chung, Junzhi Cui

Abstract: Modeling of frictional contacts is crucial for investigating mechanical performances of composite materials under varying service environments. The paper considers a linear elasticity system with strongly heterogeneous coefficients and quasistatic Tresca friction law, and studies the homogenization theories under the frameworks of H-convergence and small $ε$-periodicity. The qualitative result is… ▽ More Modeling of frictional contacts is crucial for investigating mechanical performances of composite materials under varying service environments. The paper considers a linear elasticity system with strongly heterogeneous coefficients and quasistatic Tresca friction law, and studies the homogenization theories under the frameworks of H-convergence and small $ε$-periodicity. The qualitative result is based on H-convergence, which shows the original oscillating solutions will converge weakly to the homogenized solution, while our quantitative result provides an estimate of asymptotic errors in $H^1$-norm for the periodic homogenization. This paper also designs several numerical experiments to validate the convergence rates in the quantitative analysis. △ Less

Submitted 20 July, 2023; v1 submitted 20 April, 2022; originally announced April 2022.

MSC Class: 74Q05; 74Q10; 35J86

Local multiscale model reduction using discontinuous Galerkin coupling for elasticity problems

Authors: Zhongqian Wang, Shubin Fu, Eric Chung

Abstract: In this paper, we consider the constrained energy minimizing generalized multiscale finite element method (CEM-GMsFEM) with discontinuous Galerkin (DG) coupling for the linear elasticity equations in highly heterogeneous and high contrast media. We will introduce the construction of a DG version of the CEM-GMsFEM, such as auxiliary basis functions and offline basis functions. The DG version of the… ▽ More In this paper, we consider the constrained energy minimizing generalized multiscale finite element method (CEM-GMsFEM) with discontinuous Galerkin (DG) coupling for the linear elasticity equations in highly heterogeneous and high contrast media. We will introduce the construction of a DG version of the CEM-GMsFEM, such as auxiliary basis functions and offline basis functions. The DG version of the method offers some advantages such as flexibility in coarse grid construction and sparsity of resulting discrete systems. Moreover, to our best knowledge, this is the first time where the proof of the convergence of the CEM-GMsFEM in the DG form is given. Some numerical examples will be presented to illustrate the performance of the method. △ Less

Submitted 10 October, 2022; v1 submitted 16 April, 2022; originally announced April 2022.

Effects of magnetic field orientations in dense cores on gas kinematics in protostellar envelopes

Authors: Aashish Gupta, Hsi-Wei Yen, Patrick Koch, Pierre Bastien, Tyler L. Bourke, Eun Jung Chung, Tetsuo Hasegawa, Charles L. H. Hull, Shu-ichiro Inutsuka, Jungmi Kwon, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Chin-Fei Lee, Kate Pattle, Keping Qiu, Mehrnoosh Tahani, Motohide Tamura, Derek Ward-Thompson

Abstract: Theoretically, misalignment between the magnetic field and rotational axis in a dense core is considered to be dynamically important in the star formation process, however, extent of this influence remains observationally unclear. For a sample of 32 Class 0 and I protostars in the Perseus Molecular Cloud, we analyzed gas motions using C$^{18}$O data from the SMA MASSES survey and the magnetic fiel… ▽ More Theoretically, misalignment between the magnetic field and rotational axis in a dense core is considered to be dynamically important in the star formation process, however, extent of this influence remains observationally unclear. For a sample of 32 Class 0 and I protostars in the Perseus Molecular Cloud, we analyzed gas motions using C$^{18}$O data from the SMA MASSES survey and the magnetic field structures using 850 $μ$m polarimetric data from the JCMT BISTRO-1 survey and archive. We do not find any significant correlation between the velocity gradients in the C$^{18}$O emission in the protostellar envelopes at a 1,000 au scale and the misalignment between the outflows and magnetic field orientations in the dense cores at a 4,000 au scale, and there is also no correlation between the velocity gradients and the angular dispersions of the magnetic fields. However, a significant dependence on the misalignment angles emerges after we normalize the rotational motion by the infalling motion, where the ratios increase from $\lesssim1$ to $\gtrsim1$ with increasing misalignment angles. This suggests that the misalignment could prompt angular momentum transportation to the envelope scale but is not a dominant factor in determining the envelope rotation, and other parameters, like mass accretion in protostellar sources, also play an important role. These results remain valid after taking into account projection effects. The comparison between our estimated angular momentum in the protostellar envelopes and the sizes of the known protostellar disks suggests that significant angular momentum is likely lost between radii of $\sim$1,000-100 au in protostellar envelopes. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Comments: Accepted by ApJ

Journal ref: ApJ 930 67 (2022)

Constraint energy minimizing generalized multiscale finite element method for convection diffusion equation

Authors: Lina Zhao, Eric Chung

Abstract: In this paper we present and analyze a constraint energy minimizing generalized multiscale finite element method for convection diffusion equation. To define the multiscale basis functions, we first build an auxiliary multiscale space by solving local spectral problems motivated by analysis. Then constraint energy minimization performed in oversampling domains is exploited to construct the multisc… ▽ More In this paper we present and analyze a constraint energy minimizing generalized multiscale finite element method for convection diffusion equation. To define the multiscale basis functions, we first build an auxiliary multiscale space by solving local spectral problems motivated by analysis. Then constraint energy minimization performed in oversampling domains is exploited to construct the multiscale space. The resulting multiscale basis functions have a good decay property even for high contrast diffusion and convection coefficients. Furthermore, if the number of oversampling layer is chosen properly, we can prove that the convergence rate is proportional to the coarse mesh size. Our analysis also indicates that the size of the oversampling domain weakly depends on the contrast of the heterogeneous coefficients. Several numerical experiments are presented illustrating the performances of our method. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: arXiv admin note: text overlap with arXiv:1704.03193

A conservative multiscale method for stochastic highly heterogeneous flow

Authors: Yiran Wang, Eric Chung, Shubin Fu

Abstract: In this paper, we propose a local model reduction approach for subsurface flow problems in stochastic and highly heterogeneous media. To guarantee the mass conservation, we consider the mixed formulation of the flow problem and aim to solve the problem in a coarse grid to reduce the complexity of a large-scale system. We decompose the entire problem into a training and a testing stage, namely the… ▽ More In this paper, we propose a local model reduction approach for subsurface flow problems in stochastic and highly heterogeneous media. To guarantee the mass conservation, we consider the mixed formulation of the flow problem and aim to solve the problem in a coarse grid to reduce the complexity of a large-scale system. We decompose the entire problem into a training and a testing stage, namely the offline coarse-grid multiscale basis generation stage and online simulation stage with different parameters. In the training stage, a parameter-independent and small-dimensional multiscale basis function space is constructed, which includes the media, source and boundary information. The key part of the basis generation stage is to solve some local problems defined specially. With the parameter-independent basis space, one can efficiently solve the concerned problems corresponding to different samples of permeability field in a coarse grid without repeatedly constructing a multiscale space for each new sample. A rigorous analysis on convergence of the proposed method is proposed. In particular, we consider a generalization error, where bases constructed with one source will be used to a different source. In the numerical experiments, we apply the proposed method for both single-phase and twophase flow problems. Simulation results for both 2D and 3D representative models demonstrate the high accuracy and impressive performance of the proposed model reduction techniques. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 29 pages

Stronger 3SUM-Indexing Lower Bounds

Authors: Eldon Chung, Kasper Green Larsen

Abstract: The $3$SUM-Indexing problem was introduced as a data structure version of the $3$SUM problem, with the goal of proving strong conditional lower bounds for static data structures via reductions. Ideally, the conjectured hardness of $3$SUM-Indexing should be replaced by an unconditional lower bound. Unfortunately, we are far from proving this, with the strongest current lower bound being a logarithm… ▽ More The $3$SUM-Indexing problem was introduced as a data structure version of the $3$SUM problem, with the goal of proving strong conditional lower bounds for static data structures via reductions. Ideally, the conjectured hardness of $3$SUM-Indexing should be replaced by an unconditional lower bound. Unfortunately, we are far from proving this, with the strongest current lower bound being a logarithmic query time lower bound by Golovnev et al. from STOC'20. Moreover, their lower bound holds only for non-adaptive data structures and they explicitly asked for a lower bound for adaptive data structures. Our main contribution is precisely such a lower bound against adaptive data structures. As a secondary result, we also strengthen the non-adaptive lower bound of Golovnev et al. and prove strong lower bounds for $2$-bit-probe non-adaptive $3$SUM-Indexing data structures via a completely new approach that we find interesting in its own right. △ Less

Submitted 25 March, 2023; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: SODA 2023

MSC Class: 68P05 ACM Class: E.1; F.2.0

Generalized multiscale finite element method for highly heterogeneous compressible flow

Authors: Shubin Fu, Eric Chung, Lina Zhao

Abstract: In this paper, we study the generalized multiscale finite element method (GMsFEM) for single phase compressible flow in highly heterogeneous porous media. We follow the major steps of the GMsFEM to construct permeability dependent offline basis for fast coarse-grid simulation. The offline coarse space is efficiently constructed only once based on the initial permeability field with parallel comput… ▽ More In this paper, we study the generalized multiscale finite element method (GMsFEM) for single phase compressible flow in highly heterogeneous porous media. We follow the major steps of the GMsFEM to construct permeability dependent offline basis for fast coarse-grid simulation. The offline coarse space is efficiently constructed only once based on the initial permeability field with parallel computing. A rigorous convergence analysis is performed for two types of snapshot spaces. The analysis indicates that the convergence rates of the proposed multiscale method depend on the coarse meshsize and the eigenvalue decay of the local spectral problem. To further increase the accuracy of multiscale method, residual driven online multiscale basis is added to the offline space. The construction of online multiscale basis is based on a carefully design error indicator motivated by the analysis. We find that online basis is particularly important for the singular source. Rich numerical tests on typical 3D highly heterogeneous medias are presented to demonstrate the impressive computational advantages of the proposed multiscale method. △ Less

Submitted 18 January, 2022; originally announced January 2022.