Showing 1–33 of 33 results for author: Lee, C T

Spatial modeling algorithms for reactions and transport (SMART) in biological cells

Authors: Emmet A. Francis, Justin G. Laughlin, Jørgen S. Dokken, Henrik N. T. Finsberg, Christopher T. Lee, Marie E. Rognes, Padmini Rangamani

Abstract: Biological cells rely on precise spatiotemporal coordination of biochemical reactions to control their many functions. Such cell signaling networks have been a common focus for mathematical models, but they remain challenging to simulate, particularly in realistic cell geometries. Herein, we present our software, Spatial Modeling Algorithms for Reactions and Transport (SMART), a package that takes… ▽ More Biological cells rely on precise spatiotemporal coordination of biochemical reactions to control their many functions. Such cell signaling networks have been a common focus for mathematical models, but they remain challenging to simulate, particularly in realistic cell geometries. Herein, we present our software, Spatial Modeling Algorithms for Reactions and Transport (SMART), a package that takes in high-level user specifications about cell signaling networks and molecular transport, and then assembles and solves the associated mathematical and computational systems. SMART uses state-of-the-art finite element analysis, via the FEniCS Project software, to efficiently and accurately resolve cell signaling events over discretized cellular and subcellular geometries. We demonstrate its application to several different biological systems, including YAP/TAZ mechanotransduction, calcium signaling in neurons and cardiomyocytes, and ATP generation in mitochondria. Throughout, we utilize experimentally-derived realistic cellular geometries represented by well-conditioned tetrahedral meshes. These scenarios demonstrate the applicability, flexibility, accuracy and efficiency of SMART across a range of temporal and spatial scales. △ Less

Submitted 24 May, 2024; originally announced May 2024.

SMART: Spatial Modeling Algorithms for Reaction and Transport

Authors: Justin G. Laughlin, Jørgen S. Dokken, Henrik N. T. Finsberg, Emmet A. Francis, Christopher T. Lee, Marie E. Rognes, Padmini Rangamani

Abstract: Recent advances in microscopy and 3D reconstruction methods have allowed for characterization of cellular morphology in unprecedented detail, including the irregular geometries of intracellular subcompartments such as membrane-bound organelles. These geometries are now compatible with predictive modeling of cellular function. Biological cells respond to stimuli through sequences of chemical reacti… ▽ More Recent advances in microscopy and 3D reconstruction methods have allowed for characterization of cellular morphology in unprecedented detail, including the irregular geometries of intracellular subcompartments such as membrane-bound organelles. These geometries are now compatible with predictive modeling of cellular function. Biological cells respond to stimuli through sequences of chemical reactions generally referred to as cell signaling pathways. The propagation and reaction of chemical substances in cell signaling pathways can be represented by coupled nonlinear systems of reaction-transport equations. These reaction pathways include numerous chemical species that react across boundaries or interfaces (e.g., the cell membrane and membranes of organelles within the cell) and domains (e.g., the bulk cell volume and the interior of organelles). Such systems of multi-dimensional partial differential equations (PDEs) are notoriously difficult to solve because of their high dimensionality, non-linearities, strong coupling, stiffness, and potential instabilities. In this work, we describe Spatial Modeling Algorithms for Reactions and Transport (SMART), a high-performance finite-element-based simulation package for model specification and numerical simulation of spatially-varying reaction-transport processes. SMART is based on the FEniCS finite element library, provides a symbolic representation framework for specifying reaction pathways, and supports geometries in 2D and 3D including large and irregular cell geometries obtained from modern ultrastructural characterization methods. △ Less

Submitted 12 June, 2023; originally announced June 2023.

Comments: 5 pages, 2 figures, submitted to the Journal of Open Source Software (JOSS), code available at https://github.com/RangamaniLabUCSD/smart

Mem3DG: Modeling Membrane Mechanochemical Dynamics in 3D using Discrete Differential Geometry

Authors: Cuncheng Zhu, Christopher T. Lee, Padmini Rangamani

Abstract: Biomembranes adopt varying morphologies that are vital to cellular functions. Many studies use computational modeling to understand how various mechanochemical factors contribute to membrane shape transformations. Compared to approximation-based methods (e.g., finite element method), the class of discrete mesh models offers greater flexibility to simulate complex physics and shapes in three dimens… ▽ More Biomembranes adopt varying morphologies that are vital to cellular functions. Many studies use computational modeling to understand how various mechanochemical factors contribute to membrane shape transformations. Compared to approximation-based methods (e.g., finite element method), the class of discrete mesh models offers greater flexibility to simulate complex physics and shapes in three dimensions; its formulation produces an efficient algorithm while maintaining coordinate-free geometric descriptions. However, ambiguities in geometric definitions in the discrete context have led to a lack of consensus on which discrete mesh model is theoretically and numerically optimal; a bijective relationship between the terms contributing to both the energy and forces from the discrete and smooth geometric theories remains to be established. We address this and present an extensible framework, $\texttt{Mem3DG}$, for modeling 3D mechanochemical dynamics of membranes based on Discrete Differential Geometry (DDG) on triangulated meshes. The formalism of DDG resolves the inconsistency and provides a unifying perspective on how to relate the smooth and discrete energy and forces. To demonstrate, $\texttt{Mem3DG}$ is used to model a sequence of examples with increasing mechanochemical complexity: recovering classical shape transformations such as 1) biconcave disk, dumbbell, and unduloid and 2) spherical bud on spherical, flat-patch membrane; investigating how the coupling of membrane mechanics with protein mobility jointly affects phase and shape transformation. As high-resolution 3D imaging of membrane ultrastructure becomes more readily available, we envision Mem3DG to be applied as an end-to-end tool to simulate realistic cell geometry under user-specified mechanochemical conditions. △ Less

Submitted 31 October, 2021; originally announced November 2021.

Stellar Masses of Giant Clumps in CANDELS and Simulated Galaxies Using Machine Learning

Authors: M. Huertas-Company, Y. Guo, O. Ginzburg, C. T. Lee, N. Mandelker, M. Metter, J. R. Primack, A. Dekel, D. Ceverino, S. M. Faber, D. C. Koo, A. Koekemoer, G. Snyder, M. Giavalisco, H. Zhang

