Detailed Atomistic Simulations Research Articles

Probing the Linear-to-Plastic Transition in Polymer Nanocomposites via Atomistic Simulations: The Role of Interphases.

Polymer nanocomposites have found ubiquitous use across diverse industries, attributable to their distinctive properties and enhanced mechanical performance compared to conventional materials. Elucidating the elastic-to-plastic transition in polymer nanocomposites under diverse mechanical loads is paramount for the bespoke design of materials with desired mechanical attributes. In the current work, the elastic-to-plastic transition is probed in model systems of polyethylene oxide (PEO) and silica, SiO2, nanoparticles, through detailed atomistic molecular dynamics simulations. This comprehensive, multi-scale analysis unveils pivotal markers of the elastic-to-plastic transition, highlighting the quintessential role of microstructural and regional heterogeneities in density, strain, and stress fields, featuring the polymer-nanoparticle interphase region. At the atomic level, the behavior of polymer chains interacting with nanoparticle surfaces is traced, differentiating between free and adsorbed chains, and identifying the microscopic origins of the linear-to-plastic transition. The mechanical behavior of subregions are characterized within the PEO/SiO2 nanocomposites, focusing on the interphase and bulk-like polymer areas, probing stress heterogeneities and their decomposition into various force contributions. At the inception of plasticity, a disruption is discerned in isotropy of the polymeric density field, the emergence of low-density regions, and microscopic voids/cavities within the polymer matrix concomitant with a transition of adsorbed chains to free. The yield strain also emerges as an inflection point in the local versus global strain diagram, demarcating the elastic limit, and the plastic regime shows pronounced strain heterogeneities. The decomposition of the atomic Virial stress into bonded and non-bonded interactions indicates that the rigidity of the material is primarily governed by non-bonded interactions, significantly influenced by the volume fraction of the nanoparticle. These findings emphasize the importance of the microstructural and micromechanical environment at the polymer-nanoparticle interface on the linear-to-plastic transition, which is of great importance in the design of nanocomposite materials with advanced mechanical properties.

Data-Driven Design of Electroactive Spacer Molecules to Tune Charge Carrier Dynamics in Layered Halide Perovskite Heterostructures.

Crafting rational heterojunctions with nanostructured materials is instrumental in fostering effective interfacial charge separation and transport for optoelectronics. Layered halide perovskites (LHPs) that form heterojunctions between organic spacer molecules and inorganic metal halide layers exhibit tunable photophysics owing to their customizable band alignment. However, controlling photogenerated carrier dynamics by strategically designing layered perovskite heterojunctions remains largely unexplored. We combine a data-driven approach with time-domain density functional theory (TD-DFT) and non-adiabatic molecular dynamics (NAMD) to screen and select electronically active spacer dications (A') that introduce a type-II heterojunction in the lead iodide-based Dion-Jacobson phase LHPs. The composition-structure-electronic property correlations reveal that the number of nitrogens in aromatic heterocycles is the key factor in designing electron-accepting spacers in these perovskites. The detailed atomistic simulations validate the design strategy further by modeling (A')PbI4 perovskites, which incorporate three different screened electroactive A' spacers. The computed excited charge carrier dynamics illustrate the phonon-mediated ultrafast interfacial electron transfer from the inorganic conduction band edge to the lower-lying unoccupied orbitals of spacers, exhibiting photoluminescence quenching in these (A')PbI4 perovskites. The spatially separated electrons and holes at the type-II heterojunction interface prolong the excited charge carrier lifetime, boosting the carrier transport and exciton dynamics. Our work illustrates a robust in silico approach for designing LHPs with exciting optoelectronic properties originating from their fine-tuned heterojunctions.

Predicting mechanical heterogeneity in glassy polymer nanocomposites via an inverse computational approach based on atomistic molecular simulations and homogenization methods

