@article{14425, abstract = {Water adsorption and dissociation processes on pristine low-index TiO2 interfaces are important but poorly understood outside the well-studied anatase (101) and rutile (110). To understand these, we construct three sets of machine learning potentials that are simultaneously applicable to various TiO2 surfaces, based on three density-functional-theory approximations. Here we show the water dissociation free energies on seven pristine TiO2 surfaces, and predict that anatase (100), anatase (110), rutile (001), and rutile (011) favor water dissociation, anatase (101) and rutile (100) have mostly molecular adsorption, while the simulations of rutile (110) sensitively depend on the slab thickness and molecular adsorption is preferred with thick slabs. Moreover, using an automated algorithm, we reveal that these surfaces follow different types of atomistic mechanisms for proton transfer and water dissociation: one-step, two-step, or both. These mechanisms can be rationalized based on the arrangements of water molecules on the different surfaces. Our finding thus demonstrates that the different pristine TiO2 surfaces react with water in distinct ways, and cannot be represented using just the low-energy anatase (101) and rutile (110) surfaces.}, author = {Zeng, Zezhu and Wodaczek, Felix and Liu, Keyang and Stein, Frederick and Hutter, Jürg and Chen, Ji and Cheng, Bingqing}, issn = {2041-1723}, journal = {Nature Communications}, publisher = {Springer Nature}, title = {{Mechanistic insight on water dissociation on pristine low-index TiO2 surfaces from machine learning molecular dynamics simulations}}, doi = {10.1038/s41467-023-41865-8}, volume = {14}, year = {2023}, } @article{14603, abstract = {Computing the solubility of crystals in a solvent using atomistic simulations is notoriously challenging due to the complexities and convergence issues associated with free-energy methods, as well as the slow equilibration in direct-coexistence simulations. This paper introduces a molecular-dynamics workflow that simplifies and robustly computes the solubility of molecular or ionic crystals. This method is considerably more straightforward than the state-of-the-art, as we have streamlined and optimised each step of the process. Specifically, we calculate the chemical potential of the crystal using the gas-phase molecule as a reference state, and employ the S0 method to determine the concentration dependence of the chemical potential of the solute. We use this workflow to predict the solubilities of sodium chloride in water, urea polymorphs in water, and paracetamol polymorphs in both water and ethanol. Our findings indicate that the predicted solubility is sensitive to the chosen potential energy surface. Furthermore, we note that the harmonic approximation often fails for both molecular crystals and gas molecules at or above room temperature, and that the assumption of an ideal solution becomes less valid for highly soluble substances.}, author = {Reinhardt, Aleks and Chew, Pin Yu and Cheng, Bingqing}, issn = {1089-7690}, journal = {Journal of Chemical Physics}, number = {18}, publisher = {AIP Publishing}, title = {{A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals}}, doi = {10.1063/5.0173341}, volume = {159}, year = {2023}, } @article{13231, abstract = {We study ab initio approaches for calculating x-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula that expresses the inelastic contribution in terms of the dielectric function. We study the electronic dynamic structure factor computed from the Mermin dielectric function using an ab initio electron-ion collision frequency in comparison to computations using a linear-response time-dependent density functional theory (LR-TDDFT) framework for hydrogen and beryllium and investigate the dispersion of free-free and bound-free contributions to the scattering signal. A separate treatment of these contributions, where only the free-free part follows the Mermin dispersion, shows good agreement with LR-TDDFT results for ambient-density beryllium, but breaks down for highly compressed matter where the bound states become pressure ionized. LR-TDDFT is used to reanalyze x-ray Thomson scattering experiments on beryllium demonstrating strong deviations from the plasma conditions inferred with traditional analytic models at small scattering angles.}, author = {Schörner, Maximilian and Bethkenhagen, Mandy and Döppner, Tilo and Kraus, Dominik and Fletcher, Luke B. and Glenzer, Siegfried H. and Redmer, Ronald}, issn = {2470-0053}, journal = {Physical Review E}, number = {6}, publisher = {American Physical Society}, title = {{X-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula}}, doi = {10.1103/PhysRevE.107.065207}, volume = {107}, year = {2023}, } @misc{14619, abstract = {Data underlying the publication "A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals" (DOI https://doi.org/10.1063/5.0173341).}, author = {Cheng, Bingqing}, publisher = {Zenodo}, title = {{BingqingCheng/solubility: V1.0}}, doi = {10.5281/ZENODO.8398094}, year = {2023}, } @article{13216, abstract = {Physical catalysts often have multiple sites where reactions can take place. One prominent example is single-atom alloys, where the reactive dopant atoms can preferentially locate in the bulk or at different sites on the surface of the nanoparticle. However, ab initio modeling of catalysts usually only considers one site of the catalyst, neglecting the effects of multiple sites. Here, nanoparticles of copper doped with single-atom rhodium or palladium are modeled for the dehydrogenation of propane. Single-atom alloy nanoparticles are simulated at 400–600 K, using machine learning potentials trained on density functional theory calculations, and then the occupation of different single-atom active sites is identified using a similarity kernel. Further, the turnover frequency for all possible sites is calculated for propane dehydrogenation to propene through microkinetic modeling using density functional theory calculations. The total turnover frequencies of the whole nanoparticle are then described from both the population and the individual turnover frequency of each site. Under operating conditions, rhodium as a dopant is found to almost exclusively occupy (111) surface sites while palladium as a dopant occupies a greater variety of facets. Undercoordinated dopant surface sites are found to tend to be more reactive for propane dehydrogenation compared to the (111) surface. It is found that considering the dynamics of the single-atom alloy nanoparticle has a profound effect on the calculated catalytic activity of single-atom alloys by several orders of magnitude.}, author = {Bunting, Rhys and Wodaczek, Felix and Torabi, Tina and Cheng, Bingqing}, issn = {1520-5126}, journal = {Journal of the American Chemical Society}, keywords = {Colloid and Surface Chemistry, Biochemistry, General Chemistry, Catalysis}, number = {27}, pages = {14894--14902}, publisher = {American Chemical Society}, title = {{Reactivity of single-atom alloy nanoparticles: Modeling the dehydrogenation of propane}}, doi = {10.1021/jacs.3c04030}, volume = {145}, year = {2023}, } @article{22505, abstract = {High Mountain Asia (HMA) is among the most vulnerable water towers globally and yet future projections of water availability in and from its high-mountain catchments remain uncertain, as their hydrologic response to ongoing environmental changes is complex. Mechanistic modeling approaches incorporating cryospheric, hydrological, and vegetation processes in high spatial, temporal, and physical detail have never been applied for high-elevation catchments of HMA. We use a land surface model at high spatial and temporal resolution (100 m and hourly) to simulate the coupled dynamics of energy, water, and vegetation for the 350 km2 Langtang catchment (Nepal). We compare our model outputs for one hydrological year against a large set of observations to gain insight into the partitioning of the water balance at the subseasonal scale and across elevation bands. During the simulated hydrological year, we find that evapotranspiration is a key component of the total water balance, as it causes about the equivalent of 20% of all the available precipitation or 154% of the water production from glacier melt in the basin to return directly to the atmosphere. The depletion of the cryospheric water budget is dominated by snow melt, but at high elevations is primarily dictated by snow and ice sublimation. Snow sublimation is the dominant vapor flux (49%) at the catchment scale, accounting for the equivalent of 11% of snowfall, 17% of snowmelt, and 75% of ice melt, respectively. We conclude that simulations should consider sublimation and other evaporative fluxes explicitly, as otherwise water balance estimates can be ill-quantified.}, author = {Buri, Pascal and Fatichi, Simone and Shaw, Thomas E. and Miles, Evan S. and McCarthy, Michael J. and Fyffe, Catriona L. and Fugger, Stefan and Ren, Shaoting and Kneib, Marin and Jouberton, Achille and Steiner, Jakob and Fujita, Koji and Pellicciotti, Francesca}, issn = {1944-7973}, journal = {Water Resources Research}, number = {10}, publisher = {American Geophysical Union}, title = {{Land surface modeling in the himalayas: on the importance of evaporative fluxes for the water balance of a high‐elevation catchment}}, doi = {10.1029/2022wr033841}, volume = {59}, year = {2023}, } @article{13988, abstract = {Most permissionless blockchains inherently suffer from throughput limitations. Layer-2 systems, such as side-chains or Rollups, have been proposed as a possible strategy to overcome this limitation. Layer-2 systems interact with the main-chain in two ways. First, users can move funds from/to the main-chain to/from the layer-2. Second, layer-2 systems periodically synchronize with the main-chain to keep some form of log of their activity on the main-chain - this log is key for security. Due to this interaction with the main-chain, which is necessary and recurrent, layer-2 systems impose some load on the main-chain. The impact of such load on the main-chain has been, so far, poorly understood. In addition to that, layer-2 approaches typically sacrifice decentralization and security in favor of higher throughput. This paper presents an experimental study that analyzes the current state of Ethereum layer-2 projects. Our goal is to assess the load they impose on Ethereum and to understand their scalability potential in the long-run. Our analysis shows that the impact of any given layer-2 on the main-chain is the result of both technical aspects (how state is logged on the main-chain) and user behavior (how often users decide to transfer funds between the layer-2 and the main-chain). Based on our observations, we infer that without efficient mechanisms that allow users to transfer funds in a secure and fast manner directly from one layer-2 project to another, current layer-2 systems will not be able to scale Ethereum effectively, regardless of their technical solutions. Furthermore, from our results, we conclude that the layer-2 systems that offer similar security guarantees as Ethereum have limited scalability potential, while approaches that offer better performance, sacrifice security and lead to an increase in centralization which runs against the end-goals of permissionless blockchains.}, author = {Neiheiser, Ray and Inacio, Gustavo and Rech, Luciana and Montez, Carlos and Matos, Miguel and Rodrigues, Luis}, issn = {2169-3536}, journal = {IEEE Access}, keywords = {General Engineering, General Materials Science, General Computer Science, Electrical and Electronic Engineering}, pages = {8651--8662}, publisher = {IEEE}, title = {{Practical limitations of Ethereum’s layer-2}}, doi = {10.1109/access.2023.3237897}, volume = {11}, year = {2023}, } @phdthesis{13074, abstract = {Deep learning has become an integral part of a large number of important applications, and many of the recent breakthroughs have been enabled by the ability to train very large models, capable to capture complex patterns and relationships from the data. At the same time, the massive sizes of modern deep learning models have made their deployment to smaller devices more challenging; this is particularly important, as in many applications the users rely on accurate deep learning predictions, but they only have access to devices with limited memory and compute power. One solution to this problem is to prune neural networks, by setting as many of their parameters as possible to zero, to obtain accurate sparse models with lower memory footprint. Despite the great research progress in obtaining sparse models that preserve accuracy, while satisfying memory and computational constraints, there are still many challenges associated with efficiently training sparse models, as well as understanding their generalization properties. The focus of this thesis is to investigate how the training process of sparse models can be made more efficient, and to understand the differences between sparse and dense models in terms of how well they can generalize to changes in the data distribution. We first study a method for co-training sparse and dense models, at a lower cost compared to regular training. With our method we can obtain very accurate sparse networks, and dense models that can recover the baseline accuracy. Furthermore, we are able to more easily analyze the differences, at prediction level, between the sparse-dense model pairs. Next, we investigate the generalization properties of sparse neural networks in more detail, by studying how well different sparse models trained on a larger task can adapt to smaller, more specialized tasks, in a transfer learning scenario. Our analysis across multiple pruning methods and sparsity levels reveals that sparse models provide features that can transfer similarly to or better than the dense baseline. However, the choice of the pruning method plays an important role, and can influence the results when the features are fixed (linear finetuning), or when they are allowed to adapt to the new task (full finetuning). Using sparse models with fixed masks for finetuning on new tasks has an important practical advantage, as it enables training neural networks on smaller devices. However, one drawback of current pruning methods is that the entire training cycle has to be repeated to obtain the initial sparse model, for every sparsity target; in consequence, the entire training process is costly and also multiple models need to be stored. In the last part of the thesis we propose a method that can train accurate dense models that are compressible in a single step, to multiple sparsity levels, without additional finetuning. Our method results in sparse models that can be competitive with existing pruning methods, and which can also successfully generalize to new tasks.}, author = {Peste, Elena-Alexandra}, issn = {2663-337X}, pages = {147}, publisher = {Institute of Science and Technology Austria}, title = {{Efficiency and generalization of sparse neural networks}}, doi = {10.15479/at:ista:13074}, year = {2023}, } @inproceedings{12548, abstract = {The limited exchange between human communities is a key factor in preventing the spread of COVID-19. This paper introduces a digital framework that combines an integration of real mobility data at the country scale with a series of modeling techniques and visual capabilities that highlight mobility patterns before and during the pandemic. The findings not only significantly exhibit mobility trends and different degrees of similarities at regional and local levels but also provide potential insight into the emergence of a pandemic on human behavior patterns and their likely socio-economic impacts.}, author = {Forghani, Mohammad and Claramunt, Christophe and Karimipour, Farid and Heiler, Georg}, booktitle = {2022 IEEE International Conference on Data Mining Workshops}, issn = {2375-9259}, location = {Orlando, FL, United States}, publisher = {IEEE}, title = {{Visual analytics of mobility network changes observed using mobile phone data during COVID-19 pandemic}}, doi = {10.1109/icdmw58026.2022.00093}, year = {2023}, } @inproceedings{13321, abstract = {We consider the problem of reconstructing the signal and the hidden variables from observations coming from a multi-layer network with rotationally invariant weight matrices. The multi-layer structure models inference from deep generative priors, and the rotational invariance imposed on the weights generalizes the i.i.d. Gaussian assumption by allowing for a complex correlation structure, which is typical in applications. In this work, we present a new class of approximate message passing (AMP) algorithms and give a state evolution recursion which precisely characterizes their performance in the large system limit. In contrast with the existing multi-layer VAMP (ML-VAMP) approach, our proposed AMP – dubbed multilayer rotationally invariant generalized AMP (ML-RI-GAMP) – provides a natural generalization beyond Gaussian designs, in the sense that it recovers the existing Gaussian AMP as a special case. Furthermore, ML-RI-GAMP exhibits a significantly lower complexity than ML-VAMP, as the computationally intensive singular value decomposition is replaced by an estimation of the moments of the design matrices. Finally, our numerical results show that this complexity gain comes at little to no cost in the performance of the algorithm.}, author = {Xu, Yizhou and Hou, Tian Qi and Liang, Shan Suo and Mondelli, Marco}, booktitle = {2023 IEEE Information Theory Workshop}, isbn = {9798350301496}, issn = {2475-4218}, location = {Saint-Malo, France}, pages = {294--298}, publisher = {IEEE}, title = {{Approximate message passing for multi-layer estimation in rotationally invariant models}}, doi = {10.1109/ITW55543.2023.10160238}, year = {2023}, } @article{12704, abstract = {Adversarial training (i.e., training on adversarially perturbed input data) is a well-studied method for making neural networks robust to potential adversarial attacks during inference. However, the improved robustness does not come for free but rather is accompanied by a decrease in overall model accuracy and performance. Recent work has shown that, in practical robot learning applications, the effects of adversarial training do not pose a fair trade-off but inflict a net loss when measured in holistic robot performance. This work revisits the robustness-accuracy trade-off in robot learning by systematically analyzing if recent advances in robust training methods and theory in conjunction with adversarial robot learning, are capable of making adversarial training suitable for real-world robot applications. We evaluate three different robot learning tasks ranging from autonomous driving in a high-fidelity environment amenable to sim-to-real deployment to mobile robot navigation and gesture recognition. Our results demonstrate that, while these techniques make incremental improvements on the trade-off on a relative scale, the negative impact on the nominal accuracy caused by adversarial training still outweighs the improved robustness by an order of magnitude. We conclude that although progress is happening, further advances in robust learning methods are necessary before they can benefit robot learning tasks in practice.}, author = {Lechner, Mathias and Amini, Alexander and Rus, Daniela and Henzinger, Thomas A}, issn = {2377-3766}, journal = {IEEE Robotics and Automation Letters}, number = {3}, pages = {1595--1602}, publisher = {IEEE}, title = {{Revisiting the adversarial robustness-accuracy tradeoff in robot learning}}, doi = {10.1109/LRA.2023.3240930}, volume = {8}, year = {2023}, } @inproceedings{13967, abstract = {A classic solution technique for Markov decision processes (MDP) and stochastic games (SG) is value iteration (VI). Due to its good practical performance, this approximative approach is typically preferred over exact techniques, even though no practical bounds on the imprecision of the result could be given until recently. As a consequence, even the most used model checkers could return arbitrarily wrong results. Over the past decade, different works derived stopping criteria, indicating when the precision reaches the desired level, for various settings, in particular MDP with reachability, total reward, and mean payoff, and SG with reachability.In this paper, we provide the first stopping criteria for VI on SG with total reward and mean payoff, yielding the first anytime algorithms in these settings. To this end, we provide the solution in two flavours: First through a reduction to the MDP case and second directly on SG. The former is simpler and automatically utilizes any advances on MDP. The latter allows for more local computations, heading towards better practical efficiency.Our solution unifies the previously mentioned approaches for MDP and SG and their underlying ideas. To achieve this, we isolate objective-specific subroutines as well as identify objective-independent concepts. These structural concepts, while surprisingly simple, form the very essence of the unified solution.}, author = {Kretinsky, Jan and Meggendorfer, Tobias and Weininger, Maximilian}, booktitle = {38th Annual ACM/IEEE Symposium on Logic in Computer Science}, isbn = {9798350335873}, issn = {1043-6871}, location = {Boston, MA, United States}, publisher = {IEEE}, title = {{Stopping criteria for value iteration on stochastic games with quantitative objectives}}, doi = {10.1109/LICS56636.2023.10175771}, volume = {2023}, year = {2023}, } @article{14751, abstract = {We consider zero-error communication over a two-transmitter deterministic adversarial multiple access channel (MAC) governed by an adversary who has access to the transmissions of both senders (hence called omniscient ) and aims to maliciously corrupt the communication. None of the encoders, jammer and decoder is allowed to randomize using private or public randomness. This enforces a combinatorial nature of the problem. Our model covers a large family of channels studied in the literature, including all deterministic discrete memoryless noisy or noiseless MACs. In this work, given an arbitrary two-transmitter deterministic omniscient adversarial MAC, we characterize when the capacity region: 1) has nonempty interior (in particular, is two-dimensional); 2) consists of two line segments (in particular, has empty interior); 3) consists of one line segment (in particular, is one-dimensional); 4) or only contains (0,0) (in particular, is zero-dimensional). This extends a recent result by Wang et al. (201 9) from the point-to-point setting to the multiple access setting. Indeed, our converse arguments build upon their generalized Plotkin bound and involve delicate case analysis. One of the technical challenges is to take care of both “joint confusability” and “marginal confusability”. In particular, the treatment of marginal confusability does not follow from the point-to-point results by Wang et al. Our achievability results follow from random coding with expurgation.}, author = {Zhang, Yihan}, issn = {1557-9654}, journal = {IEEE Transactions on Information Theory}, keywords = {Computer Science Applications, Information Systems}, number = {7}, pages = {4093--4127}, publisher = {IEEE}, title = {{Zero-error communication over adversarial MACs}}, doi = {10.1109/tit.2023.3257239}, volume = {69}, year = {2023}, } @article{13269, abstract = {This paper is a collection of results on combinatorial properties of codes for the Z-channel . A Z-channel with error fraction τ takes as input a length- n binary codeword and injects in an adversarial manner up to n τ asymmetric errors, i.e., errors that only zero out bits but do not flip 0’s to 1’s. It is known that the largest ( L - 1)-list-decodable code for the Z-channel with error fraction τ has exponential size (in n ) if τ is less than a critical value that we call the ( L - 1)- list-decoding Plotkin point and has constant size if τ is larger than the threshold. The ( L -1)-list-decoding Plotkin point is known to be L -1/L-1 – L -L/ L-1 , which equals 1/4 for unique-decoding with L -1 = 1. In this paper, we derive various results for the size of the largest codes above and below the list-decoding Plotkin point. In particular, we show that the largest ( L -1)-list-decodable code ε-above the Plotkin point, for any given sufficiently small positive constant ε > 0, has size Θ L (ε -3/2 ) for any L - 1 ≥ 1. We also devise upper and lower bounds on the exponential size of codes below the list-decoding Plotkin point.}, author = {Polyanskii, Nikita and Zhang, Yihan}, issn = {1557-9654}, journal = {IEEE Transactions on Information Theory}, number = {10}, pages = {6340--6357}, publisher = {IEEE}, title = {{Codes for the Z-channel}}, doi = {10.1109/TIT.2023.3292219}, volume = {69}, year = {2023}, } @article{22554, abstract = {Investigating modifications in the hydrological cycle is essential to understand the impacts of climate change on ecosystems. This study assesses the change in the velocity of the water cycle over land at the global scale, whereas previous studies have mostly focused on changes in the atmospheric water cycle. The hydrological acceleration is quantified by a decrease in average residence time (RT) of water in the first meter of soil. The soil water RT is shown to be sensitive to the soil texture and seasonality of hydroclimatic variables. Despite substantial local variability, most of the RTs are in the range of 50–300 days. The global mean soil water RT declined at a rate of −2.30 and −0.36 days decade−1 (−1.6 to 1.0 days decade−1 the range of nine models) from 2001 to 2020 as measured by reanalysis and CMIP6 simulations for the historical scenario, respectively, which corresponds to −6.8 and −1.1 days °C−1 when expressed per degree of global warming over land. This acceleration is projected to continue at a rate of −1.35 days decade−1 (−3.4 to 0.0 days decade−1 the range of nine models) or −2.2 days °C−1 during the period 2015–2100 under the most extreme emission scenario: SSP 585. Changes in precipitation dominantly drive the acceleration of the terrestrial water cycle compared to changes in evapotranspiration. Rising temperatures and increasing carbon dioxide have opposite effects on the speed of the terrestrial water cycle with compensatory roles keeping RT relatively unchanged in the absence of PR trends.}, author = {Wang, Y. and Meili, N. and Fatichi, Simone}, issn = {1944-7973}, journal = {Water Resources Research}, number = {8}, publisher = {American Geophysical Union}, title = {{Evidence and controls of the acceleration of the hydrological cycle over land}}, doi = {10.1029/2022wr033970}, volume = {59}, year = {2023}, } @inproceedings{13221, abstract = {The safety-liveness dichotomy is a fundamental concept in formal languages which plays a key role in verification. Recently, this dichotomy has been lifted to quantitative properties, which are arbitrary functions from infinite words to partially-ordered domains. We look into harnessing the dichotomy for the specific classes of quantitative properties expressed by quantitative automata. These automata contain finitely many states and rational-valued transition weights, and their common value functions Inf, Sup, LimInf, LimSup, LimInfAvg, LimSupAvg, and DSum map infinite words into the totallyordered domain of real numbers. In this automata-theoretic setting, we establish a connection between quantitative safety and topological continuity and provide an alternative characterization of quantitative safety and liveness in terms of their boolean counterparts. For all common value functions, we show how the safety closure of a quantitative automaton can be constructed in PTime, and we provide PSpace-complete checks of whether a given quantitative automaton is safe or live, with the exception of LimInfAvg and LimSupAvg automata, for which the safety check is in ExpSpace. Moreover, for deterministic Sup, LimInf, and LimSup automata, we give PTime decompositions into safe and live automata. These decompositions enable the separation of techniques for safety and liveness verification for quantitative specifications.}, author = {Boker, Udi and Henzinger, Thomas A and Mazzocchi, Nicolas Adrien and Sarac, Naci E}, booktitle = {34th International Conference on Concurrency Theory}, isbn = {9783959772990}, issn = {1868-8969}, location = {Antwerp, Belgium}, publisher = {Schloss Dagstuhl - Leibniz-Zentrum für Informatik}, title = {{Safety and liveness of quantitative automata}}, doi = {10.4230/LIPIcs.CONCUR.2023.17}, volume = {279}, year = {2023}, } @article{13965, abstract = {Many modes and mechanisms of epigenetic inheritance have been elucidated in eukaryotes. Most of them are relatively short-term, generally not exceeding one or a few organismal generations. However, emerging evidence indicates that one mechanism, cytosine DNA methylation, can mediate epigenetic inheritance over much longer timescales, which are mostly or completely inaccessible in the laboratory. Here we discuss the evidence for, and mechanisms and implications of, such long-term epigenetic inheritance. We argue that compelling evidence supports the long-term epigenetic inheritance of gene body methylation, at least in the model angiosperm Arabidopsis thaliana, and that variation in such methylation can therefore serve as an epigenetic basis for phenotypic variation in natural populations.}, author = {Hollwey, Elizabeth and Briffa, Amy and Howard, Martin and Zilberman, Daniel}, issn = {1879-0380}, journal = {Current Opinion in Genetics & Development}, number = {8}, publisher = {Elsevier}, title = {{Concepts, mechanisms and implications of long-term epigenetic inheritance}}, doi = {10.1016/j.gde.2023.102087}, volume = {81}, year = {2023}, } @article{14774, abstract = {Morphogen gradients impart positional information to cells in a homogenous tissue field. Fgf8a, a highly conserved growth factor, has been proposed to act as a morphogen during zebrafish gastrulation. However, technical limitations have so far prevented direct visualization of the endogenous Fgf8a gradient and confirmation of its morphogenic activity. Here, we monitor Fgf8a propagation in the developing neural plate using a CRISPR/Cas9-mediated EGFP knock-in at the endogenous fgf8a locus. By combining sensitive imaging with single-molecule fluorescence correlation spectroscopy, we demonstrate that Fgf8a, which is produced at the embryonic margin, propagates by diffusion through the extracellular space and forms a graded distribution towards the animal pole. Overlaying the Fgf8a gradient curve with expression profiles of its downstream targets determines the precise input-output relationship of Fgf8a-mediated patterning. Manipulation of the extracellular Fgf8a levels alters the signaling outcome, thus establishing Fgf8a as a bona fide morphogen during zebrafish gastrulation. Furthermore, by hindering Fgf8a diffusion, we demonstrate that extracellular diffusion of the protein from the source is crucial for it to achieve its morphogenic potential.}, author = {Harish, Rohit K and Gupta, Mansi and Zöller, Daniela and Hartmann, Hella and Gheisari, Ali and Machate, Anja and Hans, Stefan and Brand, Michael}, issn = {1477-9129}, journal = {Development}, keywords = {Developmental Biology, Molecular Biology}, number = {19}, publisher = {Company of Biologists}, title = {{Real-time monitoring of an endogenous Fgf8a gradient attests to its role as a morphogen during zebrafish gastrulation}}, doi = {10.1242/dev.201559}, volume = {150}, year = {2023}, } @article{14657, abstract = {Natural selection is usually studied between mutants that differ in reproductive rate, but are subject to the same population structure. Here we explore how natural selection acts on mutants that have the same reproductive rate, but different population structures. In our framework, population structure is given by a graph that specifies where offspring can disperse. The invading mutant disperses offspring on a different graph than the resident wild-type. We find that more densely connected dispersal graphs tend to increase the invader’s fixation probability, but the exact relationship between structure and fixation probability is subtle. We present three main results. First, we prove that if both invader and resident are on complete dispersal graphs, then removing a single edge in the invader’s dispersal graph reduces its fixation probability. Second, we show that for certain island models higher invader’s connectivity increases its fixation probability, but the magnitude of the effect depends on the exact layout of the connections. Third, we show that for lattices the effect of different connectivity is comparable to that of different fitness: for large population size, the invader’s fixation probability is either constant or exponentially small, depending on whether it is more or less connected than the resident.}, author = {Tkadlec, Josef and Kaveh, Kamran and Chatterjee, Krishnendu and Nowak, Martin A.}, issn = {1742-5662}, journal = {Journal of the Royal Society Interface}, number = {208}, publisher = {Royal Society}, title = {{Evolutionary dynamics of mutants that modify population structure}}, doi = {10.1098/rsif.2023.0355}, volume = {20}, year = {2023}, } @inproceedings{14953, abstract = {This paper provides statistical sample complexity bounds for score-matching and its applications in causal discovery. We demonstrate that accurate estimation of the score function is achievable by training a standard deep ReLU neural network using stochastic gradient descent. We establish bounds on the error rate of recovering causal relationships using the score-matching-based causal discovery method of Rolland et al. [2022], assuming a sufficiently good estimation of the score function. Finally, we analyze the upper bound of score-matching estimation within the scorebased generative modeling, which has been applied for causal discovery but is also of independent interest within the domain of generative models.y}, author = {Zhu, Zhenyu and Locatello, Francesco and Cevher, Volkan}, booktitle = {37th Conference on Neural Information Processing Systems}, isbn = {9781713899921}, issn = {1049-5258}, location = {New Orleans, LO, United States}, pages = {3325--3337}, publisher = {Neural Information Processing Systems Foundation}, title = {{Sample complexity bounds for score-matching: Causal discovery and generative modeling}}, doi = {10.52202/075280-0147}, volume = {36}, year = {2023}, }