Abstract: A significant fraction of high redshift star-forming disc galaxies are known to host giant clumps, whose nature and role in galaxy evolution are yet to be understood. In this work we first present a new method based on neural networks to detect clumps in galaxy images. We use this method to detect clumps in the rest-frame optical and UV images of a complete sample of $\sim1500$ star forming galaxi… ▽ More A significant fraction of high redshift star-forming disc galaxies are known to host giant clumps, whose nature and role in galaxy evolution are yet to be understood. In this work we first present a new method based on neural networks to detect clumps in galaxy images. We use this method to detect clumps in the rest-frame optical and UV images of a complete sample of $\sim1500$ star forming galaxies at $1<z<3$ in the CANDELS survey as well as in images from the VELA zoom-in cosmological simulations. We show that observational effects have a dramatic impact on the derived clump properties leading to an overestimation of the clump mass up to a factor of 10, which highlights the importance of fair comparisons between observations and simulations and the limitations of current HST data to study the resolved structure of distant galaxies. After correcting for these effects with a mixture density network, we estimate that the clump stellar mass function follows a power-law down to the completeness limit ($10^{7}$ solar masses) with the majority of the clumps being less massive than $10^9$ solar masses. This is in better agreement with recent gravitational lensing based measurements. The simulations explored in this work overall reproduce the shape of the observed clump stellar mass function and clumpy fractions when confronted under the same conditions, although they tend to lie in the lower limit of the confidence intervals of the observations. This agreement suggests that most of the observed clumps are formed in-situ. △ Less

Submitted 8 September, 2020; v1 submitted 25 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS - This is the final version

Measurement of spin dynamics in a layered nickelate using x-ray photon correlation spectroscopy: Evidence for intrinsic destabilization of incommensurate stripes at low temperatures

Authors: Alessandro Ricci, Nicola Poccia, Gaetano Campi, Shrawan Mishra, Leonard Müller, Boby Joseph, Bo Shi, Alexey Zozulya, Marcel Buchholz, Christoph Trabant, James C. T. Lee, Jens Viefhaus, Jeroen B. Goedkoop, Agustinus Agung Nugroho, Markus Braden, Sujoy Roy, Michael Sprung, Christian Schüßler-Langeheine

Abstract: We study the temporal stability of stripe-type spin order in a layered nickelate with X-ray photon correlation spectroscopy and observe fluctuations on time scales of tens of minutes over a wide temperature range. These fluctuations show an anomalous temperature dependence: they slow down at intermediate temperatures and speed up both upon heating and cooling. This behavior appears to be directly… ▽ More We study the temporal stability of stripe-type spin order in a layered nickelate with X-ray photon correlation spectroscopy and observe fluctuations on time scales of tens of minutes over a wide temperature range. These fluctuations show an anomalous temperature dependence: they slow down at intermediate temperatures and speed up both upon heating and cooling. This behavior appears to be directly connected with spatial correlations: stripes fluctuate slowly when stripe correlation lengths are large and become faster when spatial correlations decrease. A low-temperature decay of nickelate stripe correlations, reminiscent of what occurs in cuprates due to a competition between stripes and superconductivity, hence occurs via loss of both spatial and temporal correlations. △ Less

Submitted 27 May, 2021; v1 submitted 16 December, 2019; originally announced December 2019.

Comments: Accepted for publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 127, 057001 (2021)

Applications and Challenges of Machine Learning to Enable Realistic Cellular Simulations

Authors: Ritvik Vasan, Meagan P. Rowan, Christopher T. Lee, Gregory R. Johnson, Padmini Rangamani, Michael Holst

Abstract: In this perspective, we examine three key aspects of an end-to-end pipeline for realistic cellular simulations: reconstruction and segmentation of cellular structures; generation of cellular structures; and mesh generation, simulation, and data analysis. We highlight some of the relevant prior work in these distinct but overlapping areas, with a particular emphasis on current use of machine learni… ▽ More In this perspective, we examine three key aspects of an end-to-end pipeline for realistic cellular simulations: reconstruction and segmentation of cellular structures; generation of cellular structures; and mesh generation, simulation, and data analysis. We highlight some of the relevant prior work in these distinct but overlapping areas, with a particular emphasis on current use of machine learning technologies, as well as on future opportunities. △ Less

Submitted 12 November, 2019; originally announced November 2019.

The SFR-radius connection: data and implications for wind strength and halo concentration

Authors: Lin Lin, S. M. Faber, David C. Koo, Samir Salim, Aaron A. Dutton, Jerome J. Fang, Fangzhou Jiang, Cristoph T. Lee, Aldo Rodríguez-Puebla, A. van der Wel, Yicheng Guo, Guillermo Barro, Joel R. Primack, Avishai Dekel, Zhu Chen, Yifei Luo, Viraj Pandya, Rachel S. Somerville, Henry C. Ferguson, Susan Kassin, Anton M. Koekemoer, Norman A. Grogin, Audrey Galametz, P. Santini, Hooshang Nayyeri , et al. (4 additional authors not shown)

Abstract: This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend… ▽ More This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR. △ Less

Submitted 20 July, 2020; v1 submitted 24 October, 2019; originally announced October 2019.

Comments: 21 pages, 9 figures, ApJ accepted

An Open Source Mesh Generation Platform for Biophysical Modeling Using Realistic Cellular Geometries

Authors: Christopher T. Lee, Justin G. Laughlin, John B. Moody, Rommie E. Amaro, J. Andrew McCammon, Michael J. Holst, Padmini Rangamani

Abstract: Advances in imaging methods such as electron microscopy, tomography and other modalities are enabling high-resolution reconstructions of cellular and organelle geometries. Such advances pave the way for using these geometries for biophysical and mathematical modeling once these data can be represented as a geometric mesh, which, when carefully conditioned, enables the discretization and solution o… ▽ More Advances in imaging methods such as electron microscopy, tomography and other modalities are enabling high-resolution reconstructions of cellular and organelle geometries. Such advances pave the way for using these geometries for biophysical and mathematical modeling once these data can be represented as a geometric mesh, which, when carefully conditioned, enables the discretization and solution of partial differential equations. In this study, we outline the steps for a naïve user to approach GAMer 2, a mesh generation code written in C++ designed to convert structural datasets to realistic geometric meshes, while preserving the underlying shapes. We present two example cases, 1) mesh generation at the subcellular scale as informed by electron tomography, and 2) meshing a protein with structure from x-ray crystallography. We further demonstrate that the meshes generated by GAMer are suitable for use with numerical methods. Together, this collection of libraries and tools simplifies the process of constructing realistic geometric meshes from structural biology data. △ Less

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

Comments: 6 pages and 4 figures. Supplemental Movie available upon request

3D mesh processing using GAMer 2 to enable reaction-diffusion simulations in realistic cellular geometries

Authors: Christopher T. Lee, Justin G. Laughlin, Nils Angliviel de La Beaumelle, Rommie E. Amaro, J. Andrew McCammon, Ravi Ramamoorthi, Michael J. Holst, Padmini Rangamani