Probing the mechanical behavior of the region formed between a nanoparticle reinforcement and a polymer matrix in a polymer nanocomposite structure, denoted as the “interphase”, is a main challenge as such regions are difficult to investigate by experimental methods. Here, we accurately characterize the heterogeneous mechanical behavior of polymer nanocomposites, focusing on polymer/nanofiller interphases via a combination of nanomechanical simulations and numerical homogenization techniques. Initially, the global mechanical performance of a glassy poly(ethylene oxide) polymer nanocomposite reinforced with silica nanoparticles is studied using detailed atomistic molecular dynamics simulations for 1.9% and 12.7% silica volume fractions. Next, the polymer/silica interphase thickness is identified by probing the polymer atom-based density distribution profile in the vicinity of the nanofiller at equilibrium. On the basis of this thickness, the interphase is subdivided to check the position-dependent change in mechanical properties. Then, using continuum mechanics and atomistic simulations, we proceed to compute the effective Young’s modulus and Poisson’s ratio of the polymer/nanoparticle interphase as function of the distance from the nanoparticle. In the last step, an inverse numerical homogenization model is proposed to predict the mechanical properties of the interphase on the basis of a comparison criteria with the data from MD. The results were found to be acceptable, raising the possibility of accurately and efficiently predicting interfacial properties in nanostructured materials.

Natural compound plumbagin based inhibition of hIAPP revealed by Markov state models based on MD data along with experimental validations.

Human islet amyloid polypeptide (amylin or hIAPP) is a 37 residue hormone co-secreted with insulin from β cells of the pancreas. In patients suffering from type-2 diabetes, amylin self-assembles into amyloid fibrils, ultimately leading to the death of the pancreatic cells. However, a research gap exists in preventing and treating such amyloidosis. Plumbagin, a natural compound, has previously been demonstrated to have inhibitory potential against insulin amyloidosis. Our investigation unveils collapsible regions within hIAPP that, upon collapse, facilitates hydrophobic and pi-pi interactions, ultimately leading to aggregation. Intriguingly plumbagin exhibits the ability to bind these specific collapsible regions, thereby impeding the aforementioned interactions that would otherwise drive hIAPP aggregation. We have used atomistic molecular dynamics approach to determine secondary structural changes. MSM shows metastable states forming native like hIAPP structure in presence of PGN. Our in silico results concur with in vitro results. The ThT assay revealed a striking 50% decrease in fluorescence intensity at a 1:1 ratio of hIAPP to Plumbagin. This finding suggests a significant inhibition of amyloid fibril formation by plumbagin, as ThT fluorescence directly correlates with the presence of these fibrils. Further TEM images revealed disappearance of hIAPP fibrils in plumbagin pre-treated hIAPP samples. Also, we have shown that plumbagin disrupts the intermolecular hydrogen bonding in hIAPP fibrils leading to an increase in the average beta strand spacing, thereby causing disaggregation of pre-formed fibrils demonstrating overall disruption of the aggregation machinery of hIAPP. Our work is the first to report a detailed atomistic simulation of 22 μs for hIAPP. Overall, our studies put plumbagin as a potential candidate for both preventive and therapeutic candidate for hIAPP amyloidosis.

Gas Permeation through Mechanically Resistant Self-Standing Membranes of a Neat Amorphous Organic Cage.

The synthesis and characterization of a novel film-forming organic cage and of its smaller analogue are here described. While the small cage produced single crystals suitable for X-ray diffraction studies, the large one was isolated as a dense film. Due to its remarkable film-forming properties, this latter cage could be solution processed into transparent thin-layer films and mechanically stable dense self-standing membranes of controllable thickness. Thanks to these peculiar features, the membranes were also successfully tested for gas permeation, reporting a behavior similar to that found with stiff glassy polymers such as polymers of intrinsic microporosity or polyimides. Given the growing interest in the development of molecular-based membranes, for example for separation technologies and functional coatings, the properties of this organic cage were investigated by thorough analysis of their structural, thermal, mechanical and gas transport properties, and by detailed atomistic simulations.

Hierarchical Hydration Dynamics of RNA with Nano-Water-Pool at Its Core.

