The Molino suite contains 75,000 galaxy mock catalogs designed to quantify the information content of any cosmological observable for a redshift-space galaxy sample. They are constructed from the Quijote N-body simulations (Villaescusa-Navarro et al. 2020) using the standard Zheng et al. (2007) Halo Occupation Distribution (HOD) model. The fiducial HOD parameters are based on the SDSS high luminosity samples. The suite contains 15,000 mocks at the fiducial cosmology and HOD parameters for covariance matrix estimation. It also includes (500 N-body realizations) x (5 HOD realizations)=2,500 mocks at 24 other parameter values to estimate the derivative of the observable with respect to six cosmological parameters (Omega_m, Omega_b, h, n_s, sigma_8, and M_nu) and five HOD parameters (logMmin, sigma_logM, log M_0, alpha, and log M_1). Using the covariance matrix and derivatives calculated from Molino, one can derive Fisher matrix forecasts on the cosmological parameters marginalized over HOD parameters.
Recent advances in experimental techniques have allowed the simultaneous recordings of
populations of hundreds of neurons, fostering a debate about the nature of the collective
structure of population neural activity. Much of this debate has focused on the
empirical findings of a phase transition in the parameter space of maximum entropy
models describing the measured neural probability distributions, interpreting this phase
transition to indicate a critical tuning of the neural code. Here, we instead focus on the
possibility that this is a first-order phase transition which provides evidence that the
real neural population is in a `structured', collective state. We show that this collective
state is robust to changes in stimulus ensemble and adaptive state. We find that the
pattern of pairwise correlations between neurons has a strength that is well within the
strongly correlated regime and does not require fine tuning, suggesting that this state is
generic for populations of 100+ neurons. We find a clear correspondence between the
emergence of a phase transition, and the emergence of attractor-like structure in the
inferred energy landscape. A collective state in the neural population, in which neural
activity patterns naturally form clusters, provides a consistent interpretation for our
Verdoolaege, G.; Kaye, S.M.; Angioni, C.; Kardaunn, O.W.J.F.; Maslov, M.; Romanelli, M.; Ryter, F.; Thomsen, K.
The multi-machine ITPA Global H-mode Confinement Database has been upgraded with new data from JET with the ITER-like wall and ASDEX Upgrade with the full tungsten wall. This paper describes the new database and presents results of regression analysis to estimate the global energy confinement scaling in H-mode plasmas using a standard power law. Various subsets of the database are considered, focusing on type of wall and divertor materials, confinement regime (all H-modes, ELMy H or ELM-free) and ITER-like constraints. Apart from ordinary least squares, two other, robust regression techniques are applied, which take into account uncertainty on all variables. Regression on data from individual devices shows that, generally, the confinement dependence on density and the power degradation are weakest in the fully metallic devices. Using the multi-machine scalings, predictions are made of the confinement time in a standard ELMy H-mode scenario in ITER. The uncertainty on the scaling parameters is discussed with a view to practically useful error bars on the parameters and predictions. One of the derived scalings for ELMy H-modes on an ITER-like subset is studied in particular and compared to the IPB98(y,2) confinement scaling in engineering and dimensionless form. Transformation of this new scaling from engineering variables to dimensionless quantities is shown to result in large error bars on the dimensionless scaling. Regression analysis in the space of dimensionless variables is therefore proposed as an alternative, yielding acceptable estimates for the dimensionless scaling. The new scaling, which is dimensionally correct within the uncertainties, suggests that some dependencies of confinement in the multi- machine database can be reconciled with parameter scans in individual devices. This includes vanishingly small dependence of confinement on line-averaged density and normalized plasma pressure (β), as well as a noticeable, positive dependence on effective atomic mass and plasma triangularity. Extrapolation of this scaling to ITER yields a somewhat lower confinement time compared to the IPB98(y, 2) prediction, possibly related to the considerably weaker dependence on major radius in the new scaling (slightly above linear). Further studies are needed to compare more flexible regression models with the power law used here. In addition, data from more devices concerning possible ‘hidden variables’ could help to determine their influence on confinement, while adding data in sparsely populated areas of the parameter space may contribute to further disentangling some of the global confinement dependencies in tokamak plasmas.
Webb, Michael; Jacobs, William; An, Yaxin; Oliver, Wesley
This distribution compiles thermodynamic and (where available) dynamic properties of short protein sequences as obtained from coarse-grained molecular dynamics simulations. The dataset features 2114 protein sequences with sequence lengths ranging from N=20 up to N=50 amino acids. The simulation and analysis of these sequences is described in "Active learning of the thermodynamics--dynamics tradeoff in protein condensates'' by Yaxin An, Michael A. Webb*, and William M. Jacobs* (https://doi.org/10.48550/arXiv.2306.03696). Of the 2114 protein sequences, 80 are homomeric polypeptides (replicating a single amino acid for N = 20, 30, 40, and 50), 1266 are sourced from version 9.0 of the DisProt database, and the remaining 768 sequences are novel sequences generated during an active learning campaign described in the aforementioned manuscript. The simulations were performed using the LAMMPS molecular dynamics engine. The interactions used for simulation are obtained from R. M. Regy , J. Thompson , Y. C. Kim and J. Mittal , Improved coarse-grained model for studying sequence dependent phase separation of disordered proteins, Protein Sci., 2021, 1371 —1379. Properties included in this distribution include second virial coefficients, pressure-density data, expectation for phase behavior at 300 K, estimated condensed-phase densities at 300 K (if exist), and condensed-phase self-diffusion coefficients at 300 K (if exist).
In our study, we compare the three dimensional (3D) morphologic characteristics of Earth's first reef-building animals (archaeocyath sponges) with those of modern, photosynthetic corals. Within this repository are the 3D image data products for both groups of animals. The archaeocyath images were produced through serial grinding and imaging with the Grinding, Imaging, and Reconstruction Instrument at Princeton University. The images in this repository are the downsampled data products used in our study, and the full resolution (>2TB) image stacks are available upon request from the author. For the coral image data, the computed tomography (CT) images of all samples are included at full resolution. Also included in this repository are the manual and automated outline coordinates of the archaeocyath and coral branches, which can be directly used for morphological study.
Schwartz, Jacob A.; Nelson, A. O.; Kolemen, Egemen
Shaping a tokamak plasma to have a negative triangularity may allow operation in an ELM-free L-mode regime and with a larger strike-point radius, ameliorating divertor power-handling requirements. However, the shaping has a potential drawback in the form of a lower no-wall ideal beta limit, found using the MHD codes CHEASE and DCON. Using the new fusion systems code FAROES, we construct a steady-state DEMO2 reactor model. This model is essentially zero-dimensional and neglects variations in physical mechanisms like turbulence, confinement, and radiative power limits, which could have a substantial impact on the conclusions deduced herein. Keeping its shape otherwise constant, we alter the triangularity and compute the effects on the levelized cost of energy (LCOE). If the tokamak is limited to a fixed B field, then unless other means to increase performance (such as reduced turbulence, improved current drive efficiency or higher density operation) can be leveraged, a negative-triangularity reactor is strongly disfavored in the model due to lower \beta_N limits at negative triangularity, which leads to tripling of the LCOE. However, if the reactor is constrained by divertor heat fluxes and not by magnet engineering, then a negative-triangularity reactor with higher B0 could be favorable: we find a class of solutions at negative triangularity with lower peak heat flux and lower LCOE than those of the equivalent positive triangularity reactors.