Abstract: Recent advances in electron microscopy have enabled the imaging of single cells in 3D at nanometer length scale resolutions. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations requires watertight meshing of electron micrograph images into 3D volume meshes, which can then form the basis of computer sim… ▽ More Recent advances in electron microscopy have enabled the imaging of single cells in 3D at nanometer length scale resolutions. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations requires watertight meshing of electron micrograph images into 3D volume meshes, which can then form the basis of computer simulations of such processes using numerical techniques such as the Finite Element Method. In this paper, we describe the use of our recently rewritten mesh processing software, GAMer 2, to bridge the gap between poorly conditioned meshes generated from segmented micrographs and boundary marked tetrahedral meshes which are compatible with simulation. We demonstrate the application of a workflow using GAMer 2 to a series of electron micrographs of neuronal dendrite morphology explored at three different length scales and show that the resulting meshes are suitable for finite element simulations. This work is an important step towards making physical simulations of biological processes in realistic geometries routine. Innovations in algorithms to reconstruct and simulate cellular length scale phenomena based on emerging structural data will enable realistic physical models and advance discovery at the interface of geometry and cellular processes. We posit that a new frontier at the intersection of computational technologies and single cell biology is now open. △ Less

Submitted 17 December, 2019; v1 submitted 29 January, 2019; originally announced January 2019.

Comments: 39 pages, 14 figures. High resolution figures and supplemental movies available upon request

The Implementation of the Colored Abstract Simplicial Complex and its Application to Mesh Generation

Authors: C. T. Lee, J. B. Moody, R. E. Amaro, J. A. McCammon, M. Holst

Abstract: We introduce CASC: a new, modern, and header-only C++ library which provides a data structure to represent arbitrary dimension abstract simplicial complexes (ASC) with user-defined classes stored directly on the simplices at each dimension. This is accomplished by using the latest C++ language features including variadic template parameters introduced in C++11 and automatic function return type de… ▽ More We introduce CASC: a new, modern, and header-only C++ library which provides a data structure to represent arbitrary dimension abstract simplicial complexes (ASC) with user-defined classes stored directly on the simplices at each dimension. This is accomplished by using the latest C++ language features including variadic template parameters introduced in C++11 and automatic function return type deduction from C++14. Effectively CASC decouples the representation of the topology from the interactions of user data. We present the innovations and design principles of the data structure and related algorithms. This includes a metadata aware decimation algorithm which is general for collapsing simplices of any dimension. We also present an example application of this library to represent an orientable surface mesh. △ Less

Submitted 27 March, 2019; v1 submitted 3 July, 2018; originally announced July 2018.

Comments: 24 pages, 6 figures

Journal ref: ACM Trans. Math. Softw. 45, 3, Article 28 (August 2019)

Dark Matter Halo Properties vs. Local Density and Cosmic Web Location

Authors: Tze Goh, Joel Primack, Christoph T. Lee, Miguel Aragon-Calvo, Doug Hellinger, Peter Behroozi, Aldo Rodriguez-Puebla, Elliot Eckholm, Kathryn Johnston

Abstract: We study the effects of the local environmental density and the cosmic web environment (filaments, walls, and voids) on key properties of dark matter halos using the Bolshoi-Planck LCDM cosmological simulation. The z = 0 simulation is analysed into filaments, walls, and voids using the SpineWeb method and also the VIDE package of tools, both of which use the watershed transform. The key halo prope… ▽ More We study the effects of the local environmental density and the cosmic web environment (filaments, walls, and voids) on key properties of dark matter halos using the Bolshoi-Planck LCDM cosmological simulation. The z = 0 simulation is analysed into filaments, walls, and voids using the SpineWeb method and also the VIDE package of tools, both of which use the watershed transform. The key halo properties that we study are the specific mass accretion rate, spin parameter, concentration, prolateness, scale factor of the last major merger, and scale factor when the halo had half of its z = 0 mass. For all these properties, we find that there is no discernible difference between the halo properties in filaments, walls, or voids when compared at the same environmental density. As a result, we conclude that environmental density is the core attribute that affects these properties. This conclusion is in line with recent findings that properties of galaxies in redshift surveys are independent of their cosmic web environment at the same environmental density at z ~ 0. We also find that the local web environment of the Milky Way and the Andromeda galaxies near the centre of a cosmic wall does not appear to have any effect on the properties of these galaxies' dark matter halos except for their orientation, although we find that it is rather rare to have such massive halos near the centre of a relatively small cosmic wall. △ Less

Submitted 22 October, 2019; v1 submitted 13 May, 2018; originally announced May 2018.

Comments: 23 pages

Deep Learning Identifies High-z Galaxies in a Central Blue Nugget Phase in a Characteristic Mass Range

Authors: M. Huertas-Company, J. R. Primack, A. Dekel, D. C. Koo, S. Lapiner, D. Ceverino, R. C. Simons, G. F. Snyder, M. Bernardi, Z. Chen, H. Domínguez-Sánchez, Z. Chen, C. T. Lee, B. Margalef-Bentabol, D. Tuccillo

Abstract: We use machine learning to identify in color images of high-redshift galaxies an astrophysical phenomenon predicted by cosmological simulations. This phenomenon, called the blue nugget (BN) phase, is the compact star-forming phase in the central regions of many growing galaxies that follows an earlier phase of gas compaction and is followed by a central quenching phase. We train a Convolutional Ne… ▽ More We use machine learning to identify in color images of high-redshift galaxies an astrophysical phenomenon predicted by cosmological simulations. This phenomenon, called the blue nugget (BN) phase, is the compact star-forming phase in the central regions of many growing galaxies that follows an earlier phase of gas compaction and is followed by a central quenching phase. We train a Convolutional Neural Network (CNN) with mock "observed" images of simulated galaxies at three phases of evolution: pre-BN, BN and post-BN, and demonstrate that the CNN successfully retrieves the three phases in other simulated galaxies. We show that BNs are identified by the CNN within a time window of $\sim0.15$ Hubble times. When the trained CNN is applied to observed galaxies from the CANDELS survey at $z=1-3$, it successfully identifies galaxies at the three phases. We find that the observed BNs are preferentially found in galaxies at a characteristic stellar mass range, $10^{9.2-10.3} M_\odot$ at all redshifts. This is consistent with the characteristic galaxy mass for BNs as detected in the simulations, and is meaningful because it is revealed in the observations when the direct information concerning the total galaxy luminosity has been eliminated from the training set. This technique can be applied to the classification of other astrophysical phenomena for improved comparison of theory and observations in the era of large imaging surveys and cosmological simulations. △ Less

Submitted 19 April, 2018; originally announced April 2018.

Comments: Accepted for publication in ApJ

A catalog of polychromatic bulge-disk decompositions of ~ 17.600 galaxies in CANDELS

Authors: Paola Dimauro, Marc Huertas-Company, Emanuele Daddi, Pablo G. Pérez-González, Mariangela Bernardi, Guillermo Barro, Fernando Buitrago, Fernando Caro, Andrea Cattaneo, Helena Dominguez-Sánchez, Sandra M. Faber, Boris Häußler, Dale D. Kocevski, Anton M. Koekemoer, David C. Koo, Christoph T. Lee, Simona Mei, Berta Margalef-Bentabol, Joel Primack, Aldo Rodriguez-Puebla, Mara Salvato, Francesco Shankar, Diego Tuccillo