Many functional RNAs fold into a compact, roughly globular shape by minimizing the electrostatic repulsion between their negatively charged phosphodiester backbone. The fold of such close, compact RNA architecture is often so designed that its outer surface and complex core both are predominately populated by phosphate groups loosely sequestering bases in the intermediate layers. A number of helical junctions maintain the RNA core and its nano-water-pool. While the folding of RNA is manifested by its counterion environment composed of mixed mono- and divalent salts, the concerted role of ion and water in maintaining an RNA fold is yet to be explored. In this work, detailed atomistic simulations of SAM-I and Add Adenine riboswitch aptamers, and subgenomic flavivirus RNA (sfRNA) have been performed in a physiological mixed mono- and divalent salt environment. All three RNA systems have compact folds with a core diameter of range 1-1.7 nm. The spatiotemporal heterogeneity of RNA hydration was probed in a layer-wise manner by distinguishing the core, the intermediate, and the outer layers. The layer-wise decomposition of hydrogen bonds and collective single-particle reorientational dynamics reveal a nonmonotonic relaxation pattern with the slowest relaxation observed at the intermediate layers that involves functionally important tertiary motifs. The slowness of this intermediate layer is attributed to two types of long-resident water molecules: (i) water from ion-hydration layers and (ii) structurally trapped water (distant from ions). The relaxation kinetics of the core and the surface water essentially exposed to the phosphate groups show well-separated time scales from the intermediate layers. In the slow intermediate layers, site-specific ions and water control the functional dynamics of important RNA motifs like kink-turn, observed in different structure-probing experiments. Most interestingly, we find that as the size of the RNA core increases (SAM1 core < sfRNAcore < Add adenine core), its hydration tends to show faster relaxation. The hierarchical hydration and the layer-wise base-phosphate composition uniquely portray the globular RNA to act like a soft vesicle with a quasi-dynamic nano-water-pool at its core.

In Silico Demonstration of Fast Anhydrous Proton Conduction on Graphanol.

Development of new materials capable of conducting protons in the absence of water is crucial for improving the performance, reducing the cost, and extending the operating conditions for proton exchange membrane fuel cells. We present detailed atomistic simulations showing that graphanol (hydroxylated graphane) will conduct protons anhydrously with very low diffusion barriers. We developed a deep learning potential (DP) for graphanol that has near-density functional theory accuracy but requires a very small fraction of the computational cost. We used our DP to calculate proton self-diffusion coefficients as a function of temperature, to estimate the overall barrier to proton diffusion, and to characterize the impact of thermal fluctuations as a function of system size. We propose and test a detailed mechanism for proton conduction on the surface of graphanol. We show that protons can rapidly hop along Grotthuss chains containing several hydroxyl groups aligned such that hydrogen bonds allow for conduction of protons forward and backward along the chain without hydroxyl group rotation. Long-range proton transport only takes place as new Grotthuss chains are formed by rotation of one or more hydroxyl groups in the chain. Thus, the overall diffusion barrier consists of a convolution of the intrinsic proton hopping barrier and the intrinsic hydroxyl rotation barrier. Our results provide a set of design rules for developing new anhydrous proton conducting membranes with even lower diffusion barriers.

Dynamics of Polymer Chains in Poly(ethylene oxide)/Silica Nanocomposites via a Combined Computational and Experimental Approach.

The dynamics of polymer chains in poly(ethylene oxide)/silica (PEO/SiO2) nanoparticle nanohybrids have been investigated via a combined computational and experimental approach involving atomistic molecular dynamics simulations and dielectric relaxation spectroscopy (DRS) measurements. The complementarity of the approaches allows us to study systems with different polymer molecular weights, nanoparticle radii, and compositions across a broad range of temperatures. We study the effects of spatial confinement, which is induced by the nanoparticles, and chain adsorption on the polymer's structure and dynamics. The investigation of the static properties of the nanocomposites via detailed atomistic simulations revealed a heterogeneous polymer density layer at the vicinity of the PEO/SiO2 interface that exhibited an intense maximum close to the inorganic surface, whereas the bulk density was reached for distances ∼1-1.2 nm away from the nanoparticle. For small volume fractions of nanoparticles, the polymer dynamics, probed by the atomistic simulations of low-molecular-weight chains at high temperatures, are consistent with the presence of a thin adsorbed layer that exhibits slow dynamics, with the dynamics far away from the nanoparticle being similar to those in the bulk. However, for high volume fractions of nanoparticles (strong confinement), the dynamics of all polymer chains were predicted slower than that in the bulk. On the other hand, similar dynamics were found experimentally for both the local β-process and the segmental dynamics for high-molecular-weight systems measured at temperatures below the melting temperature of the polymer, which were probed by DRS. These differences can be attributed to various parameters, including systems of different molecular weights and nanoparticle states of dispersion, the different temperature range studied by the different methods, the potential presence of a reduced-mobility PEO/SiO2 interfacial layer that does not contribute to the dielectric spectrum, and the presence of amorphous-crystalline interfaces in the experimental samples that may lead to a different dynamical behaviors of the PEO chains.

