B. M.
Lowe
a,
Y.
Maekawa‡
b,
Y.
Shibuta‡
b,
T.
Sakata‡
b,
C.-K.
Skylaris‡
c and
N. G.
Green‡
*d
aInstitute for Complex Systems Simulation and the Electronics and Computer Science Department, University of Southampton, UK. E-mail: ben.lowe.uk@googlemail.com
bDepartment of Materials Engineering School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-8656, Japan. E-mail: sakata@biofet.t.u-tokyo.ac.jp
cSchool of Chemistry, University of Southampton, UK
dDepartment of Electronics and Computer Science, Nano Research Group, University of Southampton, UK. E-mail: ng2@ecs.soton.ac.uk
First published on 19th October 2016
Electronic devices are becoming increasingly used in chemical- and bio-sensing applications and therefore understanding the silica-electrolyte interface at the atomic scale is becoming increasingly important. For example, field-effect biosensors (BioFETs) operate by measuring perturbations in the electric field produced by the electrical double layer due to biomolecules binding on the surface. In this paper, explicit-solvent atomistic calculations of this electric field are presented and the structure and dynamics of the interface are investigated in different ionic strengths using molecular dynamics simulations. Novel results from simulation of the addition of DNA molecules and divalent ions are also presented, the latter of particular importance in both physiological solutions and biosensing experiments. The simulations demonstrated evidence of charge inversion, which is known to occur experimentally for divalent electrolyte systems. A strong interaction between ions and DNA phosphate groups was demonstrated in mixed electrolyte solutions, which are relevant to experimental observations of device sensitivity in the literature. The bound DNA resulted in local changes to the electric field at the surface; however, the spatial- and temporal-mean electric field showed no significant change. This result is explained by strong screening resulting from a combination of strongly polarised water and a compact layer of counterions around the DNA and silica surface. This work suggests that the saturation of the Stern layer is an important factor in determining BioFET response to increased salt concentration and provides novel insight into the interplay between ions and the EDL.
Addition of charged macromolecules to oxide surfaces results in a perturbation of the EDL which cannot be accurately described by conventional mean-field models. In this work, electrolyte and biomolecule dynamics were studied with atomistic resolution, providing a detailed description of the electric field generated at the interface. DNA was chosen as an example of a highly-charged macromolecular polyelectrolyte which is both well-characterised and has relevance to a range of biotechnology applications. In addition, divalent ions were included, which are known to have a strong influence on the structure of the EDL, important to silica dissolution processes13 and prominent in physiological solutions14 but despite this, have received surprisingly little attention in the atomistic simulation literature.
One application of this work is in improving understanding of the mechanism-of-action of a promising class of biosensors, termed biologically-sensitive field-effect transistors (BioFETs). These sensors operate by detecting changes in the electric field within the EDL as a result of biomolecule binding, as shown schematically in Fig. 1. Reliable and quantitative prediction of changes in the electric field due to biomolecule adsorption, and hence BioFET response, is currently difficult primarily due to the complexity of the EDL.
Not only has the presence of divalent ions recently been shown to increase BioFET sensitivity,15–17 but the electric field and ion-dynamics at the interface are thought to be crucial in determining field effect biosensor response.
Recent experimental work recognises the importance of ion dynamics in BioFET engineering, with deterministic information extracted from BioFET signals in the frequency domain of the response15,22,23 which has been explained as a result of adsorption–desorption noise of biomolecules24 and perturbed charge fluctuations in the EDL.25 Experiments have shown a decrease in low frequency noise with increased ionic strength due to increased screening competition between the EDL and the semiconductor device.26 Heitzinger et al. suggested a different trend for DNA-sensing, in which they calculated an increase in the standard deviation of the FET channel current with ionic concentration due to an increasingly variable orientation of the DNA.21 These studies show that addition of a biomolecule, such as DNA, can affect EDL dynamics to the extent that a response can be observed in the frequency domain that is not apparent in the time-domain, even under high ionic strength conditions.
Experiments are not able to unambiguously decouple the signal noise originating from the semiconductor device and the EDL region in the electrolyte solution. Most current EDL theories used in the BioFET engineering field are based on equilibrium, mean-field solutions of the Poisson–Boltzmann equation. These models offer the advantage of low computational cost and can be accurate for low ionic strength and low surface-charge systems. However, for BioFET systems, this is rarely the case and finite-size steric effects render the Poisson–Boltzmann equation inaccurate without modification. Modern advances in computational power have enabled the exploration of more detailed atomistic models of the structure and dynamics of the silica-water(-bio) interface27via both classical28–36 and ab initio Molecular Dynamics (MD).37–41
This work presents MD simulations of EDL structure and ion dynamics in the interfacial region to investigate how (a) increased ionic strength and (b) addition of DNA perturbs the electric field and charge density at the silica-water interface.
Divalent ions are also known to be important in the phenomenon of charge inversion, in which the first diffuse layer in the EDL contains more counterions than needed to compensate for surface charge, which is then balanced by a second co-ion layer. This phenomenon has been attributed to two (non mutually-exclusive) mechanisms. One mechanism is via ‘specific adsorption’ of ions, caused by forces such as chemical bonding or water-mediated interactions.45 The other mechanism is via many-body ion–ion correlations, in which the chemical potential near the surface is reduced due to spatial correlations between discrete ions, with the electrostatic interactions outweighing the entropic cost of forming such a highly correlated system.45,46 Atomic Force Microscopy (AFM) has measured the effects of charge inversion at the silica-water interface using trivalent and quadrivalent ions at low concentrations (≤1 μM).47 These experiments did not demonstrate charge inversion for Mg2+; however, it would be expected that divalent ions would require a higher concentration than trivalent ions. Edel and de Mello have measured the streaming current for silica nanochannels and Mg2+ counterions, observing charge inversion for concentrations exceeding approximately 400 mM.48
From the perspective of BioFET biosensing, the importance of adding divalent and multivalent salts has also been shown. Jayant et al. have recently shown a significant enhancement in DNA-hybridisation sensitivity upon addition of trace amounts of multivalent salt (Mg2+ or Co3+) to a low-concentration NaCl background buffer.15,49,50 In their work, a background of 1 mM NaCl was used for hybridisation of ssDNA, and it was found that addition of the complementary strand in 1 mM NaCl with 100 μM Mg2+ produced a 350 ± 150% increase in potential shift relative to the control. They found that this effect was only significant when the initial concentration of monovalent salt was low, which supports an ion-competition mechanism. Modeling the DNA as a membrane, they putatively assigned the signal to a combination of (a) increased ion-exclusion of multivalent ions and (b) increased DNA-condensation in the presence of multivalent ions.15
Other authors support the notion that multivalent salts can have a significant effect on BioFET signal, for example, Rica et al. have observed an increase in DNA hybridisation on addition of multivalent salt (spermidine) and observed a corresponding FET signal indicative of charge inversion at 10 μM spermidine. The signal changed polarity at higher concentrations of spermidine, which they attribute to increased screening of the charge inverted DNA molecule.16 Shul'ga et al. have reported an almost 100% increase in glucose-sensitive enzyme-FET signal on addition of 0.1 M MgCl2 which they attribute to divalent cations affecting the rate of charge transfer of the enzyme substrate oxidation.17
Despite the importance of divalent ions, a monovalent electrolyte is typically assumed in MD simulations of hydrated surface-biomolecule systems or mathematical modelling of BioFET signals. Sakata et al.28–31 have recently used MD simulations to investigate the EDL structure and dynamics for hydrated silica-water28,29 and hydrated silica-DNA systems.31 This paper extends this work to investigate the effect of divalent Mg2+ ions upon the structure and dynamics of the EDL at this technologically important interface.
![]() | ||
Fig. 2 System summary. (a) Side-on schematic of the simulation cell. (b–d) Initial configuration of the 0 mM, 200 mM and 1 M ionic strength systems with DNA. Mg2+ = purple, Na+ = green, Cl− = yellow, O = red, Si = gold. Mg-coordinated water drawn as transparent. DNA drawn as stick representation. Videos of each trajectory can be found in the ESI.† |
Solvent boxes were prepared with an initial density of 1 g cm−3 and a vertical height of approximately 73 Å. Three different solvent boxes were considered: 0 M electrolyte (salt free, corresponding to deionised water); approximately 0.2 M ionic strength electrolyte; and approximately 1 M ionic strength electrolyte, each containing a 1:
1 molar ratio of NaCl to MgCl2. The system was geometry optimised for 5000 steps and then NVT dynamics51 were performed at 300 K for 3 ns. Dynamics were performed with the Nosé–Hoover thermostat using a Q ratio of 0.01.51
For comparison, three further systems were prepared at each ionic strength, incorporating DNA neutralised with Na+. DNA was constructed and chemically bonded to the surface following the method of Maekawa et al.31 and consisted of a d(AAAAAAAAAA) decamer with a complementary base-T strand and a net charge of −19e. The DNA was superimposed onto the solvent box, the DNA was kept fixed and a 5000 step geometry optimisation was performed for the three cases. Then, water molecules within the DNA were removed and electrolyte ions from inside to outside of the DNA. A further 5000 step geometry optimisation and 100 ps of NVT dynamics were performed in order to further relax the system. Similarly to Luan et al.,52 ionic strengths were calculated using the number of electrolyte ions counted beyond those required to neutralise the silica and DNA. The volume of the liquid system (without DNA and after geometry optimisation) was used for this calculation with the result that stated ionic strengths (0.2 M/1 M) are only estimates. Systems referred to as 0 M contain no added electrolyte in the solvent box, but are electroneutral due to Na+ associated with the surface layer and DNA. A summary of each the composition of each model is given in Table 1 and images of the initial configurations of the DNA systems at varying ionic strength are shown in Fig. 2.
Element | 0 mM | 200 mM | 1000 mM | 0 mM & DNA | 200 mM & DNA | 1000 mM & DNA |
---|---|---|---|---|---|---|
O | 8902 | 8850 | 8654 | 8724 | 8653 | 8494 |
H | 13![]() |
13![]() |
12![]() |
12![]() |
12![]() |
12![]() |
Si | 1200 | 1200 | 1200 | 1200 | 1200 | 8494 |
Cl | 0 | 18 | 87 | 0 | 6 | 87 |
Mg | 0 | 6 | 29 | 0 | 6 | 29 |
Na | 0 (33) | 6 (33) | 29 (33) | 0 (33)[19] | 6 (33)[19] | 29 (33)[19] |
In the simulations, the COMPASS II 1.2 forcefield was used. This forcefield has been parameterised predominately using ab initio data for a wide range of condensed systems; both organic and inorganic systems and on a range of ionic liquids.53,54 In this forcefield, the charges are the same as in the COMPASS I forcefield (see for reference the previous work by Maekawa et al.28,30,31), with the exception of the phosphate group on DNA molecules. For this functional group, the COMPASS II defaults were used; −0.3 for the sugar-linking oxygens, −0.822 for the non-linking oxygens (O and –O−) and +0.9246 for phosphorus atoms resulting in a total charge of −1 for each [RCH2PO4C(H)R]− substructure. Ewald summation was used for the electrostatic interactions with a 4.184 J mol threshold and an atom-based summation with a 12.5 Å cutoff for the van der Waals interactions. Unless otherwise specified, all analysis was performed over the mean of the last 1 ns with 1 ps windows.
The diffusion coefficient, D, was calculated using the Einstein relation (1), from the Mean Squared Displacement (MSD) as a function of time, t:
![]() | (1) |
The residence times of molecules coordinated with ions was calculated as per the definition of Impey et al.,55 based on the rate of decay of the time-correlation function with the parameter t* = 0.
The electric field was calculated in this work from the electrostatic forces on a test charge evaluated using two methods: coulomb summation and Ewald summation. In the Ewald sum, the calculation is periodic in x, y and z, whereas in the coulomb summation method, a finite supercell is used. Details of these methodologies can be found in the ESI 2,† each of which have advantages and disadvantages. Briefly, the coulombic sum method has error from truncation of the full periodicity of the system to a finite sum, and the Ewald sum has error from the contribution of periodic interactions in the z-direction and from the introduction of a non-uniform compensating background charge for orthorhombic unit-cells.56
Fig. 3 shows a plot of the coordination number for each ion with respect to pure water (oxygen atoms), as well as snapshots of the equilibrium solvation sphere around the (a) Mg2+ and (b) Na+ cations. The calculated coordination number of Mg2+ was approximately 6, in good agreement with neutron scattering data,57 whereas for Na+ the coordination number was approximately 5 which is in agreement with neutron scattering of 4.9 ± 158 and ab initio MD of 4.6.59 The Mg2+ cation showed a more ordered solvation sphere than Na+, demonstrated by the initial steep rise and flat region of the curve for the first compared to the second. The secondary peak in the curve shows evidence of a structured secondary solvation sphere. The difference in structure between the ions is due to the smaller ionic radii and stronger coulombic attraction of Mg2+ resulting in a much tighter, more ordered solvation sphere.
Another test parameter for the forcefield is the diffusion coefficient; the values for Cl−, Na+, Mg2+ and H2O are presented in Table 2. The calculated diffusion coefficient for water and ions overestimates the experimental value by a factor of approximately 2, which is consistent with other common water forcefields such as TIP3P and is a consequence of the well-known difficulty of accurately capturing water dynamics in empirical forcefields. Agreement within 2- to 3-fold is considered reasonable for diffusion coefficients.60 It is promising that the relative diffusivity of the ions is in good agreement with experiment, suggesting qualitatively correct dynamics.
System | Simulated D 1 × 10−5 cm2 s−1 | Literature D 1 × 10−5 cm2 s−1 |
---|---|---|
a Extrapolated to infinite dilution. b Simulated using TIP4P and OPLS forcefield at very low ionic strength. c Simulated using TIP3P forcefield. | ||
Na+ 0.2 M in water | 4.1 | 1.3 (expt.a)61 |
Mg2+ 1 M MgCl2 | 1.3 | 0.71 (expt.a)61 0.50 (expt.)62 0.60–0.79b![]() |
Cl− 1 M MgCl2 | 3.6 | 2.0 (expt.a)61 1.4 (expt.)62 2.4–2.6b![]() |
H2O bulk water | 5.930 | 2.3 (expt.)64 5.2–7.0c![]() |
Lastly, the calculated radial distribution function of Mg2+/Na+ ions to H/O atoms in bulk water is shown in Fig. 4. The coordination number previously shown in Fig. 3 was calculated as the integral of this RDF, assuming a fixed box-volume based on the last frame of analysis. Excellent agreement is seen between the simulated RDF peak position for Mg2+–O and Mg2+–H versus neutron scattering experiments; however, exact agreement is not expected due to the higher concentration of the experimental data.
![]() | ||
Fig. 4 Simulated (red/blue) Mg2+–O (left) and Mg2+–H (right) Radial Distribution Function (RDF) in bulk water, compared to experimental (green) data for Mg2+. For the simulations, 268 water molecules and 1 ion (Mg2+/Na+ respectively) were prepared in an approximately 20 Å box, NPT molecular dynamics were performed with 2 ns equilibration and a 8 ns production period over which the RDF was calculated. The neutron scattering data was obtained from Bruni et al.57 and is measured at 1![]() ![]() ![]() ![]() |
For the 0 mM case, the RDF demonstrated that for the silanolate groups, the mean O⋯Hwater hydrogen bond length was 1.25–2.0 Å and each silanolate was coordinated to 2–3 water molecules. The silanolate–Na+ coordination number increased with increasing ionic strength, resulting in coordination numbers of 0.95, 1.03 and 1.10 for 0 mM, 200 mM and 1000 mM respectively, as shown in Fig. 5(a). The increased concentration of cations did not affect the ionic bond length or result in multiple sodium ions per silanolate, but did result in a slight reduction in silanolate–water (hydrogen) coordination resulting in coordination numbers of 2.83, 2.67 and 2.5 for 0 mM, 200 mM and 1000 mM ionic strength respectively, as shown in Fig. 5(b). One explanation for this is that Na+ accumulation near the interface results in increased displacement of interfacial water.
The residence time of water molecules to silanolate groups was calculated to be 80 ps (based on a silanolate to hydrogen distance cutoff of 3.5 Å) and showed no dependence on ionic strength. This indicates that with increasing interfacial ionic concentration, the water dissociation rate is not significantly affected (ESI 4†).
In conclusion, increasing ionic strength reduced the equilibrium water coordination to silanolate groups, but the water–silanolate kinetics were not significantly affected.
For the duration of all simulated systems presented in this paper, the Mg(aq)2+ ions retained their octahedral water coordination sphere, essentially remaining a single hexahydrated cluster, [Mg(H2O)6]2+, therefore supporting the notion that Mg2+ does not directly coordinate to the DNA. However, NMR experiments have measured a mean water–Mg2+ residence time of 1.5 μs69 and therefore, microsecond timescale molecular dynamics simulations would be required to sample the full configuration space of the solvent shell.
As discussed, magnesium ions are particularly important in nature and biological systems. Mg2+ is known to specifically adsorb to some oxides surface with the extent being highly surface and pH dependent.70,71 In this work, only on a few occasions did the magnesium ions remain near (<10 Å) the surface. In the 1 M simulation, a Mg2+ that was initialised near the surface remained near the surface for the first 1.5 ns and then adsorbed to a specific site, as shown in Fig. 7. This was the only example observed of stable (>500 ps residence time) adsorption of Mg2+ to the silica surface, and is the result of the formation of a hydrogen bonded network between two silanolate groups. This result contrasts with the CP/MAS NMR experiments of d'Espinose de la Caillerie et al. which suggested that Mg2+ forms direct Si–O–Mg bonds to the surface;72 this disagreement may be because their system does not contain sodium ions at the surface. These results suggest that magnesium ions do not readily displace silanolate-bound sodium ions, and that Mg2+ interacts with the surface via its solvation shell.
The surface potential boundary condition was set to be equivalent to the surface charge density in the MD simulation of 0.2 C m−2 (see ESI 1 and 7† for further discussion); for example, for the same ionic concentrations as the 1 M mixed-valency MD system, a surface potential of 190 mV was used. This surface potential is in quantitative agreement with experimental measurements at the silica-water interface.75 At this surface potential, the ions have reached their bulk concentration within ∼1 nm from the surface, independent of ionic strength (ESI 7†).
The GCS has many limitations in high concentration systems. The dielectric constant of water will not be 80 near the surface due to the strong local interactions. Also, correlated motion between ions is not incorporated, which can be particularly important for systems containing divalent ions.45 A further, well-known limitation of the GCS model is that it does not describe finite-size effects and cannot describe the adsorption of ions to specific sites on the surface or ion–water interactions.
In the GCS, high surface potentials can result in extremely high concentrations <1 Å from the surface, which is physically unrealistic¶ due to steric constraints. By increasing the Stern layer thickness, the maximum concentration in the system is reduced, however it is possible to consider extensions to this models to better treat finite-size effects, for example, Kilic and Bazant73 have presented a model which replaces the Stern layer with a layer of cmax cations, or Brown et al.76 have presented a model incorporating hydration repulsion interactions which produced a Stern-like layer. These approaches share some findings in common, namely there is expected to be a high cationic concentration within ∼1 nm of the surface under these surface conditions, followed by the smooth decay into the bulk.36
In Fig. 8, the 1 M MD simulation results are compared to the GCS model. As expected, both give the result that the bulk concentration is reached within approximately 1 nm and there is a high cationic concentration within a few angstroms of the surface in a Stern-like layer. The simulations showed that, as discussed previously, strongly favourable solvation of the Mg2+ ions resulted in them being distributed roughly evenly through the solvent and not displacing the sodium ions at the surface, with only a small accumulation of ions at the surface (Fig. 8). This result is in good agreement with experimental observations for Mg2+ ions around DNA molecules in which the bound state is characteristic by almost complete hydration and free translation and rotational mobility.77 These results provide a description that cannot be obtained from the GCS formalism.
These results might help to interpret the experimental results of Jayant et al. demonstrating an increase in DNA hybridisation sensitivity upon addition of trace amounts of divalent salt to a monovalent electrolyte system.15 In their paper, this effect was modelled using a Poisson–Boltzmann model modified to include the effect of variable ion-permittivity due to a biomolecule layer. The results presented here suggest that the double-layer structure for Mg2+ containing electrolyte may not be adequately described by the Poisson–Boltzmann model. This discussion will be extended to the effect of DNA later.
The distribution of the diffuse-layer charge (≳3 Å from the surface) due to the ions can be seen more clearly by plotting the same function but excluding charges from water atoms, as seen in Fig. 9(b). With increasing ionic strength, an increase in cationic charge at the interface was observed which overcompensates the surface charge up to around ∼1.5 nm from the surface (z = 30 Å). This was due to a significant increase of Na+ and Mg2+ ions within a few angstroms of the surface, outweighing the cumulative charge from increased chloride ion density in the bulk. This effect is sometimes referred to as charge inversion and is known to occur for multivalent electrolytes near highly charged surfaces and highly charged molecules such as DNA.44 In this work, the charge inversion 1.5 nm from the surface was observed to be roughly proportional in magnitude to the ionic strength change (five-fold increase in ionic strength showed a five-fold increase in charge).
The orientation of the water with respect to the normal of the surface is shown in Fig. 9(c). With increasing ionic strength, orientational water polarisation increased. The water orientated H-down towards the negatively charged silica surface near the surface; this oxide-surface induced water polarisation is a well-known phenomenon.78,79 As expected, the water becomes isotropic as the bulk is reached.
In the high ionic strength simulations (0.2 M & 1 M), the water reoriented at 19–30 Å so as to be H-up towards chloride ions with oxygen towards the interfacial Mg2+ and Na+. MD simulations of the wet-charged interface in the literature have demonstrated water polarisation in monovalent electrolytes.28,80 However, the simulations presented here showed that for high ionic strengths, the accumulation of negative charge in the 35 Å to 40 Å region resulted in a secondary layer of H-down orientated water.
In conclusion, with increasing ionic strength, the equilibrium interfacial charge distribution was found to be primarily determined by the water structure around the silanolate and sodium ions at the surface, rather than as a result of the diffuse layer of ions. Charge accumulation and inversion were observed, however, water polarisation lowered the electrostatic energy of the system so as to produce a charge distribution (and therefore potential profile and electric field) which was independent of ionic strength.
In contrast, in the work presented here, the surface was initialised with a layer of compensating sodium ions; which represents an electroneutral system on the length-scale of the simulation domain. At high ionic strengths this is supported by both theory (Section 4.2) and experiment.81 For the low ionic strength (“0 mM”') system, this assumption may no longer hold, however the system provides a control against which to contrast the effect of increasing ionic strength. No increase in the interfacial charge density was observed with increasing the ionic strength of the bulk (Fig. 9(a) and ESI 8†); this can be explained by the fact that the surface layer was initialised saturated with a 1:
1 ratio of counterions. The presence of counterions in the Stern-like layer meant that cation accumulation (Fig. 9(b)) was necessarily weaker than the systems of Maekawa et al., and was completely compensated by water polarisation. This observation has significance for interpreting the response of field-effect transistor-sensors which demonstrate a change in signal with changing ionic strength;30 if many silanolate groups are ion-paired with cations or sterically obstructed from cation binding e.g. proteins, an amorphous surface, then an increase in ionic strength is predicted to correspond to a weaker change in surface-charge accumulation and therefore device response.
A competing factor that could influence device response that is not considered in this work is the effect of variable surface charge. Titration experiments suggest that increasing ionic strength results in an larger apparent negative surface charge, which may be due to altering the chemical equilibrium of the silanol groups.82
For the MD simulations presented in this paper, there was significant structuring of the water surrounding the DNA when compared to the simulations in absence of DNA (ESI 9†). The RDF of the 1 M DNA–electrolyte system is shown in Fig. 10, and shows that the DNA phosphate backbone attracted a structured counterion cloud in which sodium ions were associated closely with the phosphate groups at approximately 3 Å distance, and formed a secondary more diffuse layer at approximately 6 Å. Consistent with the literature,77 the magnesium ions were bound to phosphate groups via hydrogen bonding through the solvation shell; the exception was a single magnesium ion that was initialised in contact with a phosphate group.
The percentage of ions per phosphate group was calculated by inspection of the phosphate–ion coordination number shown in the inset of Fig. 10. Values of 76% and 42% for Na+ and Mg2+ respectively, were calculated following the methodology of Young et al.,91,92 in which the second inflection point in the RDF curve was used. These values are in excellent agreement with both the Manning condensation theory predictions and experiment, for non-mixed electrolytes NaCl and MgCl2 solutions at low concentrations. It can be noted that, given the divalence of magnesium ions, each phosphate group has +1.69e counterion charge within approximately 6 Å, which is evidence of charge inversion around the DNA. This is a phenomenon that is expected for multivalent systems at high concentrations46 and is not described by Manning condensation theory, although more recent revisions of the theory have attempted to incorporate these effects.93
Continuing from before, in the work of Jayant et al.15 the enhanced FET-signal was attributed to a combination of (a) increased ion-exclusion of multivalent ions and (b) increased DNA-condensation in the presence of multivalent ions. The simulations presented here show no evidence of Mg2+ ion-exclusion in the DNA layer, suggesting that the previous observed increase in signal in multivalent salt was due to other effects, such as Mg2+ induced DNA-condensation onto the surface as supported by other MD studies, which discussed Mg2+ induced DNA aggregation.44
The full set of data from the simulations of the charge, potential and electric field in the interfacial region for all six cases, with and without DNA, is shown in ESI 8† for ease of comparison. In the simulations, DNA did not produce a strong effect on the time-average charge distribution of the systems. This might be seen as a counter-intuitive result given that the DNA has a negative charge, however the mechanical flexibility of the DNA–Na+ system means that the time-averaged charge of DNA with respect to the surface normal, is expected to be small at any bulk ionic strength.
The electrostatic potential (relative to the silica substrate at 0 V) calculated from this charge distribution, demonstrates that the potential change at the surface (ΔDNAψ = ψDNA(zsurf) − ψnoDNA(zsurf)) on the addition of DNA was on the scale of millivolts, for example, at the position of the silanolate groups (zsurf = 17 Å), Δψ = −18 mV, −0.5 mV and −37 mV for 0 mM, 200 mM and 1 M systems respectively. This is of the same order of magnitude as surface potential change measurements for biomolecule-oxide systems,96,97 however these changes were highly sensitive to the choice of surface coordinate and therefore cannot be taken as accurate prediction of surface potential change due to DNA.
Due to natural thermal fluctuations, this mean potential is variable. In order to explore this variability over time due to DNA, the long-range electric field in to the EDL was calculated by measuring the electric field on a test charge 1 Å below the base of the silica. Fig. 11 shows the z-component of the electric field as a function of time for all six cases. No significant difference was found between the electric field for the 0 mM (M = −0.00604, SD = 0.00986) and the 0 mM DNA system (M = −0.00596, SD = 0.0102) based on an independent samples t-test (t(2999) = −0.306, p = 0.759). For the higher ionic strength systems, a small but statistically significant change (200 mM system: p = 5 × 10−51; 1 M system: p = 7 × 10−27) in electric field was observed upon addition of DNA ([E(DNA) − E(noDNA)] ≈ 0.002 V Å−1). If this change was a result of the DNA itself and not a result of noise, it would be expected that this response would be strongest in the 0 M system due to lowest ionic screening. Interestingly, the control 0 mM systems showed a 30–40% greater standard deviation in the electric field than the higher ionic strength systems; this indicates that bulk electrolyte may play a role in dampening transient fluctuations in the surface potential.
As discussed previously (Section 4.3.1), increasing ionic strength is expected to increase the surface potential.28 The incorporation of an unsaturated Stern-like layer (silanolate groups at the surface without ion paired cations), produced an electric field which was screened as a result of bulk electrolyte and a compensating background charge introduced by the Ewald summation.28 In order to compare these simulations with the work presented in this paper, the mean electric field was calculated (shown in the figure in ESI 10†).
In addition, for the simulations presented in this paper, there was no compensating background charge and the system was neutralised by a Stern-like layer and the bulk electrolyte, resulting in a more compact double layer and therefore a weaker electric field by comparison. No trend in electric field change upon ionic strength increase was observed, in contrast to the strong changes in electric field for the unsaturated Stern layer systems of Maekawa et al. (ΔE[1M]–[0M] ≈ 0.02 V Å−1).28 This comparison suggests that the electrolyte structure within several angstroms of highly-charged interfaces has a far more significant effect on the electric field, and therefore BioFET response, than biomolecule net-charge/orientation. This also emphasises the importance of developing atomistic models with a realistic description of the Stern layer in order to obtain quantitative atomistic prediction of surface potential.
Finally, for nanowire BioFETs, with small cross sections in the semiconducting region, and a correspondingly high sensitivity, the spatial variation of the field is crucial in understanding the response rather than a smeared out average approximation of the behaviour of an artificial one-dimensional system.3,98,99 The spatial variation in the electric field at the surface of the silica was investigated for the control system (“0 mM”) and is shown in Fig. 12. This demonstrated that the DNA is having an effect on the local electric field at the surface; however screening from the condensed sodium ions and polarised water reduced the field such that the field perturbation due to DNA was weak compared to thermal noise.
These conclusions suggest that, with a compact neutralising Stern-like layer on the silica surface and around the DNA phosphate groups, that the electric field due to addition of DNA and ionic strength changes is negligible compared to thermal noise. The most likely explanation for BioFET response is that the double layer is more diffuse than determined by these simulations. Over timescales not currently reachable within this MD model, a more diffuse layer might be formed via sodium ions in the Stern-like layer dissociating due to water- or DNA-induced displacement.
As discussed, understanding ion-dynamics at the silica-water interface is important for a range of systems. The simulations presented here demonstrate that increased ionic strength reduces the equilibrium water coordination to silanolate groups without significantly affecting the water–silanolate kinetics. Sodium ion surface-desorption kinetics required a two-step mechanism involving hypercoordination of the silanolate followed by desorption into the bulk. Sodium ions demonstrated a more kinetically stable solvation sphere at the silica-surface relative to the bulk due to the structuring of water at the interface.
Direct experimental measurement of the distribution of interfacial charge is not available36 and therefore the MD simulation results were compared to a continuum Poisson–Boltzmann model and revealed good agreement with regard to double-layer thickness and sodium ion accumulation. Mg2+ did not accumulate significantly at the interface, instead distributing more diffusely than predicted by the Poisson–Boltzmann model. This was as a result of the strong solvation of Mg2+ meaning it could not readily displace Na+ bound to the surface not described by the Poisson–Boltzmann formalism.
Fundamental understanding of the interfacial charge distribution and electric field is vital to understanding the mechanism of action of field-effect transistor (FET)-sensors. Increasing ionic strength was shown to result in charge inversion due to cation accumulation, an effect which is observed experimentally for divalent ions. The charge inversion 1.5 nm from the surface was observed to be roughly proportional in magnitude to the ionic strength change (five-fold increase in ionic strength showed a five-fold increase in charge). This suggests that charge inversion can begin at lower ionic strengths than those measured by Edel and de Mello.48
The results demonstrate that the equilibrium interfacial charge distribution was primarily determined by the water structure around the silanolate and sodium ions at the surface, rather than as a result of the diffuse layer of ions. Furthermore, the electrolyte structure within several angstroms of highly-charged interfaces had a far more significant effect on the electric field, and therefore FET-sensor response, than biomolecule net-charge/orientation. This supports the theory that the mechanism of action of FET-sensor is via modification of the surface chemistry (e.g. altering silanol/silanolate chemical equilibria) rather than the traditional picture based on directly sensing the electric field of the biomolecule.2
As discussed, modelling of the FET-sensor response due to biomolecules is inhibited by a lack of understanding of the interfacial electric field and ion distribution in the presence of biomolecules. The results showed that with DNA present, there was minimal effect due to water polarisation and strong screening by the condensed layer of electrolyte. The first calculation of the time-varying electric field for these systems was presented and, by comparison to a low ionic strength control, showed that bulk electrolyte plays a role in dampening transient fluctuations in the electric field and therefore device response.
These results have relevance to interpreting the experimental results of Jayant et al.15 which demonstrated an increase in DNA hybridisation sensitivity upon addition of trace amounts of divalent salt to a monovalent-electrolyte system, attributed to ion-exclusion from the DNA region and/or DNA aggregation. The simulations presented here showed no evidence of ion-exclusion from the DNA region and suggests that both (a) the Poisson–Boltzmann model may not be capable of accurately describing the EDL in the presence of mixed electrolyte, and (b) a mechanism other than ion-exclusion, such as DNA aggregation, explains the observed increase in response.
The results also emphasise the role of the Stern-like layer in understanding the response of FET sensors. Changes in the surface charge density (density of silanolate and condensed ions) would be expected to alter the electric field significantly, but in this work, surface charge is effectively fixed due to the long timescale of cation desorption and no chemical reactions. Future work will address this limitation by varying the surface charge on the surface of the model, which can be compared directly with existing experimental titration data.75 This will allow quantification of the extent to which the Stern-like layer modulates the electric field and therefore improve fundamental understanding of sensor response.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c6cp04101a |
‡ These authors contributed equally to this work. |
§ However, another more recently hypothesised mechanism of detection is via detecting changes in EDL dipole moment in absence of changes to the surface concentration of ionised groups (e.g. surface charge) or free charges (e.g. electrolyte ions or charged biomolecules like DNA) at the surface. This notion is supported by the experiments of Cahen et al. which showed FET response on addition of neutral organic molecules, or simply oxygen–water vapour.19 The importance of the dipole moment of the surface has been supported by the simulations of Heitzinger et al. based on the mean-field solution of the Poisson equation.20 This effect could be caused by electronic polarisation of neutral molecules at the silica-water interface. |
¶ The maximum concentration possible (cmax = a−3 where a is the ionic radius) given steric constraints is 25 M to 207 M for 2–4 Å cations respectively. |
This journal is © the Owner Societies 2017 |