Abstract: Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge… ▽ More Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge-disk decompositions of the surface brightness profiles of ~17.600 H-band selected galaxies in the CANDELS fields (F160W<23, 0<z<2) in 4 to 7 filters covering a spectral range of 430-1600nm. This is the largest available catalog of this kind up to z = 2. By using a novel approach based on deep-learning to select the best model to fit, we manage to control systematics arising from wrong model selection and obtain less contaminated samples than previous works. We show that the derived structural properties are within $\sim10-20\%$ of random uncertainties. We then fit stellar population models to the decomposed SEDs (Spectral Energy Distribution) of bulges and disks and derive stellar masses (and stellar mass bulge-to-total ratios) as well as rest-frame colors (U,V,J) for bulges and disks separately. All data products are publicly released with this paper and through the web page https://lerma.obspm.fr/huertas/form_CANDELS and will be used for scientific analysis in forthcoming works. △ Less

Submitted 30 March, 2018; v1 submitted 27 March, 2018; originally announced March 2018.

Comments: 18 pages, 18 figures, 3 table, submitted to MNRAS after the first referee's report

Tidal Stripping and Post-Merger Relaxation of Dark Matter Halos: Causes and Consequences of Mass Loss

Authors: Christoph T. Lee, Joel R. Primack, Peter Behroozi, Aldo Rodríguez-Puebla, Doug Hellinger, Avishai Dekel

Abstract: We study the properties of distinct dark matter halos (i.e., those that are not subhalos) that have a final virial mass $M_{\mathrm{vir}}$ at $z = 0$ less than their peak mass ($M_{\mathrm{peak}}$) in the Bolshoi-Planck cosmological simulation. We identify two primary causes of halo mass loss: relaxation after a major merger and tidal stripping by a massive neighbouring halo. Major mergers initial… ▽ More We study the properties of distinct dark matter halos (i.e., those that are not subhalos) that have a final virial mass $M_{\mathrm{vir}}$ at $z = 0$ less than their peak mass ($M_{\mathrm{peak}}$) in the Bolshoi-Planck cosmological simulation. We identify two primary causes of halo mass loss: relaxation after a major merger and tidal stripping by a massive neighbouring halo. Major mergers initially boost $M_{\mathrm{vir}}$ and typically cause the final halo to become more prolate and less relaxed and to have higher spin and lower NFW concentration. As the halo relaxes, high energy material from the recent merger gradually escapes beyond the virial radius, temporarily resulting in a net negative accretion rate that reduces the halo mass by $5-15\%$ on average. Halos that experience a major merger around $z = 0.4$ typically reach a minimum mass near $z = 0$. Tidal stripping mainly occurs in dense regions, and it causes halos to become less prolate and have lower spins and higher NFW concentrations. Tidally stripped halos often lose a large fraction of their peak mass ($> 20\%$) and most never recover (or even reattain a positive accretion rate). Low mass halos can be strongly affected by both post-merger mass loss and tidal stripping, while high mass halos are predominantly influenced by post-merger mass loss and show few signs of significant tidal stripping. △ Less

Submitted 28 November, 2017; originally announced November 2017.

Textured heterogeneity in square artificial spin ice

Authors: J. C. T Lee, S. K. Mishra, V. S. Bhat, R. Streubel, B. Farmer, X. Shi, L. E. De Long, I. McNulty, P. Fischer, S. D. Kevan, S. Roy

Abstract: We report evidence of spontaneous formation of a heterogeneous network of superdomains in two-dimensional square artificial spin ice nanostructures in externally applied magnetic fields. The magnetic heterogeneity is locally disordered but has a zig-zag texture at longer length scales. Resonant coherent soft-x-ray scattering off such textures give rise to unique internal structure in Bragg peaks.… ▽ More We report evidence of spontaneous formation of a heterogeneous network of superdomains in two-dimensional square artificial spin ice nanostructures in externally applied magnetic fields. The magnetic heterogeneity is locally disordered but has a zig-zag texture at longer length scales. Resonant coherent soft-x-ray scattering off such textures give rise to unique internal structure in Bragg peaks. Our result shows that the macroscopic magnetic texture is derived from the microscopic structure of the Dirac strings. △ Less

Submitted 31 October, 2017; originally announced November 2017.

Comments: Main text: 5 pages, 4 figures. Supplemental: 13 pages, 11 figures

Journal ref: Phys. Rev. B 99, 024406 (2019)

Does the Galaxy-Halo Connection Vary with Environment?

Authors: Radu Dragomir, Aldo Rodriguez-Puebla, Joel R. Primack, Christoph T. Lee

Abstract: SubHalo Abundance Matching (SHAM) assumes that one (sub)halo property, such as mass Mvir or peak circular velocity Vpeak, determines properties of the galaxy hosted in each (sub)halo such as its luminosity or stellar mass. This assumption implies that the dependence of Galaxy Luminosity Functions (GLFs) and the Galaxy Stellar Mass Function (GSMF) on environmental density is determined by the corre… ▽ More SubHalo Abundance Matching (SHAM) assumes that one (sub)halo property, such as mass Mvir or peak circular velocity Vpeak, determines properties of the galaxy hosted in each (sub)halo such as its luminosity or stellar mass. This assumption implies that the dependence of Galaxy Luminosity Functions (GLFs) and the Galaxy Stellar Mass Function (GSMF) on environmental density is determined by the corresponding halo density dependence. In this paper, we test this by determining from an SDSS sample the observed dependence with environmental density of the ugriz GLFs and GSMF for all galaxies, and for central and satellite galaxies separately. We then show that the SHAM predictions are in remarkable agreement with these observations, even when the galaxy population is divided between central and satellite galaxies. However, we show that SHAM fails to reproduce the correct dependence between environmental density and g-r color for all galaxies and central galaxies, although it better reproduces the color dependence on environmental density of satellite galaxies. △ Less

Submitted 8 February, 2018; v1 submitted 25 October, 2017; originally announced October 2017.

Comments: 21 pages, 11 figures. Accepted for publication in MNRAS

Resonant properties of dipole skyrmions in amorphous Fe/Gd multilayers

Authors: S. A. Montoya, S. Couture, J. J. Chess, J. C. T Lee, N. Kent, M. -Y. Im, S. D. Kevan, P. Fischer, B. J. McMorran, S. Roy, V. Lomakin, E. E. Fullerton