Mechanical Behavior of Polymer Nanocomposites via Atomistic Simulations: Conformational Heterogeneity and the Role of Strain Rate.

Polymer nanohybrids with a high fraction of nanofillers have been found to exhibit improved mechanical properties compared to the neat polymer homogeneous systems. Since polymer-based materials are characterized by a broad range of relaxation times, it is expected that their response under external load would depend on the actual rate of the applied deformation. In this work, we investigate the heterogeneous mechanical behavior in glassy poly(ethylene oxide)/silica nanoparticles (PEO/SiO2) nanocomposites via detailed atomistic molecular dynamics simulations. Our goal is to unravel the effect of strain rate on the mechanism of polymer nanocomposite reinforcement, within the glassy state, by directly probing the mechanical properties at the molecular level. By applying tensile deformations with various strain rates we clearly demonstrate that the value of the applied strain rate strongly affects the mechanical properties of the PEO/SiO2 model systems, inducing a transition from a rubber-like behavior, at low strain rate, to a more brittle one, at high strain rates. Then, we further investigate the mechanical heterogeneity in glassy PEO/SiO2 systems by probing directly the stress and strain fields for various conformations of adsorbed (trains, tails, loops, and bridges), and free polymer chains. Our data emphasize the importance of both train and bridge conformations on the mechanical reinforcement of the polymer nanocomposites. Last, we also probe the mobility of various chain conformations, under different applied strain rates, and their contribution to the overall mechanical behavior of the nanocomposite, during the deformation process.

Accelerated Diffusion Following Deprotection in Chemically Amplified Resists.

Polymeric chemically amplified resists (CARs) are critical materials for high-throughput lithographic processes. A photoactivated acid-anion catalyst changes the polymer's solubility via a deprotection reaction, which enables pattern development through selective dissolution. To capture observed reaction kinetics, reaction-diffusion models employ a catalyst diffusivity that is accelerated by reaction. However, the microscopic origin and factors contributing to this phenomena remain unclear. Herein, we employ detailed atomistic molecular dynamics simulations to examine the impact of protecting group removal and material relaxation on catalyst mobility. We report data on polymer density, catalyst dispersion, excess free volume, and segmental dynamics with increasing time/extent of deprotection. We then propose simple kinetic Monte Carlo algorithms that can describe both molecular dynamics simulations of deprotection reactions and experimental data.

On the influence of Al-concentration on the fracture toughness of NiAl: Microcantilever fracture tests and atomistic simulations

While the mechanical properties of the stoichiometric B2 β-phase of NiAl are well established, the effect of off-stoichiometric composition on the fracture behavior has not yet been systematically studied over the entire composition range of 40-55% Al. Here we use microbending tests on notched cantilever beams FIB-milled from NiAl single crystals with an aluminized as well as an oxidation-induced composition gradient to determine the influence of the Al concentration on the fracture toughness. Small variations from the stoichiometric composition are shown to lead to drastic changes in fracture behavior. Ni-rich NiAl shows a gradual decrease of fracture toughness and plastic energy dissipation with Ni-concentration that follows the inverse proportionality to the yield stress postulated by the Ritchie, Knott and Rice (RKR) model. Al-rich NiAl, however, shows a concentration-independent low fracture toughness and no signs of accompanying plasticity - both which can not be described by the RKR model. Detailed atomistic simulations show that rather than affecting the theoretical fracture toughness according to the thermodynamic Griffith criterion, structural vacancies in Al-rich NiAl lower the kinetic barrier for crack propagation, explaining the observed different concentration dependence of the fracture toughness. Similar observations are also expected for fracture in other off-stoichiometric B2 aluminides.

Mechanical properties of glassy polymer nanocomposites via atomistic and continuum models: The role of interphases

