@article{9677, abstract = {Progress in the atomic-scale modeling of matter over the past decade has been tremendous. This progress has been brought about by improvements in methods for evaluating interatomic forces that work by either solving the electronic structure problem explicitly, or by computing accurate approximations of the solution and by the development of techniques that use the Born–Oppenheimer (BO) forces to move the atoms on the BO potential energy surface. As a consequence of these developments it is now possible to identify stable or metastable states, to sample configurations consistent with the appropriate thermodynamic ensemble, and to estimate the kinetics of reactions and phase transitions. All too often, however, progress is slowed down by the bottleneck associated with implementing new optimization algorithms and/or sampling techniques into the many existing electronic-structure and empirical-potential codes. To address this problem, we are thus releasing a new version of the i-PI software. This piece of software is an easily extensible framework for implementing advanced atomistic simulation techniques using interatomic potentials and forces calculated by an external driver code. While the original version of the code (Ceriotti et al., 2014) was developed with a focus on path integral molecular dynamics techniques, this second release of i-PI not only includes several new advanced path integral methods, but also offers other classes of algorithms. In other words, i-PI is moving towards becoming a universal force engine that is both modular and tightly coupled to the driver codes that evaluate the potential energy surface and its derivatives.}, author = {Kapil, Venkat and Rossi, Mariana and Marsalek, Ondrej and Petraglia, Riccardo and Litman, Yair and Spura, Thomas and Cheng, Bingqing and Cuzzocrea, Alice and Meißner, Robert H. and Wilkins, David M. and Helfrecht, Benjamin A. and Juda, Przemysław and Bienvenue, Sébastien P. and Fang, Wei and Kessler, Jan and Poltavsky, Igor and Vandenbrande, Steven and Wieme, Jelle and Corminboeuf, Clemence and Kühne, Thomas D. and Manolopoulos, David E. and Markland, Thomas E. and Richardson, Jeremy O. and Tkatchenko, Alexandre and Tribello, Gareth A. and Van Speybroeck, Veronique and Ceriotti, Michele}, issn = {0010-4655}, journal = {Computer Physics Communications}, pages = {214--223}, publisher = {Elsevier}, title = {{i-PI 2.0: A universal force engine for advanced molecular simulations}}, doi = {10.1016/j.cpc.2018.09.020}, volume = {236}, year = {2019}, } @article{9680, abstract = {Atomistic modeling of phase transitions, chemical reactions, or other rare events that involve overcoming high free energy barriers usually entails prohibitively long simulation times. Introducing a bias potential as a function of an appropriately chosen set of collective variables can significantly accelerate the exploration of phase space, albeit at the price of distorting the distribution of microstates. Efficient reweighting to recover the unbiased distribution can be nontrivial when employing adaptive sampling techniques such as metadynamics, variationally enhanced sampling, or parallel bias metadynamics, in which the system evolves in a quasi-equilibrium manner under a time-dependent bias. We introduce an iterative unbiasing scheme that makes efficient use of all the trajectory data and that does not require the distribution to be evaluated on a grid. The method can thus be used even when the bias has a high dimensionality. We benchmark this approach against some of the existing schemes on model systems with different complexity and dimensionality.}, author = {Giberti, F. and Cheng, Bingqing and Tribello, G. A. and Ceriotti, M.}, issn = {1549-9626}, journal = {Journal of Chemical Theory and Computation}, number = {1}, pages = {100--107}, publisher = {American Chemical Society}, title = {{Iterative unbiasing of quasi-equilibrium sampling}}, doi = {10.1021/acs.jctc.9b00907}, volume = {16}, year = {2019}, } @article{12600, abstract = {The snow cover dynamics of High Mountain Asia are usually assessed at spatial resolutions of 250 m or greater, but this scale is too coarse to clearly represent the rugged topography common to the region. Higher-resolution measurement of snow-covered area often results in biased sampling due to cloud cover and deep shadows. We therefore develop a Normalized Difference Snow Index-based workflow to delineate snow lines from Landsat Thematic Mapper/Enhanced Thematic Mapper+ imagery and apply it to the upper Langtang Valley in Nepal, processing 194 scenes spanning 1999 to 2013. For each scene, we determine the spatial distribution of snow line altitudes (SLAs) with respect to aspect and across six subcatchments. Our results show that the mean SLA exhibits distinct seasonal behavior based on aspect and subcatchment position. We find that SLA dynamics respond to spatial and seasonal trade-offs in precipitation, temperature, and solar radiation, which act as primary controls. We identify two SLA spatial gradients, which we attribute to the effect of spatially variable precipitation. Our results also reveal that aspect-related SLA differences vary seasonally and are influenced by solar radiation. In terms of seasonal dominant controls, we demonstrate that the snow line is controlled by snow precipitation in winter, melt in premonsoon, a combination of both in postmonsoon, and temperature in monsoon, explaining to a large extent the spatial and seasonal variability of the SLA in the upper Langtang Valley. We conclude that while SLA and snow-covered area are complementary metrics, the SLA has a strong potential for understanding local-scale snow cover dynamics and their controlling mechanisms.}, author = {Girona‐Mata, Marc and Miles, Evan S. and Ragettli, Silvan and Pellicciotti, Francesca}, issn = {1944-7973}, journal = {Water Resources Research}, keywords = {Water Science and Technology}, number = {8}, pages = {6754--6772}, publisher = {American Geophysical Union}, title = {{High‐resolution snowline delineation from Landsat imagery to infer snow cover controls in a Himalayan catchment}}, doi = {10.1029/2019wr024935}, volume = {55}, year = {2019}, } @article{12602, abstract = {This study aims at developing and applying a spatially-distributed coupled glacier mass balance and ice-flow model to attribute the response of glaciers to natural and anthropogenic climate change. We focus on two glaciers with contrasting surface characteristics: a debris-covered glacier (Langtang Glacier in Nepal) and a clean-ice glacier (Hintereisferner in Austria). The model is applied from the end of the Little Ice Age (1850) to the present-day (2016) and is forced with four bias-corrected General Circulation Models (GCMs) from the historical experiment of the CMIP5 archive. The selected GCMs represent region-specific warm-dry, warm-wet, cold-dry, and cold-wet climate conditions. To isolate the effects of anthropogenic climate change on glacier mass balance and flow runs from these GCMs with and without further anthropogenic forcing after 1970 until 2016 are selected. The outcomes indicate that both glaciers experience the largest reduction in area and volume under warm climate conditions, whereas area and volume reductions are smaller under cold climate conditions. Simultaneously with changes in glacier area and volume, surface velocities generally decrease over time. Without further anthropogenic forcing the results reveal a 3% (9%) smaller decline in glacier area (volume) for the debris-covered glacier and a 18% (39%) smaller decline in glacier area (volume) for the clean-ice glacier. The difference in the magnitude between the two glaciers can mainly be attributed to differences in the response time of the glaciers, where the clean-ice glacier shows a much faster response to climate change. We conclude that the response of the two glaciers can mainly be attributed to anthropogenic climate change and that the impact is larger on the clean-ice glacier. The outcomes show that the model performs well under different climate conditions and that the developed approach can be used for regional-scale glacio-hydrological modeling.}, author = {Wijngaard, René R. and Steiner, Jakob F. and Kraaijenbrink, Philip D. A. and Klug, Christoph and Adhikari, Surendra and Banerjee, Argha and Pellicciotti, Francesca and van Beek, Ludovicus P. H. and Bierkens, Marc F. P. and Lutz, Arthur F. and Immerzeel, Walter W.}, issn = {2296-6463}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, title = {{Modeling the response of the Langtang Glacier and the Hintereisferner to a changing climate since the Little Ice Age}}, doi = {10.3389/feart.2019.00143}, volume = {7}, year = {2019}, } @article{12601, abstract = {Ice cliffs and ponds on debris-covered glaciers have received increased attention due to their role in amplifying local melt. However, very few studies have looked at these features on the catchment scale to determine their patterns and changes in space and time. We have compiled a detailed inventory of cliffs and ponds in the Langtang catchment, central Himalaya, from six high-resolution satellite orthoimages and DEMs between 2006 and 2015, and a historic orthophoto from 1974. Cliffs cover between 1.4% (± 0.4%) in the dry and 3.4% (± 0.9%) in the wet seasons and ponds between 0.6% (± 0.1%) and 1.6% (± 0.3%) of the total debris-covered tongues. We find large variations between seasons, as cliffs and ponds tend to grow in the wetter monsoon period, but there is no obvious trend in total area over the study period. The inventory further shows that cliffs are predominately north-facing irrespective of the glacier flow direction. Both cliffs and ponds appear in higher densities several hundred metres from the terminus in areas where tributaries reach the main glacier tongue. On the largest glacier in the catchment ~10% of all cliffs and ponds persisted over nearly a decade.}, author = {STEINER, JAKOB F. and BURI, PASCAL and MILES, EVAN S. and RAGETTLI, SILVAN and Pellicciotti, Francesca}, issn = {1727-5652}, journal = {Journal of Glaciology}, number = {252}, pages = {617--632}, publisher = {Cambridge University Press}, title = {{Supraglacial ice cliffs and ponds on debris-covered glaciers: Spatio-temporal distribution and characteristics}}, doi = {10.1017/jog.2019.40}, volume = {65}, year = {2019}, } @article{12192, abstract = {Transposable elements (TEs), the movement of which can damage the genome, are epigenetically silenced in eukaryotes. Intriguingly, TEs are activated in the sperm companion cell – vegetative cell (VC) – of the flowering plant Arabidopsis thaliana. However, the extent and mechanism of this activation are unknown. Here we show that about 100 heterochromatic TEs are activated in VCs, mostly by DEMETER-catalyzed DNA demethylation. We further demonstrate that DEMETER access to some of these TEs is permitted by the natural depletion of linker histone H1 in VCs. Ectopically expressed H1 suppresses TEs in VCs by reducing DNA demethylation and via a methylation-independent mechanism. We demonstrate that H1 is required for heterochromatin condensation in plant cells and show that H1 overexpression creates heterochromatic foci in the VC progenitor cell. Taken together, our results demonstrate that the natural depletion of H1 during male gametogenesis facilitates DEMETER-directed DNA demethylation, heterochromatin relaxation, and TE activation.}, author = {He, Shengbo and Vickers, Martin and Zhang, Jingyi and Feng, Xiaoqi}, issn = {2050-084X}, journal = {eLife}, keywords = {General Immunology and Microbiology, General Biochemistry, Genetics and Molecular Biology, General Medicine, General Neuroscience}, publisher = {eLife Sciences Publications, Ltd}, title = {{Natural depletion of histone H1 in sex cells causes DNA demethylation, heterochromatin decondensation and transposon activation}}, doi = {10.7554/elife.42530}, volume = {8}, year = {2019}, } @article{12190, abstract = {Meiotic crossover frequency varies within genomes, which influences genetic diversity and adaptation. In turn, genetic variation within populations can act to modify crossover frequency in cis and trans. To identify genetic variation that controls meiotic crossover frequency, we screened Arabidopsis accessions using fluorescent recombination reporters. We mapped a genetic modifier of crossover frequency in Col × Bur populations of Arabidopsis to a premature stop codon within TBP-ASSOCIATED FACTOR 4b (TAF4b), which encodes a subunit of the RNA polymerase II general transcription factor TFIID. The Arabidopsis taf4b mutation is a rare variant found in the British Isles, originating in South-West Ireland. Using genetics, genomics, and immunocytology, we demonstrate a genome-wide decrease in taf4b crossovers, with strongest reduction in the sub-telomeric regions. Using RNA sequencing (RNA-seq) from purified meiocytes, we show that TAF4b expression is meiocyte enriched, whereas its paralog TAF4 is broadly expressed. Consistent with the role of TFIID in promoting gene expression, RNA-seq of wild-type and taf4b meiocytes identified widespread transcriptional changes, including in genes that regulate the meiotic cell cycle and recombination. Therefore, TAF4b duplication is associated with acquisition of meiocyte-specific expression and promotion of germline transcription, which act directly or indirectly to elevate crossovers. This identifies a novel mode of meiotic recombination control via a general transcription factor.}, author = {Lawrence, Emma J. and Gao, Hongbo and Tock, Andrew J. and Lambing, Christophe and Blackwell, Alexander R. and Feng, Xiaoqi and Henderson, Ian R.}, issn = {0960-9822}, journal = {Current Biology}, keywords = {General Agricultural and Biological Sciences, General Biochemistry, Genetics and Molecular Biology}, number = {16}, pages = {2676--2686.e3}, publisher = {Elsevier BV}, title = {{Natural variation in TBP-ASSOCIATED FACTOR 4b controls meiotic crossover and germline transcription in Arabidopsis}}, doi = {10.1016/j.cub.2019.06.084}, volume = {29}, year = {2019}, } @unpublished{8305, abstract = {In this paper, we present the first fully asynchronous distributed key generation (ADKG) algorithm as well as the first distributed key generation algorithm that can create keys with a dual (f,2f+1)−threshold that are necessary for scalable consensus (which so far needs a trusted dealer assumption). In order to create a DKG with a dual (f,2f+1)− threshold we first answer in the affirmative the open question posed by Cachin et al. how to create an AVSS protocol with recovery thresholds f+1