Abstract: The dynamic response of dipole skyrmions in Fe/Gd multilayer films is investigated by ferromagnetic resonance measurements and compared to micromagnetic simulations. We detail thickness and temperature dependent studies of the observed modes as well as the effects of magnetic field history on the resonant spectra. Correlation between the modes and the magnetic phase maps constructed from real-spac… ▽ More The dynamic response of dipole skyrmions in Fe/Gd multilayer films is investigated by ferromagnetic resonance measurements and compared to micromagnetic simulations. We detail thickness and temperature dependent studies of the observed modes as well as the effects of magnetic field history on the resonant spectra. Correlation between the modes and the magnetic phase maps constructed from real-space imaging and scattering patterns allows us to conclude the resonant modes arise from local topological features such as dipole skyrmions but does not depend on the collective response of a closed packed lattice of these chiral textures. Using, micromagnetic modeling, we are able to quantitatively reproduce our experimental observations which suggests the existence of localized spin-wave modes that are dependent on the helicity of the dipole skyrmion. We identify four localized spin wave excitations for the skyrmions that are excited under either in-plane or out-of-plane r.f. fields. Lastly we show that dipole skyrmions and non-chiral bubble domains exhibit qualitatively different localized spin wave modes. △ Less

Submitted 16 May, 2017; v1 submitted 15 February, 2017; originally announced February 2017.

Comments: 38 pages, 17 figures

Journal ref: Phys. Rev. B 95, 224405 (2017)

Properties of Dark Matter Halos as a Function of Local Environment Density

Authors: Christoph T Lee, Joel R. Primack, Peter Behroozi, Aldo Rodriguez-Puebla, Doug Hellinger, Avishai Dekel

Abstract: We study how properties of discrete dark matter halos depend on halo environment, characterized by the mass density around the halos on scales from 0.5 to 16 $h^{-1}{\rm Mpc}$. We find that low mass halos (those less massive than the characteristic mass $M_{\rm C}$ of halos collapsing at a given epoch) in high-density environments have lower accretion rates, lower spins, higher concentrations, and… ▽ More We study how properties of discrete dark matter halos depend on halo environment, characterized by the mass density around the halos on scales from 0.5 to 16 $h^{-1}{\rm Mpc}$. We find that low mass halos (those less massive than the characteristic mass $M_{\rm C}$ of halos collapsing at a given epoch) in high-density environments have lower accretion rates, lower spins, higher concentrations, and rounder shapes than halos in median density environments. Halos in median and low-density environments have similar accretion rates and concentrations, but halos in low density environments have lower spins and are more elongated. Halos of a given mass in high-density regions accrete material earlier than halos of the same mass in lower-density regions. All but the most massive halos in high-density regions are losing mass (i.e., being stripped) at low redshifts, which causes artificially lowered NFW scale radii and increased concentrations. Tidal effects are also responsible for the decreasing spins of low mass halos in high density regions at low redshifts $z < 1$, by preferentially removing higher angular momentum material from halos. Halos in low-density regions have lower than average spins because they lack nearby halos whose tidal fields can spin them up. We also show that the simulation density distribution is well fit by an Extreme Value Distribution, and that the density distribution becomes broader with cosmic time. △ Less

Submitted 6 January, 2017; v1 submitted 6 October, 2016; originally announced October 2016.

Comments: 27 pages, 20 figures

Magnetically Charged Superdomain Walls In Square Artificial Spin Ice

Authors: J. C. T Lee, S. K. Mishra, V. S. Bhat, R. Streubel, B. Farmer, X. Shi, L. E. De Long, I. McNulty, P. Fischer, S. D. Kevan, S. Roy

Abstract: We report direct evidence that magnetically charged superdomain walls form spontaneously in two dimensional square artificial spin ice nanostructures in response to external magnetic fields. These extended magnetic defects were revealed by the development of internal structure, which varies as a function of applied magnetic field, within the Bragg peaks of resonant soft x-ray magnetic scattering p… ▽ More We report direct evidence that magnetically charged superdomain walls form spontaneously in two dimensional square artificial spin ice nanostructures in response to external magnetic fields. These extended magnetic defects were revealed by the development of internal structure, which varies as a function of applied magnetic field, within the Bragg peaks of resonant soft x-ray magnetic scattering patterns. Magnetically charged superdomain walls extend over tens of lattice sites and do not necessarily align with the applied field. Our results illustrate a novel approach to detect hierarchical magnetic structures within spin textures. △ Less

Submitted 25 September, 2016; originally announced September 2016.

Comments: 7 pages, 7 figures

Tailoring magnetic energies to form dipole skyrmions and skyrmion lattices

Authors: S. A. Montoya, S. Couture, J. J. Chess, J. C. T. Lee, N. Kent, D. Henze, S. K. Sinha, M. -Y. Im, S. D. Kevan, P. Fischer, B. J. McMorran, V. Lomakin, S. Roy, E. E. Fullerton

Abstract: The interesting physics and potential memory technologies resulting from topologically protected spin textures such as skyrmions, has prompted efforts to discover new material systems that can host these kind of magnetic structures. Here we use the highly tunable magnetic properties of amorphous Fe/Gd multilayer films to explore the magnetic properties that lead to dipole-stabilized skyrmions and… ▽ More The interesting physics and potential memory technologies resulting from topologically protected spin textures such as skyrmions, has prompted efforts to discover new material systems that can host these kind of magnetic structures. Here we use the highly tunable magnetic properties of amorphous Fe/Gd multilayer films to explore the magnetic properties that lead to dipole-stabilized skyrmions and skyrmion lattices that form from the competition of dipolar field and exchange energy. Using both real space imaging and reciprocal space scattering techniques we determined the range of material properties and magnetic fields where skyrmions form. Micromagnetic modeling closely matches our observation of small skyrmion features (~50 to 70nm) and suggests these class of skyrmions have a rich domain structure that is Bloch like in the center of the film and more Néel like towards each surface. Our results provide a pathway to engineer the formation and controllability of dipole skyrmion phases in a thin film geometry at different temperatures and magnetic fields. △ Less

Submitted 10 December, 2016; v1 submitted 3 August, 2016; originally announced August 2016.

Comments: 34 pages, 11 figures

Journal ref: Phys. Rev. B 95, 024415 (2017)

Synthesizing Skyrmion Molecules in Fe-Gd Thin Films

Authors: J. C. T Lee, J. J. Chess, S. A. Montoya, X. Shi, N. Tamura, S. K. Mishra, P. Fischer, B. J. McMorran, S. K. Sinha, E. E. Fullerton, S. D. Kevan, S. Roy

Abstract: We show that properly engineered amorphous Fe-Gd alloy thin films with perpendicular magnetic anisotropy exhibit room-temperature skyrmion molecules, or a pair of like-polarity, opposite-helicity skyrmions. Magnetic mirror symmetry planes present in the stripe phase, instead of chiral exchange, determine the internal skyrmion structure and the net achirality of the skyrmion phase. Our study shows… ▽ More We show that properly engineered amorphous Fe-Gd alloy thin films with perpendicular magnetic anisotropy exhibit room-temperature skyrmion molecules, or a pair of like-polarity, opposite-helicity skyrmions. Magnetic mirror symmetry planes present in the stripe phase, instead of chiral exchange, determine the internal skyrmion structure and the net achirality of the skyrmion phase. Our study shows that stripe domain engineering in amorphous alloy thin films may enable the creation of skyrmion phases with technologically desirable properties. △ Less