The mechanical properties of glassy polymer nanocomposites (PNCs) are investigated via a new hierarchical computational methodology, which combines atomistic molecular dynamics (MD) simulations and homogenization approaches. The homogenization methodology is based on a systematic nano/micro/macro coupling between detailed atomistic MD simulations and a variational approach through the Hill–Mandel lemma. The proposed methodology is applied in model glassy polybutadiene/silica nanocomposites for different nanoparticle (NP) volume fractions. Initially, using MD simulations, the polymer/NP interphase in PNCs is directly examined by probing the density distribution and the stress profile at equilibrium. By using a continuum mechanics based approach, we can compute effective deformation gradients for each atom, allowing us to probe the distribution of the (local) stress and strain fields in the atomistic model. With this new approach, the effective Young modulus and the Poisson ratio of the polymer/NP interphases are directly calculated, exhibiting a higher rigidity compared to the polymer matrix. In the last part of the proposed approach, the mechanical properties at the interphases and the polymer are used together with the homogenization approach to develop a continuum model for predicting the mechanical properties of the PNCs, which are found to be in very good agreement with the effective mechanical properties calculated through atomistic MD simulations.

Atomic insights on intermixing of nanoscale nitride multilayer triggered by nanoindentation

Mechanical properties of nanoscale multilayer coatings are to a large extent governed by the number of interfaces and their characteristics. While for a reduced layer thickness, increasing strength and toughness values have been reported, properties degrade for layer thicknesses of just several nanometers. Here, we report on an entirely overlooked phenomenon occurring during the indentation of nanolayers, presumably explaining the degradation of properties. Nanoindentation, commonly used to determine properties of hard coatings, is found to disrupt and intermix the multilayer structure due to the deformation imposed. Detailed electron microscopy studies and atomistic simulations provide evidence for intermixing in an epitaxial transition metal nitride superlattice thin film induced by nanoindentation. The formation of a solid solution reduces the interfacial density and leads to a sharp drop in the dislocation density. Our results confirm that plastic deformation causes the microstructure instability of nitride multilayer, which may further improve our understanding of multilayer strength mechanisms.

Salt ions induced transport properties of poly(methacrylic acid) PMA in aqueous solutions studied by molecular dynamics simulations

Abstract The detailed atomistic molecular dynamics simulations studies of multiple chain anionic polyelectrolyte poly(methacrylic acid) PMA was carried out to investigate the effect of monovalent salt i.e. NaCl on the transport properties such as self-diffusion coefficient of PMA, salt-ions and water molecules in dilute aqueous solutions. The concentration of NaCl was varied from dilute (0.01 M) to concentrated regime (1 M) and the self-diffusion coefficient of PMA, salt-ions and water molecules were calculated. The self-diffusion coefficient of PMA decreases with increase in NaCl concentration in agreement with experiments and previous reported simulations works. The salt ions i.e. Na+ and Cl- too shows a decrease in their self-diffusion coefficient values with NaCl concentration. The self-diffusivity of water molecules was found to be decrease in a linear manner with salt concentration. The salt-ions exhibits faster movement as compared to PMA and water on account of a higher self-diffusion coefficient.

Molecular dynamics simulations studies of transport properties of sodium-(polymethacrylate) (Na+-PMA) in aqueous solutions – Effect of salt concentration

Abstract The detailed atomistic molecular dynamics simulations studies of multiple chain anionic polyelectrolyte sodium-(polymethacrylate) (Na+-PMA) was carried out to investigate the effect of monovalent salt i.e. NaCl on the transport properties such as self-diffusivity of sodium salt i.e., Na+-PMA, salt-ions and water molecules in dilute aqueous solutions. The NaCl concentration was varied from dilute (0.01 M) regime to concentrated regime (1 M) and the self-diffusivity of Na+-PMA, salt-ions and water molecules were obtained. The Na+-PMA diffusivity shows a decrease with added salt concentration in excellent agreement with related experiments and previous reported atomistic and model simulation works. The salt ions i.e. Na+ and Cl− too shows a decrease in their self-diffusivity values with NaCl concentration. The self-diffusivity of water molecules was found to be decrease in a linear manner with salt concentration. The salt-ions exhibits faster movement as compared to Na+-PMA and water on account of a higher self-diffusion coefficient.

A multiscale coarse-grained model of the SARS-CoV-2 virion

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and ongoing development of a largely “bottom-up” coarse-grained (CG) model of the SARS-CoV-2 virion. Data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data become publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion.

Phase Boundary and Salt Partitioning in Coacervate Complexes Formed between Poly(acrylic acid) and Poly(N,N-dimethylaminoethyl methacrylate) from Detailed Atomistic Simulations Combined with Free Energy Perturbation and Thermodynamic Integration Calculations

