Search

The experimentally observed effect of selective deuterium substitution on the open circuit voltage for a blend of poly(3-hexylthiophene) (P3HT) and [6,6]-phenyl-C61-butyric acid methyl ester (PCBM; Nat. Commun. 2014, 5, 3180) is explored using a 221-atom model of a polymer-wrapped PCBM molecule. The protonic and deuteronic wave functions for the H/D isotopologues of the hexyl side chains are described within a quantum trajectory/electronic structure approach where the dynamics is performed with newly developed nonlinear corrections to the quantum forces, necessary to describe the nuclear wave functions; the classical forces are generated with a density functional tight binding method. The resulting protonic and deuteronic time-dependent wave functions are used to assess the effects of isotopic substitution (deuteration) on the energy gaps relevant to the charge transfer for the donor and acceptor electronic states. While the isotope effect on the electronic energy levels is found negligible, the quantum-induced fluctuations of the energy gap between the charge transfer and charge separated states due to nuclear wave functions may account for experimental trends by promoting charge transfer in P3HT:PCBM and increasing charge recombination on the donor in the deuterium substituted P3HT:PCBM.

Self-assembly of block copolymers provides numerous opportunities to create functional materials, utilizing self-assembled microdomains with a variety of morphology and periodic architectures as templates for functional nanofillers. Here new progress is reported toward the fabrication of thermally responsive and electrically conductive polymeric self-assemblies made from a water-soluble poly(thiophene) derivative with short poly(ethylene oxide) side chains and Pluronic L62 block copolymer solution in water. The structural and electrical properties of conjugated polymer-embedded self-assembled architectures are investigated by combining small-angle neutron and X-ray scattering, coarse-grained molecular dynamics simulations, and impedance spectroscopy. The L62 solution template organizes the conjugated polymers by stably incorporating them into the hydrophilic domains thus inhibiting aggregation. The changing morphology of L62 during the micellar-to-lamellar phase transition defines the embedded conjugated polymer network. As a result, the conductivity is strongly coupled to the structural change of the templating L62 phase and exhibits thermally reversible behavior with no signs of quenching of the conductivity at high temperature. This study shows promise for enabling more flexibility in processing and utilizing water-soluble conjugated polymers in aqueous solutions for self-assembly based fabrication of stimuli-responsive nanostructures and sensory materials., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Bulk molecular ionic solids exhibit fascinating electronic properties, including electron correlations, phase transitions, and superconducting ground states. In contrast, few of these phenomena have been observed in low-dimensional molecular structures, including thin films, nanoparticles, and molecular blends, not in the least because most of such structures have been composed of nearly closed-shell molecules. It is therefore desirable to develop low-dimensional ionic molecular structures that can capture potential applications. Here, we present detailed analysis of monolayer-thick structures of the canonical TTF-TCNQ (tetrathiafulvalene 7,7,8,8-tetracyanoquinodimethane) system grown on low-index gold and silver surfaces. The most distinctive property of the epitaxial growth is the wide abundance of stable TTF/TCNQ ratios, in sharp contrast to the predominance of a 1:1 ratio in the bulk. We propose the existence of the surface phase diagram that controls the structures of TTF-TCNQ on the surfaces and demonstrate phase transitions that occur upon progressively increasing the density of TCNQ while keeping the surface coverage of TTF fixed. Based on direct observations, we propose the binding motif behind the stable phases and infer the dominant interactions that enable the existence of the rich spectrum of surface structures. Finally, we also show that the surface phase diagram will control the epitaxy beyond monolayer coverage. Multiplicity of stable surface structures, the corollary rich phase diagram, and the corresponding phase transitions present an interesting opportunity for low-dimensional molecular systems, particularly if some of the electronic properties of the bulk can be preserved or modified in the surface phases.

Defect engineering has been a critical step in controlling the transport characteristics of electronic devices, and the ability to create, tune, and annihilate defects is essential to enable the range of next-generation devices. Whereas defect formation has been well-demonstrated in three-dimensional semiconductors, similar exploration of the heterogeneity in atomically thin two-dimensional semiconductors and the link between their atomic structures, defects, and properties has not yet been extensively studied. Here, we demonstrate the growth of MoSe2-x single crystals with selenium (Se) vacancies far beyond intrinsic levels, up to ∼20%, that exhibit a remarkable transition in electrical transport properties from n- to p-type character with increasing Se vacancy concentration. A new defect-activated phonon band at ∼250 cm(-1) appears, and the A1g Raman characteristic mode at 240 cm(-1) softens toward ∼230 cm(-1) which serves as a fingerprint of vacancy concentration in the crystals. We show that post-selenization using pulsed laser evaporated Se atoms can repair Se-vacant sites to nearly recover the properties of the pristine crystals. First-principles calculations reveal the underlying mechanisms for the corresponding vacancy-induced electrical and optical transitions.

The realization of controllable morphologies of bulk heterojunctions (BHJ) in organic photovoltaics (OPVs) is one of the key factors enabling high-efficiency devices. We provide new insights into the fundamental mechanisms essential for the optimization of power conversion efficiencies (PCEs) with additive processing to PBDTTT-CF:PC71BM system. We have studied the underlying mechanisms by monitoring the 3D nanostructural modifications in BHJs and correlated the modifications with the optical analysis and theoretical modeling of charge transport. Our results demonstrate profound effects of diiodooctane (DIO) on morphology and charge transport in the active layers. For small amounts of DIO (<3 vol %), DIO promotes the formation of a well-mixed donor-acceptor compact film and augments charge transfer and PCE. In contrast, for large amounts of DIO (>3 vol %), DIO facilitates a loosely packed mixed morphology with large clusters of PC71BM, leading to deterioration in PCE. Theoretical modeling of charge transport reveals that DIO increases the mobility of electrons and holes (the charge carriers) by affecting the energetic disorder and electric field dependence of the mobility. Our findings show the implications of phase separation and carrier transport pathways to achieve optimal device performances.