Submitted 30 June, 2016; v1 submitted 25 March, 2016; originally announced March 2016.

Comments: 15 pages, 6 figures. Accepted for publication in Applied Physics Letters

Known Structure, Unknown Function: An Inquiry-based Undergraduate Biochemistry Laboratory Course

Authors: Cynthia Gray, Carol W. Price, Christopher T. Lee, Alison H. Dewald, Matthew A. Cline, Charles E. McAnany, Linda Columbus, Cameron Mura

Abstract: Undergraduate biochemistry laboratory courses often do not provide students with an authentic research experience, particularly when the express purpose of the laboratory is purely instructional. However, an instructional laboratory course that is inquiry- and research-based could simultaneously impart scientific knowledge and foster a student's research expertise and confidence. We have developed… ▽ More Undergraduate biochemistry laboratory courses often do not provide students with an authentic research experience, particularly when the express purpose of the laboratory is purely instructional. However, an instructional laboratory course that is inquiry- and research-based could simultaneously impart scientific knowledge and foster a student's research expertise and confidence. We have developed a year-long undergraduate biochemistry laboratory curriculum wherein students determine, via experiment and computation, the function of a protein of known three-dimensional structure. The first half of the course is inquiry-based and modular in design; students learn general biochemical techniques while gaining preparation for research experiments in the second semester. Having learned standard biochemical methods in the first semester, students independently pursue their own (original) research projects in the second semester. This new curriculum has yielded an improvement in student performance and confidence as assessed by various metrics. To disseminate teaching resources to students and instructors alike, a freely accessible Biochemistry Laboratory Education resource is available at http://biochemlab.org. △ Less

Submitted 8 July, 2015; originally announced July 2015.

Comments: 75 pages, 8 figures, 2 tables, to appear in the journal "Biochemistry & Molecular Biology Education" (open-access license)

Anatomy-specific classification of medical images using deep convolutional nets

Authors: Holger R. Roth, Christopher T. Lee, Hoo-Chang Shin, Ari Seff, Lauren Kim, Jianhua Yao, Le Lu, Ronald M. Summers

Abstract: Automated classification of human anatomy is an important prerequisite for many computer-aided diagnosis systems. The spatial complexity and variability of anatomy throughout the human body makes classification difficult. "Deep learning" methods such as convolutional networks (ConvNets) outperform other state-of-the-art methods in image classification tasks. In this work, we present a method for o… ▽ More Automated classification of human anatomy is an important prerequisite for many computer-aided diagnosis systems. The spatial complexity and variability of anatomy throughout the human body makes classification difficult. "Deep learning" methods such as convolutional networks (ConvNets) outperform other state-of-the-art methods in image classification tasks. In this work, we present a method for organ- or body-part-specific anatomical classification of medical images acquired using computed tomography (CT) with ConvNets. We train a ConvNet, using 4,298 separate axial 2D key-images to learn 5 anatomical classes. Key-images were mined from a hospital PACS archive, using a set of 1,675 patients. We show that a data augmentation approach can help to enrich the data set and improve classification performance. Using ConvNets and data augmentation, we achieve anatomy-specific classification error of 5.9 % and area-under-the-curve (AUC) values of an average of 0.998 in testing. We demonstrate that deep learning can be used to train very reliable and accurate classifiers that could initialize further computer-aided diagnosis. △ Less

Submitted 15 April, 2015; originally announced April 2015.

Comments: Presented at: 2015 IEEE International Symposium on Biomedical Imaging, April 16-19, 2015, New York Marriott at Brooklyn Bridge, NY, USA

Journal ref: Biomedical Imaging (ISBI), 2015 IEEE 12th International Symposium on Year: 2015 Pages: 101 - 104

Coupled skyrmion sublattices in Cu2OSeO3

Authors: M. C. Langner, S. Roy, S. K. Mishra, J. C. T. Lee, X. W. Shi, M. A. Hossain, Y. -D. Chuang, S. Seki, Y. Tokura, S. D. Kevan, R. W. Schoenlein

Abstract: We report the observation of the skyrmion lattice in the chiral multiferroic insulator Cu2OSeO3 using Cu L3-edge resonant soft x-ray diffraction. We observe the unexpected existence of two distinct skyrmion sublattices that arise from inequivalent Cu sites with chemically identical coordination numbers but different magnetically active orbitals. The skyrmion sublattices are rotated with respect to… ▽ More We report the observation of the skyrmion lattice in the chiral multiferroic insulator Cu2OSeO3 using Cu L3-edge resonant soft x-ray diffraction. We observe the unexpected existence of two distinct skyrmion sublattices that arise from inequivalent Cu sites with chemically identical coordination numbers but different magnetically active orbitals. The skyrmion sublattices are rotated with respect to each other implying a long wavelength modulation of the lattice. The modulation vector could be controlled with an applied magnetic field, associating this Moire-like phase with a continuous phase transition. Our findings will open a new class of science involving manipulation of quantum topological states. △ Less

Submitted 19 April, 2014; v1 submitted 12 February, 2014; originally announced February 2014.

A supersonic turbulence origin of Larson's laws

Authors: Alexei G. Kritsuk, Christoph T. Lee, Michael L. Norman

Abstract: We revisit the origin of Larson's scaling laws describing the structure and kinematics of molecular clouds. Our analysis is based on recent observational measurements and data from a suite of six simulations of the interstellar medium, including effects of self-gravity, turbulence, magnetic field, and multiphase thermodynamics. Simulations of isothermal supersonic turbulence reproduce observed slo… ▽ More We revisit the origin of Larson's scaling laws describing the structure and kinematics of molecular clouds. Our analysis is based on recent observational measurements and data from a suite of six simulations of the interstellar medium, including effects of self-gravity, turbulence, magnetic field, and multiphase thermodynamics. Simulations of isothermal supersonic turbulence reproduce observed slopes in linewidth-size and mass-size relations. Whether or not self-gravity is included, the linewidth-size relation remains the same. The mass-size relation, instead, substantially flattens below the sonic scale, as prestellar cores start to form. Our multiphase models with magnetic field and domain size 200 pc reproduce both scaling and normalization of the first Larson law. The simulations support a turbulent interpretation of Larson's relations. This interpretation implies that: (i) the slopes of linewidth-size and mass-size correlations are determined by the inertial cascade; (ii) none of the three Larson laws is fundamental; (iii) instead, if one is known, the other two follow from scale invariance of the kinetic energy transfer rate. It does not imply that gravity is dynamically unimportant. The self-similarity of structure established by the turbulence breaks in star-forming clouds due to the development of gravitational instability in the vicinity of the sonic scale. The instability leads to the formation of prestellar cores with the characteristic mass set by the sonic scale. The high-end slope of the core mass function predicted by the scaling relations is consistent with the Salpeter power-law index. △ Less