With the help of detailed atomistic simulations, we have determined the phase boundary of a complex coacervate system resulting from the complexation of two oppositely and fully charged weak polyel...

Photoinduced Dynamics of Charge Carriers in Metal Halide Perovskites from an Atomistic Perspective.

Perovskite solar cells have attracted intense attention over the past decade because of their low cost, abundant raw materials, and rapidly growing power conversion efficiency (PCE). However, nonradiative charge carrier losses still constitute a major factor limiting the PCE to well below the Shockley-Queisser limit. This Perspective summarizes recent atomistic quantum dynamics studies on the photoinduced excited-state processes in metal halide perovskites (MHPs), including both hybrid organic-inorganic and all-inorganic MHPs and three- and two-dimensional MHPs. The simulations, performed using a combination of time-domain ab initio density functional theory and nonadiabatic molecular dynamics, allow emphasis on various intrinsic and extrinsic features, such as components, structure, dimensionality and interface engineering, control and exposure to various environmental factors, defects, surfaces, and their passivation. The detailed atomistic simulations advance our understanding of electron-vibrational dynamics in MHPs and provide valuable guidelines for enhancing the performance of perovskite solar cells.

Aggregation of methacrylate-based ternary biomimetic antimicrobial polymers in solution

Using detailed atomistic simulations, we explore the morphological characteristics of aggregates formed in solution phase by ternary biomimetic antimicrobial (AM) methacrylate polymers, composed of hydrophobic, charged cationic and polar functional groups and compare it with aggregate morphologies of binary methacrylate polymers, composed only of hydrophobic and charged cationic functional groups. The effect of sequence of the constituent functional groups on aggregate conformation is also studied by considering random and block sequences along the polymer backbone. Our results show that while binary polymers tend to form robust aggregates, replacing some of the hydrophobic groups with overall charge neutral polar groups weakens the aggregate considerably, leading to increased conformational fluctuations and formation of loose-packed, open aggregates, particularly in the case of random ternary polymers. Interaction energy calculations clearly suggest that the role of inclusion of polar groups in ternary polymers is two-fold: (1) to reduce possible strong local concentration of hydrophobic groups and ‘smear’ the overall hydrophobicity along the polymer backbone to increase the solubility of the polymers (2) to compensate the loss of attractive hydrophobic interactions by forming attractive electrostatic interactions with the charged groups and contribute to aggregation formation, albeit weak. Given that most of the naturally occurring AM peptides have contributions from all the three functional groups, this study elucidates the functionally tuneable role of inclusion of polar groups in the way AM agents interact with each other in solution phase, which can eventually dictate their partitioning behaviour into bacterial and mammalian membranes.

Chalcogen-Assisted Enhanced Atomic Orbital Interaction at TMD–Metal Interface and Sulfur Passivation for Overall Performance Boost of 2-D TMD FETs

Metal-semiconductor interface is a bottleneck for efficient transport of charge carriers through Transition Metal Dichalcogenide (TMD) based field-effect transistors (FETs). Injection of charge carriers across such interfaces is mostly limited by Schottky barrier at the contacts which must be reduced to achieve highly efficient contacts for carrier injection into the channel. Here we introduce a universal approach involving dry chemistry to enhance atomic orbital interaction between various TMDs (MoS2, WS2, MoSe2 and WSe2) & metal contacts has been experimentally demonstrated. Quantum chemistry between TMDs, Chalcogens and metals has been explored using detailed atomistic (DFT & NEGF) simulations, which is then verified using Raman, PL and XPS investigations. Atomistic investigations revealed lower contact resistance due to enhanced orbital interaction and unique physics of charge sharing between constituent atoms in TMDs with introduced Chalcogen atoms which is subsequently validated through experiments. Besides contact engineering, which lowered contact resistance by 72, 86, 1.8, 13 times in MoS2, WS2, MoSe2 and WSe2 respectively, a novel approach to cure / passivate dangling bonds present at the 2D TMD channel surface has been demonstrated. While the contact engineering improved the ON-state performance (ION, gm, mobility and RON) of 2D TMD FETs by orders of magnitude, Chalcogen based channel passivation was found to improve gate control (IOFF, SS, & VTH) significantly. This resulted in an overall performance boost. The engineered TMD FETs were shown to have performance on par with best reported till date.