Development of devices and structures based on the layered 2D materials critically hinges on the capability to induce, control, and tailor the electronic, transport, and optoelectronic properties via defect engineering, much like doping strategies have enabled semiconductor electronics and forging enabled introduction the of iron age. Here, we demonstrate the use of a scanning helium ion microscope (HIM) for tailoring the functionality of single layer MoSe2 locally, and decipher associated mechanisms at the atomic level. We demonstrate He(+) beam bombardment that locally creates vacancies, shifts the Fermi energy landscape and increases the Young's modulus of elasticity. Furthermore, we observe for the first time, an increase in the B-exciton photoluminescence signal from the nanoforged regions at the room temperature. The approach for precise defect engineering demonstrated here opens opportunities for creating functional 2D optoelectronic devices with a wide range of customizable properties that include operating in the visible region.

The recent synthesis of organic molecular liquids with permanent porosity opens up exciting new avenues for gas capture, storage, and separation. Using molecular simulations, we study the thermodynamics and kinetics for the storage of CH4, CO2, and N2 molecules in porous liquids consisting of crown-ether-substituted cage molecules in a 15-crown-5 solvent. It is found that the intrinsic gas storage capacity per cage molecule follows the order CH4 > CO2 > N2, which does not correlate simply with the size of gas molecules. Different gas molecules are stored inside the cage differently; e.g., CO2 molecules prefer the cage's core whereas CH4 molecules favor both the core and the branch regions. All gas molecules considered can enter the cage essentially without energy barriers and leave the cage on a nanosecond time scale by overcoming a modest energy penalty. The molecular mechanisms of these observations are clarified.

Understanding how additives interact and segregate within bulk heterojunction (BHJ) thin films is critical for exercising control over structure at multiple length scales and delivering improvements in photovoltaic performance. The morphological evolution of poly(3-hexylthiophene) (P3HT) and phenyl-C61-butyric acid methyl ester (PCBM) blends that are commensurate with the size of a BHJ thin film is examined using petascale coarse-grained molecular dynamics simulations. Comparisons between two-component and three-component systems containing short P3HT chains as additives undergoing thermal annealing demonstrate that the short chains alter the morphology in apparently useful ways: they efficiently migrate to the P3HT/PCBM interface, increasing the P3HT domain size and interfacial area. Simulation results agree with depth profiles determined from neutron reflectometry measurements that reveal PCBM enrichment near substrate and air interfaces but a decrease in that PCBM enrichment when a small amount of short P3HT chains are integrated into the BHJ blend. Atomistic simulations of the P3HT/PCBM blend interfaces show a nonmonotonic dependence of the interfacial thickness as a function of number of repeat units in the oligomeric P3HT additive, and the thiophene rings orient parallel to the interfacial plane as they approach the PCBM domain. Using the nanoscale geometries of the P3HT oligomers, LUMO and HOMO energy levels calculated by density functional theory are found to be invariant across the donor/acceptor interface. These connections between additives, processing, and morphology at all length scales are generally useful for efforts to improve device performance.

It is generally believed that the strength of the polymer-nanoparticle interaction controls the modification of near-interface segmental mobility in polymer nanocomposites (PNCs). However, little is known about the effect of covalent bonding on the segmental dynamics and glass transition of matrix-free polymer-grafted nanoparticles (PGNs), especially when compared to PNCs. In this article, we directly compare the static and dynamic properties of poly(2-vinylpyridine)/silica-based nanocomposites with polymer chains either physically adsorbed (PNCs) or covalently bonded (PGNs) to identical silica nanoparticles (RNP = 12.5 nm) for three different molecular weight (MW) systems. Interestingly, when the MW of the matrix is as low as 6 kg/mol (RNP/Rg = 5.4) or as high as 140 kg/mol (RNP/Rg= 1.13), both small-angle X-ray scattering and broadband dielectric spectroscopy show similar static and dynamic properties for PNCs and PGNs. However, for the intermediate MW of 18 kg/mol (RNP/Rg = 3.16), the difference between physical adsorption and covalent bonding can be clearly identified in the static and dynamic properties of the interfacial layer. We ascribe the differences in the interfacial properties of PNCs and PGNs to changes in chain stretching, as quantified by self-consistent field theory calculations. These results demonstrate that the dynamic suppression at the interface is affected by the chain stretching; that is, it depends on the anisotropy of the segmental conformations, more so than the strength of the interaction, which suggests that the interfacial dynamics can be effectively tuned by the degree of stretching-a parameter accessible from the MW or grafting density.

Enabling memristive, neuromorphic, and quantum-based computing as well as efficient mainstream energy storage and conversion technologies requires the next generation of materials customized at the atomic scale. This requires full control of atomic arrangement and bonding in three dimensions. The last two decades witnessed substantial industrial, academic, and government research efforts directed toward this goal through various lithographies and scanning-probe-based methods. These technologies emphasize 2D surface structures, with some limited 3D capability. Recently, a range of focused electron- and ion-based methods have demonstrated compelling alternative pathways to achieving atomically precise manufacturing of 3D structures in solids, liquids, and at interfaces. Electron and ion microscopies offer a platform that can simultaneously observe dynamic and static structures at the nano- and atomic scales and also induce structural rearrangements and chemical transformation. The addition of predictive modeling or rapid image analytics and feedback enables guiding these in a controlled manner. Here, we review the recent results that used focused electron and ion beams to create free-standing nanoscale 3D structures, radiolysis, and the fabrication potential with liquid precursors, epitaxial crystallization of amorphous oxides with atomic layer precision, as well as visualization and control of individual dopant motion within a 3D crystal lattice. These works lay the foundation for approaches to directing nanoscale level architectures and offer a potential roadmap to full 3D atomic control in materials. In this paper, we lay out the gaps that currently constrain the processing range of these platforms, reflect on indirect requirements, such as the integration of large-scale data analysis with theory, and discuss future prospects of these technologies.

Orthorhombic black phosphorus (BP) and other layered materials, such as gallium telluride (GaTe) and tin selenide (SnSe), stand out among two-dimensional (2D) materials owing to their anisotropic in-plane structure. This anisotropy adds a new dimension to the properties of 2D materials and stimulates the development of angle-resolved photonics and electronics. However, understanding the effect of anisotropy has remained unsatisfactory to date, as shown by a number of inconsistencies in the recent literature. We use angle-resolved absorption and Raman spectroscopies to investigate the role of anisotropy on the electron-photon and electron-phonon interactions in BP. We highlight, both experimentally and theoretically, a nontrivial dependence between anisotropy and flake thickness and photon and phonon energies. We show that once understood, the anisotropic optical absorption appears to be a reliable and simple way to identify the crystalline orientation of BP, which cannot be determined from Raman spectroscopy without the explicit consideration of excitation wavelength and flake thickness, as commonly used previously.