Submitted 9 October, 2013; v1 submitted 23 September, 2013; originally announced September 2013.

Comments: 17 pages, 8 figures, MNRAS, in press

Journal ref: MNRAS, 436, 3247-3261, 2013

Emergence of charge density wave domain walls above the superconducting dome in TiSe2

Authors: Y. I. Joe, X. M. Chen, P. Ghaemi, K. D. Finkelstein, G. A. de la Peña, Y. Gan, J. C. T. Lee, S. Yuan, J. Geck, G. J. MacDougall, T. C. Chiang, S. L. Cooper, E. Fradkin, P. Abbamonte

Abstract: Superconductivity (SC) in so-called "unconventional superconductors" is nearly always found in the vicinity of another ordered state, such as antiferromagnetism, charge density wave (CDW), or stripe order. This suggests a fundamental connection between SC and fluctuations in some other order parameter. To better understand this connection, we used high-pressure x-ray scattering to directly study t… ▽ More Superconductivity (SC) in so-called "unconventional superconductors" is nearly always found in the vicinity of another ordered state, such as antiferromagnetism, charge density wave (CDW), or stripe order. This suggests a fundamental connection between SC and fluctuations in some other order parameter. To better understand this connection, we used high-pressure x-ray scattering to directly study the CDW order in the layered dichalcogenide TiSe2, which was previously shown to exhibit SC when the CDW is suppressed by pressure [1] or intercalation of Cu atoms [2]. We succeeded in suppressing the CDW fully to zero temperature, establishing for the first time the existence of a quantum critical point (QCP) at Pc = 5.1 +/- 0.2 GPa, which is more than 1 GPa beyond the end of the SC region. Unexpectedly, at P = 3 GPa we observed a reentrant, weakly first order, incommensurate phase, indicating the presence of a Lifshitz tricritical point somewhere above the superconducting dome. Our study suggests that SC in TiSe2 may not be connected to the QCP itself, but to the formation of CDW domain walls. △ Less

Submitted 16 September, 2013; originally announced September 2013.

Comments: 10 pages, 3 figures

Distinct oxygen hole doping in different layers of $\rm Sr_{2}CuO_{4-δ}-La_{2}CuO_4$ superlattices

Authors: S. Smadici, J. C. T. Lee, P. Abbamonte, A. Rusydi, G. Logvenov, I. Bozovic

Abstract: X-ray absorption in $\rm Sr_{2}CuO_{4-δ}-La_{2}CuO_4$ (SCO-LCO) superlattices shows a variable occupation with doping of a hole state different from holes doped for $x \lesssim x_{optimal}$ in bulk $\rm La_{2-x}Sr_{x}CuO_4$ and suggests that this hole state is on apical oxygen atoms and polarized in the $a-b$ plane. Considering the surface reflectivity gives a good qualitative description of the l… ▽ More X-ray absorption in $\rm Sr_{2}CuO_{4-δ}-La_{2}CuO_4$ (SCO-LCO) superlattices shows a variable occupation with doping of a hole state different from holes doped for $x \lesssim x_{optimal}$ in bulk $\rm La_{2-x}Sr_{x}CuO_4$ and suggests that this hole state is on apical oxygen atoms and polarized in the $a-b$ plane. Considering the surface reflectivity gives a good qualitative description of the line shapes of resonant soft X-ray scattering. The interference between superlattice and surface reflections was used to distinguish between scatterers in the SCO and the LCO layers, with the two hole states maximized in different layers of the superlattice. △ Less

Submitted 10 November, 2011; originally announced November 2011.

Resonant soft x-ray scattering from La(1-x)Sr(x)MnO(3) quantum wire arrays

Authors: X. M. Chen, E. M. Spanton, S. Wang, J. C. T. Lee, S. Smadici, X. Zhai, T. Naibert, J. N. Eckstein, A. Bhattacharya, T. Santos, R. Budakian, P. Abbamonte

Abstract: We describe a strategy for using resonant soft x-ray scattering (RSXS) to study the electronic structure of transition metal oxide quantum wires. Using electron beam lithography and ion milling, we have produced periodic, patterned arrays of colossal magnetoresistance (CMR) phase La(1-x)Sr(x)MnO(3) consisting of ~ 5000 wires, each of which is 80 nm in width. The scattered intensity exhibits a seri… ▽ More We describe a strategy for using resonant soft x-ray scattering (RSXS) to study the electronic structure of transition metal oxide quantum wires. Using electron beam lithography and ion milling, we have produced periodic, patterned arrays of colossal magnetoresistance (CMR) phase La(1-x)Sr(x)MnO(3) consisting of ~ 5000 wires, each of which is 80 nm in width. The scattered intensity exhibits a series of peaks that can be interpreted as Bragg reflections from the periodic structure or, equivalently, diffraction orders from the grating-like structure. RSXS measurements at the Mn L(2,3) edge, which has a large magnetic cross section, show clear evidence for a magnetic superstructure with a commensurate period of five wires, which we interpret as commensurately modulated antiferromagnetism. This superstructure, which is accompanied by non-trivial reorganization of the magnetization within each wire, likely results from classical dipole interactions among the wires. We introduce a simple, exactly soluble, analytic model of the scattering that captures, semi-quantitatively, the primary features in the RSXS data; this model will act as a foundation for forthcoming, detailed studies of the magnetic structure in these systems. △ Less

Submitted 4 September, 2011; originally announced September 2011.

Comments: 8 pages, 10 figures

Graded Orbital Occupation near Interfaces in a La2NiO4 - La2CuO4 Superlattice

Authors: S. Smadici, J. C. T. Lee, J. Morales, G. Logvenov, O. Pelleg, I. Bozovic, Y. Zhu, P. Abbamonte

Abstract: X-ray absorption spectroscopy and resonant soft x-ray reflectivity show a non-uniform distribution of oxygen holes in a La2NiO4 - La2CuO4 (LNO-LCO) superlattice, with excess holes concentrated in the LNO layers. Weak ferromagnetism with Tc = 160 K suggests a coordinated tilting of NiO6 octahedra, similar to that of bulk LNO. Ni d3z2-r2 orbitals within the LNO layers have a spatially variable occup… ▽ More X-ray absorption spectroscopy and resonant soft x-ray reflectivity show a non-uniform distribution of oxygen holes in a La2NiO4 - La2CuO4 (LNO-LCO) superlattice, with excess holes concentrated in the LNO layers. Weak ferromagnetism with Tc = 160 K suggests a coordinated tilting of NiO6 octahedra, similar to that of bulk LNO. Ni d3z2-r2 orbitals within the LNO layers have a spatially variable occupation. This variation of the Ni valence near LNO-LCO interfaces is observed with resonant soft x-ray reflectivity at the Ni L edge, at a reflection suppressed by the symmetry of the structure, and is possible through graded doping with holes, due to oxygen interstitials taken up preferentially by inner LNO layers. Since the density of oxygen atoms in the structure can be smoothly varied with standard procedures, this orbital occupation, robust up to at least 280 K, is tunable. △ Less

Submitted 11 November, 2011; v1 submitted 3 June, 2011; originally announced June 2011.

Comments: 11 pages, 8 figures

Journal ref: Physical Review B 84, 155411 (2011)

Atomic-scale patterning of hydrogen terminated Ge(001) by scanning tunneling microscopy

Authors: G. Scappucci, G. Capellini, W. C. T. Lee, M. Y. Simmons

Abstract: In this paper we demonstrate atomic-scale lithography on hydrogen terminated Ge(001. The lithographic patterns were obtained by selectively desorbing hydrogen atoms from a H resist layer adsorbed on a clean, atomically flat Ge(001) surface with a scanning tunneling microscope tip operating in ultra-high vacuum. The influence of the tip-to-sample bias on the lithographic process have been investi… ▽ More In this paper we demonstrate atomic-scale lithography on hydrogen terminated Ge(001. The lithographic patterns were obtained by selectively desorbing hydrogen atoms from a H resist layer adsorbed on a clean, atomically flat Ge(001) surface with a scanning tunneling microscope tip operating in ultra-high vacuum. The influence of the tip-to-sample bias on the lithographic process have been investigated. Lithographic patterns with feature-sizes from 200 nm to 1.8 nm have been achieved by varying the tip-to-sample bias. These results open up the possibility of a scanning-probe lithography approach to the fabrication of future atomic-scale devices in germanium. △ Less

Submitted 3 December, 2009; originally announced December 2009.

Journal ref: Nanotechnology 20, 295302 (2009)

Ultra-dense phosphorus in germanium delta-doped layers

Authors: G. Scappucci, G. Capellini, W. C. T. Lee, M. Y. Simmons

Abstract: Phosphorus (P) in germanium (Ge) delta-doped layers are fabricated in ultra-high vacuum by adsorption of phosphine molecules onto an atomically flat clean Ge(001) surface followed by thermal incorporation of P into the lattice and epitaxial Ge overgrowth by molecular beam epitaxy. Structural and electrical characterizations show that P atoms are confined, with minimal diffusion, into an ultra-na… ▽ More Phosphorus (P) in germanium (Ge) delta-doped layers are fabricated in ultra-high vacuum by adsorption of phosphine molecules onto an atomically flat clean Ge(001) surface followed by thermal incorporation of P into the lattice and epitaxial Ge overgrowth by molecular beam epitaxy. Structural and electrical characterizations show that P atoms are confined, with minimal diffusion, into an ultra-narrow 2-nm-wide layer with an electrically-active sheet carrier concentration of 4x10^13 cm-2 at 4.2 K. These results open up the possibility of ultra-narrow source/drain regions with unprecedented carrier densities for Ge n-channel field effect transistors. △ Less

Submitted 3 December, 2009; originally announced December 2009.

Journal ref: Appl. Phys. Lett. 94, 162106 (2009)

Two-step stabilization of orbital order and the dynamical frustration of spin in the model charge-transfer insulator KCuF3

Authors: James C. T. Lee, Shi Yuan, Siddhartha Lal, Young Il Joe, Yu Gan, Serban Smadici, Ken Finkelstein, Yejun Feng, Andrivo Rusydi, Paul M. Goldbart, S. Lance Cooper, Peter Abbamonte

Abstract: We report a combined experimental and theoretical study of KCuF3, which offers - because of this material's relatively simple lattice structure and valence configuration (d9, i.e., one hole in the d-shell) - a particularly clear view of the essential role of the orbital degree of freedom in governing the dynamical coupling between the spin and lattice degrees of freedom. We present Raman and x-r… ▽ More We report a combined experimental and theoretical study of KCuF3, which offers - because of this material's relatively simple lattice structure and valence configuration (d9, i.e., one hole in the d-shell) - a particularly clear view of the essential role of the orbital degree of freedom in governing the dynamical coupling between the spin and lattice degrees of freedom. We present Raman and x-ray scattering evidence that the phase behaviour of KCuF3 is dominated above the Neel temperature (T_N = 40 K) by coupled orbital/lattice fluctuations that are likely associated with rotations of the CuF6 octahedra, and we show that these orbital fluctuations are interrupted by a static structural distortion that occurs just above T_N. A detailed model of the orbital and magnetic phases of KCuF3 reveals that these orbital fluctuations - and the related frustration of in-plane spin-order-are associated with the presence of nearly degenerate low-energy spin-orbital states that are highly susceptible to thermal fluctuations over a wide range of temperatures. A striking implication of these results is that the ground state of KCuF3 at ambient pressure lies near a quantum critical point associated with an orbital/spin liquid phase that is obscured by emergent Neel ordering of the spins; this exotic liquid phase might be accessible via pressure studies. △ Less

Submitted 3 November, 2009; originally announced November 2009.

Comments: 13 pages, 3 figures

Journal ref: Nature Physics, Vol. 8, 63 (2012)

Superconducting transition at 38 K in insulating-overdoped La2CuO4-La1.64Sr0.36CuO4 superlattices: Evidence for interface electronic redistribution from resonant soft x-ray scattering

Authors: S. Smadici, J. C. T. Lee, S. Wang, P. Abbamonte, G. Logvenov, A. Gozar, C. D. Cavellin, I. Bozovic

Abstract: We use resonant soft x-ray scattering (RSXS) to quantify the hole distribution in a superlattice of insulating La2CuO4 (LCO) and overdoped La_{2-x}Sr_{x}CuO4 (LSCO). Despite its non-superconducting constituents, this structure is superconducting with T_{c}=38 K. We found that the conducting holes redistribute electronically from LSCO to the LCO layers. The LCO layers were found to be optimally d… ▽ More We use resonant soft x-ray scattering (RSXS) to quantify the hole distribution in a superlattice of insulating La2CuO4 (LCO) and overdoped La_{2-x}Sr_{x}CuO4 (LSCO). Despite its non-superconducting constituents, this structure is superconducting with T_{c}=38 K. We found that the conducting holes redistribute electronically from LSCO to the LCO layers. The LCO layers were found to be optimally doped, suggesting they are the main drivers of superconductivity. Our results demonstrate the utility of RSXS for separating electronic from structural effects at oxide interfaces. △ Less

Submitted 13 February, 2009; v1 submitted 20 May, 2008; originally announced May 2008.

Comments: 4 pages, 3 figures, to appear in Physical Review Letters