The structure and dynamics of a model branched polymer was investigated through molecular dynamics simulations and neutron scattering experiments. The polymer confinement, monomer concentration, and solvent quality were varied in the simulations and detailed comparisons between the calculated structural and dynamical properties of the unconfined polymer and those confined within an adsorbing and nonadsorbing cylindrical pore, representing the silica based structural support of the composite, were made. The simulations show a direct relationship in the structure of the polymer and the nonmonotonic dynamics as a function of monomer concentration within an adsorbing cylindrical pore. However, the nonmonotonic behavior disappears for the case of the branched polymer within a nonadsorbing cylindrical pore. Overall, the simulation results are in good agreement with quasi-elastic neutron scattering (QENS) studies of branched poly(ethylenimine) in mesoporous silica (SBA-15) of comparable size, suggesting an approach that can be a useful guide for understanding how to tune porous polymer composites for enhancing desired dynamical and structural behavior targeting carbon dioxide adsorption.

Unique twisted bilayers of MoSe2 with multiple stacking orientations and interlayer couplings in the narrow range of twist angles, 60 ± 3°, are revealed by low-frequency Raman spectroscopy and theoretical analysis. The slight deviation from 60° allows the concomitant presence of patches featuring all three high-symmetry stacking configurations (2H or AA', AB', and A'B) in one unique bilayer system. In this case, the periodic arrangement of the patches and their size strongly depend on the twist angle. Ab initio modeling predicts significant changes in frequencies and intensities of low-frequency modes versus stacking and twist angle. Experimentally, the variable stacking and coupling across the interface are revealed by the appearance of two breathing modes, corresponding to the mixture of the high-symmetry stacking configurations and unaligned regions of monolayers. Only one breathing mode is observed outside the narrow range of twist angles. This indicates a stacking transition to unaligned monolayers with mismatched atom registry without the in-plane restoring force required to generate a shear mode. The variable interlayer coupling and spacing in transition metal dichalcogenide bilayers revealed in this study may provide an interesting platform for optoelectronic applications of these materials.

van der Waals homo- and heterostructures assembled by stamping monolayers together present optoelectronic properties suitable for diverse applications. Understanding the details of the interlayer stacking and resulting coupling is crucial for tuning these properties. We investigated the low-frequency interlayer shear and breathing Raman modes (<50 cm(-1)) in twisted bilayer MoS2 by Raman spectroscopy and first-principles modeling. Twisting significantly alters the interlayer stacking and coupling, leading to notable frequency and intensity changes of low-frequency modes. The frequency variation can be up to 8 cm(-1) and the intensity can vary by a factor of ∼5 for twisting angles near 0° and 60°, where the stacking is a mixture of high-symmetry stacking patterns and is thus sensitive to twisting. For twisting angles between 20° and 40°, the interlayer coupling is nearly constant because the stacking results in mismatched lattices over the entire sample. It follows that the Raman signature is relatively uniform. Note that for some samples, multiple breathing mode peaks appear, indicating nonuniform coupling across the interface. In contrast to the low-frequency interlayer modes, high-frequency intralayer Raman modes are much less sensitive to interlayer stacking and coupling. This research demonstrates the effectiveness of low-frequency Raman modes for probing the interfacial coupling and environment of twisted bilayer MoS2 and potentially other two-dimensional materials and heterostructures.

Molecular simulations of the diffusion of EMIM(+) and TFSI(-) ions in slit-shaped micropores under conditions similar to those during charging show that in pores that accommodate only a single layer of ions, ions diffuse increasingly faster as the pore becomes charged (with diffusion coefficients even reaching ∼5 × 10(-9) m(2)/s), unless the pore becomes very highly charged. In pores wide enough to fit more than one layer of ions, ion diffusion is slower than in the bulk and changes modestly as the pore becomes charged. Analysis of these results revealed that the fast (or slow) diffusion of ions inside a micropore during charging is correlated most strongly with the dense (or loose) ion packing inside the pore. The molecular details of the ions and the precise width of the pores modify these trends weakly, except when the pore is so narrow that the ion conformation relaxation is strongly constrained by the pore walls.

Self-assembly of a nucleoside on Au(111) was studied to ascertain whether polymerization on well-defined substrates constitutes a promising approach for making sequence-controlled polymers. Scanning tunneling microscopy and density functional theory were used to investigate the self-assembly on Au(111) of (RS)-N(9)-(2,3-dihydroxypropyl)adenine (DHPA), a plausibly prebiotic nucleoside analog of adenosine. It is found that DHPA molecules self-assemble into a hydrogen-bonded polymer that grows almost exclusively along the herringbone reconstruction pattern, has a two component sequence that is repeated over hundreds of nanometers, and is erasable with electron-induced excitation. Although the sequence is simple, more complicated ones are envisioned if two or more nucleoside types are combined. Because polymerization occurs on a substrate in a dry environment, the success of each combination can be gauged with high-resolution imaging and accurate modeling techniques. These characteristics make nucleoside self-assembly on a substrate an attractive approach for designing sequence-controlled polymers. Further, by choosing plausibly prebiotic nucleosides, insights may be provided into how nature created the first sequence-controlled polymers capable of storing information. Such insights, in turn, can inspire new ways of synthesizing sequence-controlled polymers.

The isolation of graphene in 2004 from graphite was a defining moment for the "birth" of a field: two-dimensional (2D) materials. In recent years, there has been a rapidly increasing number of papers focusing on non-graphene layered materials, including transition-metal dichalcogenides (TMDs), because of the new properties and applications that emerge upon 2D confinement. Here, we review significant recent advances and important new developments in 2D materials "beyond graphene". We provide insight into the theoretical modeling and understanding of the van der Waals (vdW) forces that hold together the 2D layers in bulk solids, as well as their excitonic properties and growth morphologies. Additionally, we highlight recent breakthroughs in TMD synthesis and characterization and discuss the newest families of 2D materials, including monoelement 2D materials (i.e., silicene, phosphorene, etc.) and transition metal carbide- and carbon nitride-based MXenes. We then discuss the doping and functionalization of 2D materials beyond graphene that enable device applications, followed by advances in electronic, optoelectronic, and magnetic devices and theory. Finally, we provide perspectives on the future of 2D materials beyond graphene.

Surface-enhanced Raman scattering (SERS) on two-dimensional (2D) layered materials has provided a unique platform to study the chemical mechanism (CM) of the enhancement due to its natural separation from electromagnetic enhancement. The CM stems from the charge interactions between the substrate and molecules. Despite the extensive studies of the energy alignment between 2D materials and molecules, an understanding of how the electronic properties of the substrate are explicitly involved in the charge interaction is still unclear. Lately, a new group of 2D layered materials with anisotropic structures, including orthorhombic black phosphorus (BP) and triclinic rhenium disulfide (ReS2), has attracted great interest due to their unique anisotropic electrical and optical properties. Herein, we report a unique anisotropic Raman enhancement on few-layered BP and ReS2 using copper phthalocyanine (CuPc) molecules as a Raman probe, which is absent on isotropic graphene and h-BN. According to detailed Raman tensor analysis and density functional theory calculations, anisotropic charge interactions between the 2D materials and molecules are responsible for the angular dependence of the Raman enhancement. Our findings not only provide new insights into the CM process in SERS, but also open up new avenues for the exploration and application of the electronic properties of anisotropic 2D layered materials.

The atomic-level sculpting of 3D crystalline oxide nanostructures from metastable amorphous films in a scanning transmission electron microscope (STEM) is demonstrated. Strontium titanate nanostructures grow epitaxially from the crystalline substrate following the beam path. This method can be used for fabricating crystalline structures as small as 1-2 nm and the process can be observed in situ with atomic resolution. The fabrication of arbitrary shape structures via control of the position and scan speed of the electron beam is further demonstrated. Combined with broad availability of the atomic resolved electron microscopy platforms, these observations suggest the feasibility of large scale implementation of bulk atomic-level fabrication as a new enabling tool of nanoscience and technology, providing a bottom-up, atomic-level complement to 3D printing., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The detailed nature of spatially heterogeneous dynamics of glycerol-silica nanocomposites is unraveled by combining dielectric spectroscopy with atomistic simulation and statistical mechanical theory. Analysis of the spatial mobility gradient shows no "glassy" layer, but the α-relaxation time near the nanoparticle grows with cooling faster than the α-relaxation time in the bulk and is ∼20 times longer at low temperatures. The interfacial layer thickness increases from ∼1.8 nm at higher temperatures to ∼3.5 nm upon cooling to near bulk Tg. A real space microscopic description of the mobility gradient is constructed by synergistically combining high temperature atomistic simulation with theory. Our analysis suggests that the interfacial slowing down arises mainly due to an increase of the local cage scale barrier for activated hopping induced by enhanced packing and densification near the nanoparticle surface. The theory is employed to predict how local surface densification can be manipulated to control layer dynamics and shear rigidity over a wide temperature range.

Rationally encoding molecular interactions that can control the assembly structure and functional expression in a solution of conjugated polymers hold great potential for enabling optimal organic optoelectronic and sensory materials. In this work, we show that thermally-controlled and surfactant-guided assembly of water-soluble conjugated polymers in aqueous solution is a simple and effective strategy to generate optoelectronic materials with the desired molecular ordering. We have studied a conjugated polymer consisting of a hydrophobic thiophene backbone and hydrophilic, thermo-responsive ethylene oxide side groups, which shows a step-wise, multi-dimensional assembly in water. By incorporating the polymer into phase-segregated domains of an amphiphilic surfactant in solution, we demonstrate that both chain conformation and degree of molecular ordering of the conjugated polymer can be tuned in hexagonal, micellar and lamellar phases of the surfactant solution. The controlled molecular ordering in conjugated polymer assembly is demonstrated as a key factor determining the electronic interaction and optical function.

Harnessing big data, deep data, and smart data from state-of-the-art imaging might accelerate the design and realization of advanced functional materials. Here we discuss new opportunities in materials design enabled by the availability of big data in imaging and data analytics approaches, including their limitations, in material systems of practical interest. We specifically focus on how these tools might help realize new discoveries in a timely manner. Such methodologies are particularly appropriate to explore in light of continued improvements in atomistic imaging, modelling and data analytics methods.

In both research and industry, materials will be exposed to stresses, be it during fabrication, normal use, or mechanical failure. The response to external stress will have an important impact on properties, especially when atomic details govern the functionalities of the materials. This review aims at summarizing current research involving the responses of graphene and graphene materials to applied stress at the nanoscale, and to categorize them by stress-strain behavior. In particular, we consider the reversible functionalization of graphene and graphene materials by way of elastic deformation and strain engineering, the plastic deformation of graphene oxide and the emergence of such in normally brittle graphene, the formation of defects as a response to stress under high temperature annealing or irradiation conditions, and the properties that affect how, and mechanisms by which, pristine, defective, and polycrystalline graphene fail catastrophically during fracture. Overall we find that there is significant potential for the use of existing knowledge, especially that of strain engineering, as well as potential for additional research into the fracture mechanics of polycrystalline graphene and device functionalization by way of controllable plastic deformation of graphene.

Developing practical approaches to effectively reduce the amount of deep defect levels in semiconductors is critical for their use in electronic and optoelectronic devices, but this still remains a very challenging task. In this Letter, we propose that specific alloying can provide an effective means to suppress the deep defect levels in semiconductors while maintaining their basic electronic properties. Specifically, we demonstrate that for transition-metal dichalcogenides, such as MoSe&#95;{2} and WSe&#95;{2}, where anion vacancies are the most abundant defects that can induce deep levels, the deep levels can be effectively suppressed in Mo&#95;{1-x}W&#95;{x}Se&#95;{2} alloys at low W concentrations. This surprising phenomenon is associated with the fact that the band edge energies can be substantially tuned by the global alloy concentration, whereas the defect level is controlled locally by the preferred locations of Se vacancies around W atoms. Our findings illustrate a concept of alloy engineering and provide a promising approach to control the defect properties of semiconductors.

Structural characteristics of the active layers in organic photovoltaic (OPV) devices play a critical role in charge generation, separation and transport. Here we report on morphology and structural control of p-DTS(FBTTh2)2:PC71BM films by means of thermal annealing and 1,8-diiodooctane (DIO) solvent additive processing, and correlate it to the device performance. By combining surface imaging with nanoscale depth-sensitive neutron reflectometry (NR) and X-ray diffraction, three-dimensional morphologies of the films are reconstituted with information extending length scales from nanometers to microns. DIO promotes the formation of a well-mixed donor-acceptor vertical phase morphology with a large population of small p-DTS(FBTTh2)2 nanocrystals arranged in an elongated domain network of the film, thereby enhancing the device performance. In contrast, films without DIO exhibit three-sublayer vertical phase morphology with phase separation in agglomerated domains. Our findings are supported by thermodynamic description based on the Flory-Huggins theory with quantitative evaluation of pairwise interaction parameters that explain the morphological changes resulting from thermal and solvent treatments. Our study reveals that vertical phase morphology of small-molecule based OPVs is significantly different from polymer-based systems. The significant enhancement of morphology and information obtained from theoretical modeling may aid in developing an optimized morphology to enhance device performance for OPVs.

Two-dimensional (2D) van der Waals (vdW) heterostructures are a family of artificially structured materials that promise tunable optoelectronic properties for devices with enhanced functionalities. Compared to transferring, direct epitaxy of vdW heterostructures is ideal for clean interlayer interfaces and scalable device fabrication. Here we report the synthesis and preferred orientations of 2D GaSe atomic layers on graphene (Gr) by vdW epitaxy. GaSe crystals are found to nucleate predominantly on random wrinkles or grain boundaries of graphene, share a preferred lattice orientation with underlying graphene, and grow into large (tens of micrometers) irregularly shaped, single-crystalline domains. The domains are found to propagate with triangular edges that merge into the large single crystals during growth. Electron diffraction reveals that approximately 50% of the GaSe domains are oriented with a 10.5 ± 0.3° interlayer rotation with respect to the underlying graphene. Theoretical investigations of interlayer energetics reveal that a 10.9° interlayer rotation is the most energetically preferred vdW heterostructure. In addition, strong charge transfer in these GaSe/Gr vdW heterostructures is predicted, which agrees with the observed enhancement in the Raman E(2)1g band of monolayer GaSe and highly quenched photoluminescence compared to GaSe/SiO2. Despite the very large lattice mismatch of GaSe/Gr through vdW epitaxy, the predominant orientation control and convergent formation of large single-crystal flakes demonstrated here is promising for the scalable synthesis of large-area vdW heterostructures for the development of new optical and optoelectronic devices.

The neighboring layers in bilayer (and few-layer) graphenes of both AA and AB stacking motifs are known to be separated at a distance corresponding to van der Waals (vdW) interactions. In this Letter, we present for the first time a new aspect of graphene chemistry in terms of a special chemical bonding between the giant graphene "molecules". Through rigorous theoretical calculations, we demonstrate that the N-doped graphenes (NGPs) with various doping levels can form an unusual two-dimensional (2D) π-π bonding in bilayer NGPs bringing the neighboring NGPs to significantly reduced interlayer separations. The interlayer binding energies can be enhanced by up to 50% compared to the pristine graphene bilayers that are characterized by only vdW interactions. Such an unusual chemical bonding arises from the π-π overlap across the vdW gap while the individual layers maintain their in-plane π-conjugation and are accordingly planar. The existence of the resulting interlayer covalent-like bonding is corroborated by electronic structure calculations and crystal orbital overlap population (COOP) analyses. In NGP-based graphite with the optimal doping level, the NGP layers are uniformly stacked and the 3D bulk exhibits metallic characteristics both in the in-plane and along the stacking directions.

Bilayer graphene (BLG) with a tunable band gap appears interesting as an alternative to graphene for practical applications; thus, its transport properties are being actively pursued. Using density functional theory and perturbation analysis, we investigated, under an external electric field, the electronic properties of BLG in various stackings relevant to recently observed complex structures. We established the first phase diagram summarizing the stacking-dependent gap openings of BLG for a given field. We further identified high-density midgap states, localized on grain boundaries, even under a strong field, which can considerably reduce the overall transport gap.

Nanocomposites based on graphene dispersed in matrices of soft materials are promising thermal management materials. Their effective thermal conductivity depends on both the thermal conductivity of graphene and the conductance of the thermal transport across graphene-matrix interfaces. Here, we report on molecular dynamics simulations of the thermal transport across the interfaces between defected graphene and soft materials in two different modes: in the "across" mode, heat enters graphene from one side of its basal plane and leaves through the other side; in the "non-across" mode, heat enters or leaves graphene simultaneously from both sides of its basal plane. We show that as the density of vacancy defects in graphene increases from 0% to 8%, the conductance of the interfacial thermal transport in the "across" mode increases from 160.4 ± 16 to 207.8 ± 11 MW/m(2) K, while that in the "non-across" mode increases from 7.2 ± 0.1 to 17.8 ± 0.6 MW/m(2) K. The molecular mechanisms for these variations of thermal conductance are clarified using the phonon density of states and structural characteristics of defected graphene. On the basis of these results and effective medium theory, we show that it is possible to enhance the effective thermal conductivity of thermal nanocomposites by tuning the density of vacancy defects in graphene despite the fact that graphene's thermal conductivity always decreases as vacancy defects are introduced.

The tunable optoelectronic properties of stacked two-dimensional (2D) crystal monolayers are determined by their stacking orientation, order, and atomic registry. Atomic-resolution Z-contrast scanning transmission electron microscopy (AR-Z-STEM) and electron energy loss spectroscopy (EELS) can be used to determine the exact atomic registration between different layers, in few-layer 2D stacks; however, fast optical characterization techniques are essential for rapid development of the field. Here, using two- and three-layer MoSe2 and WSe2 crystals synthesized by chemical vapor deposition, we show that the generally unexplored low frequency (LF) Raman modes (<50 cm(-1)) that originate from interlayer vibrations can serve as fingerprints to characterize not only the number of layers, but also their stacking configurations. Ab initio calculations and group theory analysis corroborate the experimental assignments determined by AR-Z-STEM and show that the calculated LF mode fingerprints are related to the 2D crystal symmetries.

Chemical bonding interactions are the main driving force for the formation of molecules and materials from atoms. The two-electron/multicenter pancake π-π bonding found in phenalenyl (PLY, 1) radical π-dimers is intriguing due to its unconventional nature of covalent bonding for molecular aggregations and its propensity to induce unique optical, electronic, and magnetic properties. By using high-level quantum chemistry calculations, we show that the B- or N-doped PLYs (2 and 4), usually considered closed-shell and therefore trifling, can be rendered open-shell singlet by proper edge substitutions (3 and 5). The resulting two unpaired valence electrons on each molecular unit contribute to the formation of a genuine pancake-shaped 4e/all-sites double π-π bonding upon intermolecular π-dimerization, in contrast to the 2e/half-sites single π-π bonding in the parent PLY π-dimers. The unusual double π-π bonding motif discovered in these PLY analogues may broaden the landscape of, and find new applications for, intermolecular covalent bonding interactions.

As a new two-dimensional layered material, black phosphorus (BP) is a very promising material for nanoelectronics and optoelectronics. We use Raman spectroscopy and first-principles theory to characterize and understand the low-frequency (LF) interlayer breathing modes (<100 cm(-1)) in few-layer BP for the first time. Using a laser polarization dependence study and group theory analysis, the breathing modes are assigned to Ag symmetry. Compared to the high-frequency (HF) Raman modes, the LF breathing modes are considerably more sensitive to interlayer coupling and, thus, their frequencies show a stronger dependence on the number of layers. Hence, they constitute an effective means to probe both the crystalline orientation and thickness of few-layer BP. Furthermore, the temperature dependence shows that in the temperature range -150 to 30 °C, the breathing modes have a weak anharmonic behavior, in contrast to the HF Raman modes that exhibit strong anharmonicity.

Solid polymer electrolytes, such as polyethylene oxide (PEO) based systems, have the potential to replace liquid electrolytes in secondary lithium batteries with flexible, safe, and mechanically robust designs. Previously reported PEO nanocomposite electrolytes routinely use metal oxide nanoparticles that are often 5-10 nm in diameter or larger. The mechanism of those oxide particle-based polymer nanocomposite electrolytes is under debate and the ion transport performance of these systems is still to be improved. Herein we report a 6-fold ion conductivity enhancement in PEO/lithium bis(trifluoromethanesulfonyl) imide (LiTFSI)-based solid electrolytes upon the addition of fullerene derivatives. The observed conductivity improvement correlates with nanometer-scale fullerene crystallite formation, reduced crystallinities of both the (PEO)6:LiTFSI phase and pure PEO, as well as a significantly larger PEO free volume. This improved performance is further interpreted by enhanced decoupling between ion transport and polymer segmental motion, as well as optimized permittivity and conductivity in bulk and grain boundaries. This study suggests that nanoparticle induced morphological changes, in a system with fullerene nanoparticles and no Lewis acidic sites, play critical roles in their ion conductivity enhancement. The marriage of fullerene derivatives and solid polymer electrolytes opens up significant opportunities in designing next-generation solid polymer electrolytes with improved performance.

Modulation of intermolecular interactions in response to external electric fields could be fundamental to the formation of unusual forms of water, such as water whiskers. However, a detailed understanding of the nature of intermolecular interactions in such systems is lacking. In this paper, we present novel theoretical results based on electron correlation calculations regarding the nature of H-bonds in water whiskers, which is revealed by studying their evolution under external electric fields with various field strengths. We find that the water whiskers consisting of 2-7 water molecules all have a chain-length dependent critical electric field. Under the critical electric field, the most compact chain structures are obtained, featuring very strong H-bonds, herein referred to as covalent H-bonds. In the case of a water dimer whisker, the bond length of the novel covalent H-bond shortens by 25%, the covalent bond order increases by 9 times, and accordingly the H-bond energy is strengthened by 5 times compared to the normal H-bond in a (H2O)2 cluster. Below the critical electric field, it is observed that, with increasing field strength, H-bonding orbitals display gradual evolutions in the orbital energy, orbital ordering, and orbital nature (i.e., from typical π-style orbital to unusual σ-style double H-bonding orbital). We also show that, beyond the critical electric field, a single water whisker may disintegrate to form a loosely bound zwitterionic chain due to a relay-style proton transfer, whereas two water whiskers may undergo intermolecular cross-linking to form a quasi-two-dimensional water network. Overall, these results help shed new insight on the effects of electric fields on water whisker formation.

Molecular dynamics simulations are used to understand the self-assembly and structural relaxation in ionomer melts containing less than 10% degree of ionization on the backbone. The self-assembly of charged sites and counterions shows structural ordering and agglomeration with a range of structures that can be achieved by changing the dielectric constant of the medium. The intermediate scattering function shows a decoupling of charge and counterion relaxation at longer length scales for only high dielectric constant and at shorter length scales for all dielectric constants. Overall, the slow structural decay of counterions in the strongly correlated ionomer system closely resembles transport properties of semi-flexible polymers.

The recently synthesized freestanding four-atom-thick double-layer sheet of ZnSe holds great promise as an ultraflexible and transparent photoelectrode material for solar water splitting. In this work, we report theoretical studies on a novel three-atom-thick single-layer sheet of ZnSe that demonstrates a strong quantum confinement effect by exhibiting a large enhancement of the band gap (2.0 eV) relative to the zinc blende (ZB) bulk phase. Theoretical optical absorbance shows that the largest absorption of this ultrathin single-layer sheet of ZnSe occurs at a wavelength similar to its four-atom-thick double-layer counterpart, suggesting a comparable behavior on incident photon-to-current conversion efficiency for solar water splitting, among a wealth of potential applications. The results presented herein for ZnSe may be generalized to other group II-VI analogues.

Frictional forces arise whenever objects around us are set in motion. Controlling them in a rational manner means gaining leverage over mechanical energy losses and wear. This paper presents a way of manipulating nanoscale friction by means of in situ lubrication and interfacial electrochemistry. Water lubricant is directionally condensed from the vapor phase at a moving metal-ionic crystal interface by a strong confined electric field, thereby allowing friction to be tuned up or down via an applied bias. The electric potential polarity and ionic solid solubility are shown to strongly influence friction between the atomic force microscope (AFM) tip and salt surface. An increase in friction is associated with the AFM tip digging into the surface, whereas reducing friction does not influence its topography. No current flows during friction variation, which excludes Joule heating and associated electrical energy losses. The demonstrated novel effect can be of significant technological importance for controlling friction in nano- and micro-electromechanical systems.

Polymerized ionic liquids (PolyILs) are promising materials for various solid state electronic applications such as dye-sensitized solar cells, lithium batteries, actuators, field-effect transistors, light emitting electrochemical cells, and electrochromic devices. However, fundamental understanding of interconnection between ionic transport and mechanical properties in PolyILs is far from complete. In this work, local charge transport and structural changes in films of a PolyIL are studied using an integrated experiment-theory based approach. Experimental data for the kinetics of charging and steady state current-voltage relations can be explained by taking into account the dissociation of ions under an applied electric field (known as the Wien effect). Onsager's theory of the Wien effect coupled with the Poisson-Nernst-Planck formalism for the charge transport is found to be in excellent agreement with the experimental results. The agreement between the theory and experiments allows us to predict structural properties of the PolyIL films. We have observed significant softening of the PolyIL films beyond certain threshold voltages and formation of holes under a scanning probe microscopy (SPM) tip, through which an electric field was applied. The observed softening is explained by the theory of depression in glass transition temperature resulting from enhanced dissociation of ions with an increase in applied electric field.

Understanding the dynamic charge storage in nanoporous electrodes with room-temperature ionic liquid electrolytes is essential for optimizing them to achieve supercapacitors with high energy and power densities. Herein, we report coarse-grained molecular dynamics simulations of the cyclic voltammetry of supercapacitors featuring subnanometer pores and model ionic liquids. We show that the cyclic charging and discharging of nanopores are governed by the interplay between the external field-driven ion transport and the sloshing dynamics of ions inside of the pore. The ion occupancy along the pore length depends strongly on the scan rate and varies cyclically during charging/discharging. Unlike that at equilibrium conditions or low scan rates, charge storage at high scan rates is dominated by counterions while the contribution by co-ions is marginal or negative. These observations help explain the perm-selective charge storage observed experimentally. We clarify the mechanisms underlying these dynamic phenomena and quantify their effects on the efficiency of the dynamic charge storage in nanopores.

Thermal transport across the interfaces between few-layer graphene sheets and soft materials exhibits intriguing anomalies when interpreted using the classical Kapitza model, e.g. the conductance of the same interface differs greatly for different modes of interfacial thermal transport. Using atomistic simulations, we show that such thermal transport follows a nonlocal flux-temperature drop constitutive law and is characterized jointly by a quasi-local conductance and a nonlocal conductance instead of the classical Kapitza conductance. The nonlocal model enables rationalization of many anomalies of the thermal transport across embedded few-layer graphene sheets and should be used in studies of interfacial thermal transport involving few-layer graphene sheets or other ultra-thin layered materials.

CONSPECTUS: The vital importance of energy to society continues to demand a relentless pursuit of energy responsive materials that can bridge fundamental chemical structures at the molecular level and achieve improved functionality and performance. This demand can potentially be realized by harnessing the power of self-assembly, a spontaneous process where molecules or much larger entities form ordered aggregates as a consequence of predominately noncovalent (weak) interactions. Self-assembly is the key to bottom-up design of molecular devices, because the nearly atomic-level control is very difficult to realize in a top-down, for example, lithographic, approach. However, while function in simple systems such as single crystals can often be evaluated a priori, predicting the function of the great variety of self-assembled molecular architectures is complicated by the lack of understanding and control over nanoscale interactions, mesoscale architectures, and macroscale order. To establish a foundation toward delivering practical solutions, it is critical to develop an understanding of the chemical and physical mechanisms responsible for the self-assembly of molecular and hybrid materials on various support substrates. Typical molecular self-assembly involves noncovalent intermolecular and substrate-molecule interactions. These interactions remain poorly understood, due to the combination of many-body interactions compounded by local or collective influences from the substrate atomic lattice and electronic structure. Progress toward unraveling the underlying physicochemical processes that control the structure and macroscopic physical, chemical, mechanical, electrical, and transport properties of materials increasingly requires tight integration of theory, modeling, and simulation with precision synthesis, advanced experimental characterization, and device measurements. Theory, modeling, and simulation can accelerate the process of materials understanding and design by providing atomic level understanding of the underlying physicochemical phenomena (illuminating connections between experiments). It can also provide the ability to explore new materials and conditions before they are realized in the laboratory. With tight integration and feedback with experiment, it becomes feasible to identify promising materials or processes for targeted energy applications. In this Account, we highlight recent advances and success in using an integrated approach based on electronic structure simulations and scanning probe microscopy techniques to study and design functional materials formed from the self-assembly of molecules into supramolecular or polymeric architectures on substrates.

Single-layer black phosphorus (BP), or phosphorene, is a highly anisotropic two-dimensional elemental material possessing promising semiconductor properties for flexible electronics. However, the direct bandgap of single-layer black phosphorus predicted theoretically has not been directly measured, and the properties of its edges have not been considered in detail. Here we report atomic scale electronic variation related to strain-induced anisotropic deformation of the puckered honeycomb structure of freshly cleaved black phosphorus using a high-resolution scanning tunneling spectroscopy (STS) survey along the light (x) and heavy (y) effective mass directions. Through a combination of STS measurements and first-principles calculations, a model for edge reconstruction is also determined. The reconstruction is shown to self-passivate most dangling bonds by switching the coordination number of phosphorus from 3 to 5 or 3 to 4.

Two-dimensional interfaces between crystalline materials have been shown to generate unusual interfacial electronic states in complex oxides. Recently, a one-dimensional interface has been realized in hexagonal boron nitride and graphene planar heterostructures, where a polar-on-nonpolar one-dimensional boundary is expected to possess peculiar electronic states associated with edge states of graphene and the polarity of boron nitride. Here we present a combined scanning tunnelling microscopy and first-principles theory study of the graphene-boron nitride boundary to provide a first glimpse into the spatial and energetic distributions of the one-dimensional boundary states down to atomic resolution. The revealed boundary states are about 0.6 eV below or above the Fermi level depending on the termination of the boron nitride at the boundary, and are extended along but localized at the boundary. These results suggest that unconventional physical effects similar to those observed at two-dimensional interfaces can also exist in lower dimensions.

Structural evolution from poly(lactide) (PLA) macromonomer to resultant PLA molecular bottlebrush during ring opening metathesis polymerization (ROMP) was investigated for the first time by combining size exclusion chromatography (SEC), small-angle neutron scattering (SANS), and coarse-grained molecular dynamics (CG-MD) simulations. Multiple aliquots were collected at various reaction times during ROMP and subsequently analyzed by SEC and SANS. These complementary techniques enable the understanding of systematic changes in conversion, molecular weight and dispersity as well as structural details of PLA molecular bottlebrushes. CG-MD simulation not only predicts the experimental observations, but it also provides further insight into the analysis and interpretation of data obtained in SEC and SANS experiments. We find that PLA molecular bottlebrushes undergo three conformational transitions with increasing conversion (i.e., increasing the backbone length): (1) from an elongated to a globular shape due to longer side chain at low conversion, (2) from a globular to an elongated shape at intermediate conversion caused by excluded volume of PLA side chain, and (3) the saturation of contour length at high conversion due to chain transfer reactions.

Coarse-grained molecular dynamics simulations are used to probe the dynamic phenomena of polymer melts confined in nanopores. The simulation results show excellent agreement in the values obtained for the normalized coherent single chain dynamic structure factor, S(Q,Δt)/S(Q,0). In the bulk configuration, both simulations and experiments confirm that the polymer chains follow Rouse dynamics. However, under confinement, the Rouse modes are suppressed. The mean-square radius of gyration ⟨R(g)(2)⟩ and the average relative shape anisotropy ⟨κ(2)⟩ of the conformation of the polymer chains indicate a pancake-like conformation near the surface and a bulk-like conformation near the center of the confining cylinder. This was confirmed by direct visualization of the polymer chains. Despite the presence of these different conformations, the average form factor of the confined chains still follows the Debye function which describes linear ideal chains, which is in agreement with small angle neutron scattering experiments (SANS). The experimentally inaccessible mean-square displacement (MSD) of the confined monomers, calculated as a function of radial distance from the pore surface, was obtained in the simulations. The simulations show a gradual increase of the MSD from the adsorbed, but mobile layer, to that similar to the bulk far away from the surface.

We report detailed simulation results on the formation dynamics of an electrical double layer (EDL) inside an electrochemical cell featuring room-temperature ionic liquids (RTILs) enclosed between two planar electrodes. Under relatively small charging currents, the evolution of cell potential from molecular dynamics (MD) simulations during charging can be suitably predicted by the Landau-Ginzburg-type continuum model proposed recently (Bazant et al 2011 Phys. Rev. Lett. 106 046102). Under very large charging currents, the cell potential from MD simulations shows pronounced oscillation during the initial stage of charging, a feature not captured by the continuum model. Such oscillation originates from the sequential growth of the ionic space charge layers near the electrode surface. This allows the evolution of EDLs in RTILs with time, an atomistic process difficult to visualize experimentally, to be studied by analyzing the cell potential under constant-current charging conditions. While the continuum model cannot predict the potential oscillation under such far-from-equilibrium charging conditions, it can nevertheless qualitatively capture the growth of cell potential during the later stage of charging. Improving the continuum model by introducing frequency-dependent dielectric constant and density-dependent ion diffusion coefficients may help to further extend the applicability of the model. The evolution of ion density profiles is also compared between the MD and the continuum model, showing good agreement.

The performance of graphene oxide framework (GOF) membranes for water desalination is assessed using classical molecular dynamics (MD) simulations. The coupling between water permeability and salt rejection of GOF membranes is studied as a function of linker concentration n, thickness h and applied pressure ΔP. The simulations reveal that water permeability in GOF-(n,h) membranes can be tuned from ∼5 (n = 32 and h = 6.5 nm) to 400 L cm(-2) day(-1) MPa(-1) (n = 64 and h = 2.5 nm) and follows a Cnh(-αn) law. For a given pore size (n = 16 or 32), water permeability of GOF membranes increases when the pore spacing decreases, whereas for a given pore spacing (n = 32 or 64), water permeability increases by up to two orders of magnitude when the pore size increases. Furthermore, for linker concentrations n ≤ 32, the high water permeability corresponds to a 100% salt rejection, elevating this type of GOF membrane as an ideal candidate for water desalination. Compared to experimental performance of reverse osmosis membranes, our calculations suggest that under the same conditions of applied pressure and characteristics of membranes (ΔP ∼ 10 MPa and h ∼ 100 nm), one can expect a perfect salt rejection coupled to a water permeability two orders of magnitude higher than existing technologies, i.e., from a few cL cm(-2) day(-1) MPa(-1) to a few L cm(-2) day(-1) MPa(-1).

Crystalline polymorphism of organic semiconductors is among the critical factors in determining the structure and properties of the resultant organic electronic devices. Herein we report for the first time a solvent-type-dependent polymorphism of a long fused-ring organic semiconductor and its crucial effects on charge transport. A new polymorph of 5,11-bis(triethylsilylethynyl)anthradithiophene (TES ADT) is obtained using solvent-assisted crystallization, and the crystalline polymorphism of TES ADT thin films is correlated with their measured hole mobilities. The best-performing organic thin film transistors of the two TES ADT polymorphs show subthreshold slopes close to 1 V dec(-1), and threshold voltages close to zero, indicating that the density of traps at the semiconductor-dielectric interface is negligible in these devices and the observed up to 10-fold differences in hole mobilities of devices fabricated with different solvents are largely resultant from the presence of two TES ADT polymorphs. Moreover, our results suggest that the best-performing TES ADT devices reported in the literature correspond to the new polymorph identified in this study, which involves crystallization from a weakly polar solvent (such as toluene and chloroform).