the cluster of galaxies is one of the hottest , most x - ray luminous clusters known . since its discovery by @xcite , it has been the subject of intense and ongoing research @xcite . in particular ,
chandra observations by @xcite revealed the cluster to be a supersonic merger in the plane of the sky with a textbook example of a bow shock , making this cluster a unique case in which to study hydrodynamical properties of interacting systems .
the optical images show that the cluster has two distinct components , and the x - ray analysis reveals that the lower mass sub - cluster has recently exited the core of the main cluster with a relative velocity of @xmath1 .
although this relative velocity appears to be unusually large , an analysis of cosmological simulations demonstrates that it is well within the predicted range of the currently favored cosmological model @xcite .    due to its unique geometry and physical state
, this cluster is the best known system in which to test the dark matter hypothesis @xcite .
the observed offsets between the weak gravitational lensing mass peaks and the x - ray gas component give the most direct evidence for the presence of dark matter yet available . using the same observations , @xcite placed upper limits on the dark matter self - interaction cross section .
the goal of this work is to obtain a high - resolution , absolutely calibrated mass map with no assumptions on the physical properties of the underlying cluster potential . for this purpose
we use hst acs data , incorporating the gravitational lensing information from both multiple image systems ( strong lensing ) and from distortions of background sources ( weak lensing ) . the superb spatial resolution delivered by acs
both increases the number density of background sources that can be used for weak lensing and reveals new strong lensing candidates , especially around the subcluster .
the joint strong and weak lensing analysis thus benefits greatly from the acs data , enabling us to increase the spatial resolution and signal strength of the mass map in the core of the cluster .    encouraged by the success of the combined strong and weak lensing reconstruction method developed in @xcite ( hereafter ) and the results from applying the method to cluster rx  j1347.5@xmath21145 @xcite
, we proceed to apply this method to . because the strong lensing data are richer in this case ( rx
j1347.5@xmath21145 has not been observed with either the acs or wfpc2 cameras prior to this study ) , we improve the method as described in the text .
we perform the reconstruction in the following sequence . using acs ( multi - color where available ) hst images
, we identify the multiply imaged systems . having a spectroscopic redshift for one lensed system ( from @xcite )
, we then use the predictive power of simple strong lens modeling to estimate the redshifts of the other systems .
using positions and redshifts of the strongly lensed images and the shape measurements of the `` weakly '' lensing sources , we perform a combined strong and weak lensing mass reconstruction , thereby significantly improving the constraints on the mass and positions of the main cluster and the colliding subcluster of . in terms of the structure of this paper , in section  [ sec : data ] we describe the optical images used in this analysis , the basic image processing and the extraction of the strong  (  [ sec : datasl ] ) and weak  (  [ sec : datawl ] ) gravitational lensing data .
we then infer the mass distribution of cluster from these data in section  [ sec : swunited_bullet ] , following a demonstration of our methodology on suitable simulated data in section  [ sec : swunited ] .
we discuss the possible sources of error in section  [ sec : systematics ] and summarize our conclusions in section  [ sec : conclusions ] .
acs / wfc imaging of the cluster was carried out in cycle 13 ( proposal 10200 , pi jones ) on 2004 october 21 in two pointings with one and three different filters respectively .
the two pointings are centered on the main cluster and the subcluster with a small overlap between them .
the subcluster was observed in three different filters ( three orbits with f814w , one with f606w , and one with f435w ) , while the main cluster was observed only in the f606w filter ( one orbit ) .
the demands placed by the lensing analysis require special care when reducing the images .
we use the multidrizzle @xcite routine to align the images . to register the images with the astrometric accuracy needed for lensing analysis
, we determine the offsets among the images by extracting high @xmath3 objects in the individual , distortion corrected exposures .
we use sextractor @xcite and the iraf routine geomap to identify the objects and calculate the residual shifts and rotation of individual exposures , which were then fed back into multidrizzle .
we use `` square '' as the final drizzling kernel and an output pixel scale of @xmath4 ; this is smaller than the original pixel scale of the acs ccd to reduce the impact of resampling on the shape measurements .
owing to its complex structure , the strong lensing analysis of the cluster is complicated . among the multiply - imaged systems , only the giant arc on the nw side of the main cluster s cd has a measured redshift ( @xmath5 * ? ? ?
corresponding images for multiply - imaged sources were identified by matching both morphologies and surface brightnesses in each of the available acs bands ( f435w , f606w , and f814w ) and ground - based filters , bvr data from magellan and i - band from vlt @xcite .
we confirm ( based on photometry and morphology ) that the @xmath6 systems ( labeled a - f in table  [ tab : arcs ] ) previously identified by @xcite are indeed multiply imaged .
we identify four additional systems ( g - j ) in the subcluster region , where none were previously known .
all of these are identified in fig .
[ fig : arcslabel ] , and image positions ( corresponding to the peak surface brightness ) are presented in table  [ tab : arcs ] .
unfortunately , for the greater part of the main cluster , only a single band acs image is available , which makes identification of additional candidates significantly more ambiguous .    using these identifications
, we perform a parametrized strong lensing reconstruction . at this stage
we are not generating a detailed strong lensing model , rather we use this parametrized model only to predict the redshifts of the systems where spectroscopic redshifts are not available ( all systems but a ) and as the initial model guess for the subsequent strong and weak lensing reconstruction . for this reconstruction ,
we use only the image positions as constraints .
the parametrized model consists of two non - singular isothermal ellipses ( nie ) @xcite , for which the scaled surface mass density @xmath7 is given by @xmath8 @xmath9 is defined counterclockwise from due west and @xmath10 is related to the line - of - sight velocity dispersion @xmath11 through @xmath12 .
the components are centered on the southern cd galaxy of the main cluster and the brightest cluster galaxy ( bcg ) of the sub cluster respectively ( denoted with white crosses in fig .
[ fig : arcslabel ] ; it is not straightforward to identify the bcg of the main cluster , therefore we refer to the two main galaxies as the southern and northern cd ) .
we allow the scaling @xmath10 , core radii @xmath13 , ellipticities @xmath14 , and position angles @xmath15 of these to vary . throughout this paper as in @xcite with @xmath16 , @xmath17 being the axes ratio . ]
following the prescription of @xcite , we also include the 20 brightest cluster members from the f606w - band ( selected using color information from magellan data ) in the mass model .
the galaxies are modeled as non - singular isothermal spheres with a line - of - sight velocity dispersion @xmath18 and core radius @xmath19 following @xmath20 the proportionality constants are allowed to vary as well and below we quote the values of @xmath18 and @xmath21 for a fiducial galaxy with a f606w - band magnitude @xmath22 .
the best fit model for this system has values of @xmath23 for the main cluster , @xmath24 for the subcluster , and @xmath25 for the cluster members .
we also allow the redshifts of systems b - j to vary .
the resulting best fit redshifts are given in table  [ tab : arcs ] . to evaluate the angular diameter distances throughout the paper we assume the @xmath26cdm cosmology with @xmath27 , @xmath28 , and hubble constant @xmath29 .    for objects close in projection to the main cluster , the redshift estimates are fairly reliable because the well - measured arc a is located there ( we estimate an error of @xmath30 ) . however , the redshift estimates for those close to the subcomponent are less well - determined , and are somewhat degenerate with the adopted mass of the subcluster .
we postpone an extension of the methodology described in the following sections that incorporates simultaneous solution for the unknown source redshifts to future work , but here note that : 1 ) the parametrized model allows the generation of a set of source redshifts that , while model - dependent , are self - consistent ; 2 ) the predicted redshifts are consistent with the colors of the objects .
we adopt the redshifts estimated here for the combined strong and weak lensing reconstruction , and where possible estimate the systematic error introduced as a result of this simplification .
rrrr & 104.63316 & -55.941395 & + & 104.62988 & -55.943798 & + & 104.62954 & -55.941844 & + & 104.63042 & -55.941474 & + & 104.63775 & -55.941851 & + & 104.63338 & -55.945324 & + & 104.64709 & -55.943575 & + & 104.63528 & -55.951836 & + & 104.64008 & -55.950620 & + & 104.64232 & -55.948784 & + & 104.65155 & -55.956671 & + & 104.64778 & -55.958284 & + & 104.56568 & -55.939832 & + g & 104.56402 & -55.942113 & 1.3 + & 104.56417 & -55.944131 & + & 104.56293 & -55.939764 & + h & 104.56133 & -55.942430 & 1.9 + & 104.56189 & -55.947724 & + & 104.56186 & -55.946114 & + i & 104.56052 & -55.942930 & 2.1 + & 104.56141 & -55.944264 & + & 104.56909 & -55.946016 & + & 104.57025 & -55.944050 & [ tab : arcs ]      in the weak lensing analysis , we use the f435w , f606w , and f814w exposures for the subcluster , and f606w wherever multicolor acs data are not available .
we will only outline the main steps used for this analysis , full details on how the weak lensing catalogs were generated will be given in a future paper dealing with the weak lensing mass measurements of at large radius ( * ? ? ?
* in prep . ) .
we correct galaxy images for the psf anisotropy and psf smearing closely following the technique described in @xcite .
the procedure is based on the ksb algorithm @xcite , in particular we use the modified imcat implementation .
the ksb method is formally valid only in the weak lensing regime ( i.e. at radii much larger than where multiply imaged systems form and where the distortions are small ) .
however recent simulations from the step project @xcite show that this approach is also valid in the non - weak lensing regime ( i.e. with shear values of @xmath31 ) to the accuracy needed for a single cluster reconstruction .
lastly , in the vicinity of the highly elongated arcs we are dominated by the strong lensing signal , thereby minimizing the effects of signal dilution due to imperfect psf correction .
we select stars from the half - light - radius vs. magnitude diagram and fit a fifth order polynomial to their measured ellipticities as a function of their position on the drizzled image . after the correction the rms of the stellar ellipticity components @xmath32 and @xmath33 ( as defined in @xcite ) changes from @xmath34 to @xmath35 and from @xmath36 to @xmath37 , respectively , and
the mean is shifted to @xmath38 in both cases .
the shear is then measured independently in all available filters and a weighted average is calculated from those to arrive at our final , psf corrected ellipticity estimates .
we use the @xmath39 overlap between the two f606w pointings to test for systematics in the psf correction .
the differences in the shears derived from the two exposures is consistent given the noise level .
the acs camera has a psf that varies both spatially and temporally .
while the approach described above accounts for the spatial dependence , it does not account for time variations .
we investigate the psf time dependence by measuring stellar shapes in individual exposures ( four exposures were taken within each orbit ) .
we conclude that the psf is _ not _ stable ( i.e. the global pattern of the psf changes ) even within a single orbit .
however our detected signal is much stronger than the typical residuals due to imperfect psf removal caused by the temporal variation ; therefore we use the combined images to measure shapes .
in addition , the psf shape of an individual star changes as a function of its radius ( see e.g. @xcite ) ; we therefore match the radius used to measure the stellar correction factors with that for the galaxies .
using simulated weak lensing data we identified a constant bias between the measured and estimated shear values ; this bias was calculated for each image set separately and was introduced as a multiplicative shear calibration factor when producing final catalogs .
the color cuts applied to each catalog to remove cluster and foreground galaxies were color - color cuts estimated using photometric - redshift templates of galaxies at redshifts @xmath40 .
we used templates of all but extreme starburst galaxies , under the assumption that there are not likely to be many faint , low - redshift starburst dominated galaxies , particularly in a cluster . as such , we expect that unless there is an unknown population of dwarf galaxies in the cluster with colors much different from known populations , we will have removed all but a handful of outliers ( extreme starbursts or dust obscured ) .
we only use objects which were detected in more than a single band ; in particular we reject all the objects that are detected only in the f606w pointing of the main cluster and were too faint for detection in the magellan images .
unfortunately the observations at hand do not allow us to determine the photometric redshifts of the sources used for weak and strong lensing .
we would need at least nir data to reliably determine the redshifts , as most of the sources are located at @xmath41 .
therefore following @xcite , we estimate the redshifts of the background sources by applying the color and magnitude cuts that remove the likely foreground population to the hdf - s photometric redshift catalog of @xcite . for the remaining sources , we average the ratios @xmath42 ( where @xmath43 , @xmath44 , and @xmath45 are the angular diameter distances to the cluster , source and between the cluster and the source respectively ) to generate a very crude estimate of the individual redshifts , but good average redshifts per population ( provided the redshift distribution of the hdf - s is representative of that of our background sources ) . as a result
all of the galaxies in a certain magnitude bin will be given the same redshift .
the main reason of this approach , instead of calculating a mean for the whole catalog is the fact that the images have different depths in different parts of the reconstructed area , and this method allows us to account for the sudden jumps in the mean background galaxy redshift across the image boundaries , which would otherwise introduce spurious features in the reconstructions .
this method is only approximate , therefore we test the reconstruction method using mock catalogs with only limited knowledge of the redshift of the source population in sect .  [
sec : swunited ] .
the weak and strong lensing reconstruction performs well under these conditions .
in addition , any significant dilution of the signal will be compensated by the strong lensing signal ( if more than one strong lensing system is used ) .
the cluster is also at a moderate redshift , therefore changing the sources from redshifts of @xmath46 to @xmath47 would change the lensing strength by @xmath48 . in conclusion
we expect the contributions due to errors in the assumed distribution of @xmath42 to be smaller than that expected from the galaxy ellipticity shot noise .
finally , we fit the non - singular isothermal sphere nis model to the the averaged tangential ellipticities @xmath49 .
both components generate large tangential shear signal .
we center circular bins on the southern cd of the main cluster and the bcg of the subcluster .
the resulting nis fit for the main cluster results in @xmath50 , @xmath51 , and @xmath52 , @xmath53 for the subcluster .
this is an oversimplified model because the cluster is morphologically disrupted and we do not correct for the contamination of the main cluster profile in the sub cluster region ( and vice versa ) .
our strong and weak lensing reconstruction is performed using the method described in .
the main idea behind the method is to describe the cluster mass distribution using a set of generic model parameters , which we choose to be the projected gravitational potential @xmath54 measured on a regular square grid .
the total @xmath55 is defined as @xmath56 where @xmath57 is the contribution from weak lensing , @xmath58 from strong lensing and @xmath59 is the regularization function with regularization parameter @xmath60 .
we closely follow the methodology presented in , with the few differences described below . in finite differencing formula for @xmath7
all the calculations in all three papers do however use the correct factor of @xmath61 . ]
the most important difference is that we extend the formalism to more than a single set of multiple images .
this may appear straightforward because we only need to sum up the contributions to @xmath62 from each system .
given @xmath63 multiply imaged systems having @xmath64 images each , we define @xmath62 as @xmath65 where @xmath66 is the covariance matrix and @xmath67 .
@xmath68 is the average position of the @xmath69-th source ( calculated from the image positions and given the deflection angle @xmath70 ) , and @xmath71 is the so - called `` cosmological weight '' function as defined in @xcite .    in we assumed that the covariance matrix is diagonal and independent of the lens model parameters . as discussed in e.g. @xcite ,
this approach is not optimal and can lead to solutions that are biased toward models that predict high magnification at the image positions . with more sources at differing redshifts , this problem becomes more apparent because the method likely converges toward a solution with steep gradients in the surface mass density close to image positions ( allowing for the `` fine tuning '' of the reconstruction ) . to overcome such unphysical solutions
, we use a better approximation for the error containing the magnification at the image position @xmath72 and obtain the covariance matrix @xmath73 .
such a procedure is approximately correct , provided we are close enough to the true model ( see @xcite ) .
the covariance matrix is still diagonal ; however we argue that assuming it to be such in the source plane is in fact a better approximation than assuming the covariance matrix to be diagonal in the image plane , as sources are on average more circular than their lensed images .
in addition , in practice multiple image constraints are satisfied nearly perfectly and exact values of errors on image positions are of lesser importance .
we adapt the method to accommodate data regions of arbitrary size and shape , but the potential is reconstructed only in grid cells containing weak and/or strong lensing data .
we tile the observed field with a regular grid and determine which cells contain either weak lensing galaxies or multiple images .
this can result in some cells being sparsely populated ( holes in the data , edge effects ) , but the regularization ensures that potential is properly reconstructed . once the data - grid is determined , we ensure that each of the data grid points is surrounded by a @xmath74 grid of points that are included in the reconstruction to calculate the scaled surface mass density @xmath7 , shear @xmath75 , and the deflection angle @xmath76 , using finite differencing , as described in .    to find the @xmath55-function minimum , we search for a solution of the following system of equations @xmath77 .
this is in general a non - linear set of equations , which we solve in an iterative manner .
we again linearize this system as described in , calculating all of the non - linear terms ( for the strong lensing term in particular these are now the ones containing @xmath68 and @xmath72 ) using the information from the previous iteration .
we regularize the solution using a `` moving prior '' , gradually updating our knowledge of the model by increasing its complexity and computing its likelihood relative to the previous ( simpler ) model .
we begin with an initial model ( @xmath78 , @xmath79 , @xmath80 ) computed on a relatively coarse grid , and then gradually increase the number of grid points , comparing the resulting @xmath7 map with that from the previous iteration , linearly interpolated onto the finer grid .
in addition to @xmath7 we now also use the components of the shear @xmath81 and @xmath82 in the regularizing scheme , to avoid the numerical effects described in even more efficiently .
the regularization constant was set after experimentation with simulated data .
when using many multiply imaged systems , the final reconstruction is even less sensitive to these choices than when only a single system is used ( as in ) ; we discuss this finding in more detail in sect .
[ sec : systematics ] .    with
many multiply imaged systems , we must choose a sufficiently fine grid to be able to calculate the difference in the deflection angle between images of different systems ( if images from more than one system lie within the same grid cell , it is difficult to find a good solution ) .
ideally one would like to use adaptive grids for this purpose ( having a relatively coarse grid where only weak lensing information is available and increasing the resolution in the vicinity of the multiple images ) .
unfortunately it is computationally difficult to track such a problem and is therefore beyond the scope of this paper , but will be a subject of future work .
we tested all the improvements described above using a high - resolution n - body simulation of a galaxy cluster by @xcite .
the catalogs are generated as in , although adapted to the higher - resolution data that we use here .
the weak lensing simulated data are obtained by placing @xmath83 galaxies ( giving a density of @xmath84 ) randomly on a @xmath85 field .
the intrinsic ellipticities @xmath86 are drawn from a gaussian distribution , with each component characterized by @xmath87 .
we then draw the redshifts of the background sources , following @xcite , from a gamma distribution @xmath88 with @xmath89 , giving a mean redshift of @xmath90 , and a mode of @xmath91 .
we reject the foreground objects , and use the background object redshifts when predicting the lensed shapes .
however , when using the catalogs for the reconstruction , we assume _ no _ knowledge of the individual redshifts : instead all sources are assumed to have @xmath46 , in agreement with the poor estimates available for the actual data .
the multiple - image systems are generated as described in .
four quadruply - imaged systems are generated at source redshifts of @xmath92 .
the initial model for the regularization is set to @xmath93 .    the resulting reconstruction is given in fig .
[ fig : sims ] .
it is very encouraging to observe how well the overall shape ( the agreement in the ellipticity and position angle of the main component , and the overall profile ) and its mass are reconstructed .
these successes would not be possible with weak lensing data alone , owing to the mass - sheet degeneracy ( i.e. the reconstruction would suffer from the degeneracy of the form @xmath94 , where @xmath95 is an arbitrary constant , see e.g. @xcite ) .
in addition with the acs - quality data , we are starting to resolve the small - scale substructure of the potential .
we now apply the above mentioned reconstruction method to the strong and weak lensing data of . as described in sect .
[ sec : datasl ] , we use 10 distinct multiply imaged systems and @xmath96 `` weakly '' lensed galaxies ( @xmath97 ) .
we start with a @xmath98 grid for a @xmath99 field oriented with respect to the negative ra coordinate ( not all of the grid cells contain data and those that do not are excluded from the reconstruction ) .
we gradually increase the number of grid points ( for reasons described in sect .  [
sec : swunited ] ) in steps of 1 , with the final reconstruction performed on a @xmath100 grid . in our initial model , we use a two - component nie model with parameters as given in sect .
[ sec : datasl ] to set the regularization .
the possible dependence of the reconstruction on initial conditions is discussed in the next section .    the resulting reconstruction is shown in fig .
[ fig : swunited ] , together with x - ray surface brightness contours from the 500 ks chandra observation . the two cluster components are clearly detected .
in addition , with superb strong lensing data ( due to exquisite acs resolution and multi - color information for the subcluster ) we are starting to resolve multiple components within each individual system .    from this reconstruction ,
we measure the enclosed ( projected within a cylinder ) mass for this system .
the main cluster has @xmath101 , while the sub cluster has @xmath102
( @xmath103 at the cluster redshift ) .
the mass errors represent combined errors obtained by bootstrap resampling of the weak lensing catalogs , performing a similar experiment for the strong lensing data , and estimating systematic uncertainties from the mass differences obtained using different initial models ( see discussion in the next section ) .
we plot the integrated , azimuthally - averaged mass profiles @xmath104 for the two components in fig .
[ fig : profile ] .
we fit power - law models ( @xmath105 ) to these profiles , and find logarithmic slopes of @xmath106 for the main cluster and @xmath107 for the subcluster . in terms of the surface mass density @xmath108 , this translates to @xmath109 , and @xmath110 . adopting an isothermal profile ( @xmath111 ) , we find the line of sight velocity dispersions for the main and the sub cluster are @xmath112 and @xmath113 respectively ; in good agreement with the analysis in sect .
[ sec : datawl ] .
however , as already mentioned in that section , given the size of the errors , the asphericity of the system , and the complicated geometry ( the relative contributions of the main and sub cluster are difficult to disentangle ) , the velocity dispersion and slope measurements should only be used as guidelines .
because of this cluster s violent history , we can not expect the assumption of isothermality to hold .
therefore ( and for the reasons stated above ) it is also non - trivial to compare our velocity dispersion measurements with the ones from @xcite .
their estimate for the main cluster agrees with ours , while the velocity dispersion of a subcluster is much lower in their case .
an isolated lens with the low velocity dispersion measured by @xcite for the subcluster would have multiply imaged systems only at radii @xmath114 , a factor of @xmath115 below that observed .
the influence of the main cluster can increase the radial range over which multiple images exist , but not to the extent observed here .
we conclude that the true velocity dispersion of the subcluster is likely to be smaller than that derived from our mass model ( due to the inclusion of mass from the main cluster ) , but larger than that presented by @xcite .
one of the important results of the work of @xcite was the measurement of the offset between dark matter and baryons .
we confirm these results : our strong and weak lensing offsets for the mass peaks ( by measuring the peaks in @xmath7-distribution ) and the corresponding errors are given in table  [ tab : offset ] . to estimate the significance of this offset we also measured the positions of the two peaks ( for the main and the subcluster ) in each of the 1000 reconstructions using the bootstrap resampled weak lensing catalogs ( described below ) . in none of these realizations
are the mass peaks coincident with the corresponding x - ray peaks .
the 1-@xmath11 error bars are a factor of @xmath116 and @xmath6 smaller than the offsets for the main and subcluster respectively .
if these errors are accurate and gaussian distributed , this result translates into 10-@xmath11 and 6-@xmath11 discrepancies between the corresponding gas and the total mass peaks .
these offsets directly demonstrate that dark matter is present and is the dominant mass component of this cluster .
lastly , we compare the projected gas mass ( the details of modeling will be given in @xcite ) vs. the total mass within elliptical regions centered on the two peaks of the x - ray emission and the southern cd of the main and bcg of the subcluster , respectively ( see table  [ tab : massgas ] ) .
the ratio of gas mass to total mass is smaller in the regions centered on the massive galaxies rather than the x - ray peaks , indicating that the gas has been ( to some extent ) stripped , and is no longer coincident with the total mass . in addition
, we detect an extension in the total mass map to the nw of arc a. this detection can be partly attributed to the hot x - ray gas ( see fig .  [
fig : swunited ] ) .
rrr main component & @xmath117 & @xmath118 + main component ( gas ) & 1.06 & 0.84 + sub component & @xmath119 & @xmath120 + sub component ( gas ) & -0.73 & 0.10 [ tab : offset ]
to assess the reliability of the reconstructed map , we generate 1000 bootstrap resampled weak lensing catalogs and perform reconstructions on each of these .
to further test the reliability of the strong lensing data we create 10 different reconstructions each time removing one of the multiply imaged systems we use .
the resulting @xmath7-maps do not change substantially , the main features ( i.e. the ellipticity of the main and the subcluster ) seen in fig .
[ fig : swunited ] remain in all of the reconstructions .
we also perform a reconstruction using the weak lensing data only ( bottom - left panel of fig .  [
fig : weakvsstrong ] ) . from only weak lensing
, we clearly detect both the main and sub clusters , but are unable to resolve individual substructures within each with high significance .
this failure can in part be attributed to the limited field size used here in contrast to @xcite ( our method at present does not allow us to use larger fields for the reconstruction , adaptive grid methods need to be employed then ) .
in addition , the mass estimates from weak lensing reconstruction are unreliable ( due to the the mass - sheet degeneracy ) ; the resulting potential does not reproduce any of the multiply imaged systems observed .
traditionally , to overcome the problem of the mass - sheet degeneracy and measure the mass from the weak lensing mass reconstruction , one assumes a vanishing surface mass density at the edge of the field or adopts parametrized model fitting .
the former is not applicable here because the field used is too small .
the latter will yield reliable mass estimates if the assumed profile is correct , but it is difficult to guess the appropriate model for such a morphologically disturbed cluster . combining strong and weak lensing helps us to break this degeneracy .
however , the lack of counter arc candidates and redshifts of multiply imaged systems limits our ability to provide even stronger constraints on the mass distribution .
the errors on mass estimates quoted throughout the paper are obtained by summing in quadrature the errors obtained from the reconstructions using bootstrapping of the weak lensing catalogs and removing individual strong lens systems .
in addition , we estimate and add systematic uncertainties resulting from differences in the reconstruction when using different initial models ( as explained below ) . in principle
, we should incorporate all of these effects simultaneously in the reconstructions , however this is too computational intensive and beyond the scope of this analysis .      methods involving strong lensing mass reconstruction with many multiple imaged systems are subject to degeneracies in parameter space .
it is difficult to search the parameter space because the @xmath55 function has many local extrema . in our case , there are pronounced , narrow secondary minima where individual systems are reconstructed well by the model ( and the others less so ) and high `` walls '' in @xmath55 space where images happen to form close to the corresponding critical curves . because weak lensing data are noisy , additional , though less pronounced , minima in @xmath55 are created . given that our reconstruction method does not guarantee that we find a global minimum , in principle we would need to check the whole parameter space .
this search is unfeasible with any available algorithms , given that the number of parameters ( i.e. , pixels in the mass map ) is @xmath121 .
instead , we opt to investigate the use of many initial models .
we try both extremes ; i.e. models that are clearly in disagreement with the data ( zero initial conditions , single mass peak centered on the position of the gas ) and models that use alternative methods to reconstruct the data ( two - nie model , two - nie plus cluster members model from strong lensing modeling , weak lensing mass model from @xcite ) .
although the details of the reconstruction depend on the initial conditions ( as can be seen from comparing all but bottom - right @xmath7-maps in fig .
[ fig : weakvsstrong ] ) , the main features ( se - nw elongation of the main cluster , two mass peaks close to the main cluster cds , an elongation in the mass map toward the gas peak , e - w elongation of the subcluster ) are independent of the particular initial model we use .
the resulting @xmath55 contribution of the weak lensing data @xmath122 remains the same for all sets of initial conditions ( @xmath123 for @xmath124 data points ) .
the strong lensing contribution is lowest when the two component nie model is used ( @xmath125 ) and highest when the zero model initial conditions ( @xmath126 ) are used .
although the absolute values of the @xmath55 have little meaning ( i.e. one can change the regularization parameter to obtain a lower value for @xmath122 because the weak lensing data is noisy ) , the relative values tell us which resulting reconstruction is indeed a better fit to the data ( as the same regularization parameters are used in all cases and it is the fit to the high signal to noise strong lensing data that improves ) .
we estimate the systematic uncertainty on our mass estimates by calculating the maximum difference in the resulting mass estimates when using different initial models .
we exclude the zero model in this calculation as it does not predict the strong lensing data reliably .
the resulting @xmath7-map ( see bottom - left panel in fig .  [
fig : weakvsstrong ] ) for this model is suppressed where only weak lensing data are available ( due to the mass - sheet degeneracy ) and produce steep gradients close to arc a and a @xmath127 lower overall mass .
however , other features in the reconstruction ( positions and shapes of the mass peaks for the main and the sub - cluster ) remain unchanged . in all cases ,
the resulting mass peaks `` move '' away from the initial guess ; showing that the regularization is sufficiently flexible to reproduce the features in the data .
therefore , we are confident that our results are not biased as a result of the particular choices of regularization and are independent of the initial conditions to the extent presented in fig .  [
fig : weakvsstrong ] .      finally , the mass measurements presented above depend on the reliability of the redshifts and identification used for strong lensing reconstruction .
these ( except for system a ) are not secure and depend upon the model used in sect .
[ sec : datasl ] .
in particular , the redshifts of the systems around the subcomponents ( and consequently the measurements of its mass ) can change because the subcomponent mass reconstruction does not depend strongly on the modeling around the main component ( where the arc with the known redshift is located ) .    the parametrized model used in sect .
[ sec : datasl ] to estimate the redshifts of the strong lensing systems likely deviates from the truth model .
we tested the dependence of the assumed redshifts by performing 100 reconstructions in which we modified all of the strong lensing system redshifts except for system a. for each reconstruction these uncertain redshifts were changed by assuming a gaussian scatter with @xmath128 and requiring that all sources lie behind the cluster .
the resulting mass and position errors were smaller ( by a factor of @xmath1292 ) than the ones quoted throughout the paper , and the mass and position measurements were unbiased ( as expected from such an analysis ) .
the mass measurements would change , however , if we systematically increase ( decrease ) the redshifts of all the systems .
for example , if we decrease the redshifts of the subcluster systems such that the lowest redshift system of the subcluster ( g ) is at @xmath130 ( and scale the others accordingly ) , then the mass increases by a factor of @xmath131 . if we instead increase all the redshifts , pushing the highest redshift system ( i ) to @xmath132 ) , then the mass of the subcluster will decrease by @xmath133 .
note however that these are just estimates , since the redshift estimates of the weak lensing sources also plays a role , and the changes will therefore be even smaller .    as a further test
, we also perform 10 reconstructions where we remove one strong lens system at a time from the analysis . in fig .
[ fig : slsyst ] we show the two that differ the most from the reconstruction that uses all the strong lensing systems . when we remove system i from the reconstruction ( right panel in fig .
[ fig : slsyst ] ) the double peaked structure of the subcluster disappears , while the east - west elongation remains . for the main cluster
we see the largest deviations from the original map when we remove system d from the reconstruction ( left panel in fig .  [
fig : slsyst ] ) .
however , these results do not mean that these structures are not real , rather they illustrate some possible systematics when dealing with sparse strong lensing data .
we are confident that system d is indeed multiply imaged , and images i2-i3 , at least , are also probably multiply imaged .
these tests demonstrate that even though the strong lensing system identification and their redshift determinations are important to the details of the reconstruction , the errors arising from misidentifications / redshift misestimation are within those quoted in the previous section .      the predicted redshifts and strong lensing system identification can be tested using spectroscopy .
we are planning just such a survey and will therefore be able to test the fidelity of our model in the future .
in addition , using spitzer irac data we identify a candidate high redshift object . while bright in all irac bands , the faintness of this galaxy in optical and near infrared bands indicates that it could lie at @xmath134 .
it has two components and our predicted critical curve for this redshift ( indistinguishable from the @xmath135 critical curve ) is exactly where it should be if indeed these sources are merging images of the same source ; they are denoted with crosses and plotted together with the critical curve in fig .
[ fig : ccurve ] .
furthermore , we identify a candidate counter image with similar colors in the irac data just outside the critical curve .
all three need to be further analyzed and confirmed ( using spectroscopy and/or nir data ) before we include them in the modeling .
we also identify possible counter - images for systems b , j , h , and i as predicted from our model , they are indicated in fig .
[ fig : predict ] .
if these and high - redshift candidate are confirmed , the predictive power of our model would further indicate that the strong lensing data we use is reliable enough for this work .
, the results are shown for the main ( solid line ) and for the subcluster ( dashed line).,scaledwidth=50.0% ]    rrrrr main ( gas peak ) & & 2.0 @xmath136 0.2 & 10.8@xmath136 0.6 & 0.19@xmath136 0.03 + main ( s cd ) & & 1.4 @xmath136 0.1 & 16.3 @xmath136 0.9 & 0.09 @xmath136 0.01 + sub ( gas peak ) & & 0.42 @xmath136 0.04 & 2.1@xmath136 0.2 & 0.20 @xmath136 0.06 + sub ( bcg ) & & 0.18 @xmath136 0.02 & 4.3 @xmath136 0.2 & 0.04 @xmath136 0.01 + acs field & 1300 ( 4.9 ) @xmath137 830 ( 3.2 ) & 13 @xmath136 1 & 91 @xmath136 10 & 0.14 @xmath136 0.03 [ tab : massgas ]
massive and interacting clusters , while quite rare , are remarkably well - suited to addressing outstanding issues in both galaxy evolution and fundamental physics . however , in order to study the mass distribution , methods relying on hydrostatic ( x - rays ) or dynamical equilibrium are ill - suited for such systems .    we have applied a mass reconstruction method based on strong and weak gravitational lensing to the cluster .
we use deep , high resolution optical data to identify objects belonging to the same multiple - image systems .
the same data are used to obtain weak lensing catalogs allowing us to obtain a strong+weak lensing mass map of the cluster core .
our main conclusions are the following :    1 .   using the combined strong and weak lensing mass reconstruction we derive a high - resolution , absolutely calibrated mass map ; we get projected , enclosed mass @xmath138 around the main and @xmath139 around the subcluster .
we detect the main cluster peak and a distinct mass concentration at the subcluster position , both clearly offset from the location of the x - ray gas in the system ( at 10-@xmath11 and 6-@xmath11 significance for the main and the subcluster respectively - see fig .  [
fig : swunited ] and tab .
[ tab : offset ] ) .
the majority of the mass is spatially coincident with the galaxies , which implies that the cluster mass must be dominated by a relatively collisionless form of dark matter .
the high resolution data allow us to significantly detect the shapes of both the main mass component and the subcluster with no prior assumptions on their positions or profiles .
we show that the cluster is a very efficient strong lens ; the area enclosed within the critical curve for highest redshift sources is comparable to that of a1689 @xcite .
consequently , we plan to conduct a pencil beam survey for high - redshift ( @xmath140 ) objects . at present
this is not possible , because we lack both reliable redshift information for the multiple image systems , and multi - color high - resolution imaging data in the main cluster region .
the latter is needed to identify counter images for the strong lens systems , and so improve the reliability of the mass map to the south - east of the main cluster .
once these data are obtained ( hst cycle 15 ) , we will improve our methodology to include an adaptive grid approach , which will greatly improve the constraints on the location of critical curves for these high - redshift sources .
the prospects for such a search are however very promising , as we already have two strong candidate high redshift objects .
one is the red arc b , which is constrained to lie at @xmath141 based upon the redshift of arc a and detection in f814w .
the other ( see fig .  [ fig : ccurve ] ) is spitzer irac - selected with a likely redshift of @xmath134 .
using this cluster as a gravitational telescope will therefore allow us to study galaxy properties near the epoch of reionization .
the high resolution mass map of this cluster , in addition to enabling the high - z galaxy search , will be valuable for several other ongoing science programs .
most significantly , we are currently running simulations to reconstruct the dynamical history of the cluster merger , and provide an independent constraint on the dark matter self - interaction cross - section .
the relative masses of the main and the subcluster , and the central mass distribution , two critical parameters for the simulations , were the limiting source of uncertainty prior to this work .
we would like to thank marco lombardi , roger blandford , peter schneider , and steve allen for many useful discussions that helped to improve the paper , and mischa schirmer for producing the color image used for fig .
[ fig : predict ] .
we are grateful to the anonymous referee for valuable comments and suggestions that helped to improve this manuscript , and for pointing out a typo in our previous paper of the series .
further we would like to thank volker springel for providing us with the simulations used in the paper .
support for this work was provided by nasa through grant number hst - go-10200.05-a from the space telescope science institute , which is operated by aura , inc .
, under nasa contract nas 5 - 26555 .
mb acknowledges support from the nsf grant ast-0206286 .
this project was partially supported by the department of energy contract de - ac3 - 76sf00515 to slac .
dc and dz acknowledge support from nasa through grant number ltsa04 - 0000 - 0041 .
mm acknowledges support from chandra grant go4 - 5152x .                          ,
d. , schneider , p. , aragn - salamanca , a. , bremer , m. , de lucia , g. , halliday , c. , jablonka , p. , milvang - jensen , b. , pell , r. , poggianti , b. , rudnick , g. , saglia , r. , simard , l. , white , s. , & zaritsky , d. 2006 , , 451 , 395        , p. , romer , a.  k. , peterson , j.  b. , chase , w. , runyan , m. , holzapfel , w. , kuo , c.  l. , et  al .
2004 , in aip conf .
proc . 703 : plasmas in the laboratory and in the universe : new insights and new challenges , 361366      , c. , brown , m.  l. , barden , m. , caldwell , j.  a.  r. , jahnke , k. , peng , c.  y. , rix , h .- w . , taylor , a. , beckwith , s.  v.  w. , bell , e.  f. , borch , a. , huler , b. , jogee , s. , mcintosh , d.  h. , meisenheimer , k. , snchez , s.  f. , somerville , r. , wisotzki , l. , & wolf , c. 2005 , , 361 , 160    , c. , van waerbeke , l. , bacon , d. , berge , j. , bernstein , g. , bertin , e. , bridle , s. , brown , m.  l. , clowe , d. , dahle , h. , erben , t. , gray , m. , hetterscheidt , m. , hoekstra , h. , hudelot , p. , jarvis , m. , kuijken , k. , margoniner , v. , massey , r. , mellier , y. , nakajima , r. , refregier , a. , rhodes , j. , schrabback , t. , & wittman , d. 2006 , , 368 , 1323


for the past two decades there has been a revival of interest in issues related to charge migration in dna based oligomers .
on one side , charge transfer ( ct ) in dna is supposed to play a key role in self - repair processes of dna damage in natural conditions , i.e. via oxidative stress  @xcite . on the other side , dna may have a huge potential application field in the development of complex nano - scale electronic devices with self - assembling properties , since self - assembly , recognition and fully - automatic synthesis of oligonucleotides  @xcite can allow for building almost any imaginable two and three dimensional dna motifs ( e.g. dna - origami )  @xcite .
however , besides considerable progress in elucidating the possible charge transfer mechanisms efficient in dna , there is still disagreement about the conduction properties of dsdna , for several experiments obtained quite controversial results ranging from insulating  @xcite , over semi - conducting  @xcite to even metallic - like behavior  @xcite .
this broad variability is related not only to the poor control of the experimental conditions in earlier transport experiments , but is also due to the intrinsic structural and electronic complexity of dna oligomers ( e.g. specific base pair sequence ) , which makes the formulation of a general theoretical framework very challenging .
nevertheless , several experiments appeared recently  @xcite which have shown electrical currents in the na range despite the differences in the studied base sequences and in the experimental setups .
one central issue which has emerged by the theoretical treatment of charge transport is the necessity to consider the conformational fluctuations of the molecular frame not as a weak perturbation but rather as a crucial factor promoting or hindering charge motion  @xcite .
it has been meanwhile shown that it is possible to synthesize dna derivatives , which do not necessarily have the double - strand structure of natural dna .
thus , c and g - rich dna strands are able to form four - stranded quadruplex structures , like the i - motif structure composed of two parallel hemiprotonated duplexes intercalated into each other into a head - to - tail orientation , and the g4-quadruplex ( also known as g4-dna ) which is formed by either one , two or four g - rich dna strands in a parallel or anti - parallel orientation  @xcite .
the latter is supposed to play an important role in some biological processes such as in telomeric dna regions , where it inhibits telomerase and hiv integrase  @xcite . additionally , g4-dna is known to interact with various cell proteins that cause diseases like bloom s and werner s syndromes  @xcite .
moreover , g4-dna is cytotoxic towards tumor cells and therefore , might be a key for the design of anti - cancer drugs  @xcite . eventually , g - quadruplexes are found to be thermally more stable than dsdna  @xcite .
gilbert and sen resolved the x - ray structure of the first g - quadruplex and proposed the formation of these unique structures in presence of monovalent alkali ions  @xcite , see also fig .
[ fig : meth - tetrad](a ) .
consequently , g4-dna can be regarded as stacks composed of individual planar g - tetrads ( or g - quartets ) as shown in fig .
[ fig : meth - tetrad](b ) , each formed by 4 guanines bound together by 8 hydrogen bonds via hoogsteen pairing  @xcite . for detailed reviews on the structural diversity and special mechanical properties of g4-dna based molecules see ref .
@xcite .
recently , porath _ et al . _
@xcite have shown the larger polarizability of g4-dna compared to dsdna , when attached to gold surfaces .
this might be an indicator of improved ct properties and the potential advantage of using g4-dna rather than dsdna in molecular scale devices .
eventually , the higher conductance in g4-dna could be attributed to the increased structural stability and higher number of overlapping @xmath0-orbitals .
furthermore , experimental work by kotlyar _
et al . _
@xcite revealed that up to 300 nm long g4-wires could be synthesized revealing a considerable stability even in the absence of metal ions .
also , theoretical studies have been performed on the stability and rigidity of g4-dna .
molecular dynamics ( md ) simulations by pakov _
et al . _
@xcite could confirm the experimental results .
for one thing , g4-dna is more rigid than dsdna and for another monovalent ions within the quadruplex cavity are necessary for the stability of short g4-dna .
on the other hand , md simulations  @xcite pointed out the stability of longer g4 molecules ( with 24 g - tetrads ) in absence of metal ions which is consistent with the synthetic procedure of ref  @xcite .
electronic structure calculations by di felice and co - workers  @xcite showed a higher degree of delocalization of the electronic states in g4 than in dsdna .
moreover , the presence of metal ions may contribute additional states supporting ct , though this is a still unresolved issue .
guo _ et al . _
@xcite have recently shown using the landauer theory ( coherent transport ) that g4-dna exhibits much larger delocalization lengths at the band center compared to double - stranded poly(g ) dna .
surprisingly , it was also found out that the delocalization length can be even enhanced via environment - induced disorder through the backbones .
though these studies assume a static atomic structure , they nevertheless suggest that disorder plays an important , non - trivial role in mediating ct in dna .
similarly , it is known from other biomolecules such as proteins that electron transport is dominated by non - equilibrium fluctuations which has lately been analyzed by balabin _
et al . _
as demonstrated in previous studies , idealized static structures are not representative when considering ct properties , since dynamical as well as environmental effects were shown to be too important  @xcite .
it was indicated that dynamic disorder has dramatic effects , since it suppresses ct in homogeneous sequences on the one hand , but can enhance ct in heterogeneous sequences , on the other .
for instance , in our previous work  @xcite , the conductance of the dickerson dodecamer ( sequence : 5-cgcgaattcgcg-3 ) was found to be almost one order of magnitude larger in solution ( qm / mm ) than in vacuo .
interestingly , similar findings have recently been shown by scheer and co - workers , for they obtained for a heterogeneous sequence ( i.e. a and g bases are present ) with 31 base pairs  @xcite a two order of magnitude larger conductance in solution compared to the in vacuo measurements  @xcite .
furthermore , it could be shown that only a minor part of the conformations is ct - active  @xcite . therefore , neglecting these significant factors or assuming a purely random disorder distribution can lead to a considerable loss of a vital part of ct - relevant structural and electronic information .
relying on our previously developed methodologies to deal with ct in different dna oligomers  @xcite , we present here a detailed investigation of charge transport in g4-dna .
the paper is organized as follows , in sec .
[ sec : meth ] details of the md simulations and of the electronic structure approach are described .
subsequently , the md and the electronic structure data are analyzed in sec .
[ subsec : res - md ] and [ subsec : res - es ] , respectively . finally in sec .
[ subsec : res - landauer ] , ct results for g4 molecules are compared with those of dsdna in solution and in vacuo .
more importantly , we point out significant factors which are responsible for the enhanced conductance in g4 molecular wires .
the molecules used in this work are based on the x - ray crystal structure of a tetrameric parallel - stranded quadruplex ( 244d ) formed by the hexanucleotide sequence d(tg@xmath1 t ) in the presence of sodium ions .
the structure has been resolved at 1.2  resolution by laughlan _
et al . _
@xcite . as shown in fig .
[ fig : meth - tetrad](a ) , it contains 2 pairs of quadruplexes ( tg@xmath1t)@xmath1 , for which each pair is stacked coaxially with opposite polarity at the 5 ends .
moreover , 9 sodium ions are located inside each stack of quadruplexes illustrating well ordered g4-dna constructs .
however , the terminal thymine residues , shown in gray , were not completely resolved because of high thermal disorder  @xcite . as indicated in fig .
[ fig : meth - tetrad ] only one of these 4 quadruplexes is taken as our basis structure , which corresponds to the parallel strands a , b , c and d from the pdb file 244d .
this quadruplex , from now on denoted as ( tg@xmath1t)@xmath1 , will be used for simulations and calculations .
furthermore , two longer g4 quadruplexes with 12 and 30 tetrads , denoted as ( g@xmath2)@xmath1 and ( g@xmath3)@xmath1 , are generated by omitting the terminal thymine residues and adding subsequently g4 tetrads with a distance of 3.4  and twisted by 30 . for comparison ,
corresponding double - stranded b - dnas with base sequence 5-tggggt-3 , poly(g ) denoted as g@xmath2 and g@xmath3 , respectively , and a sequence containing 31 base pairs  @xcite ( scheer ) are built with the make - na server  @xcite .
backbones are indicated as light blue ribbons , central coordinated sodium ions as dark blue spheres , terminal thymine and guanine residues are shown in gray and orange , respectively .
the highlighted quadruplex is used as starting structure and also to generate longer g4 molecules ( g@xmath2)@xmath1 and ( g@xmath3)@xmath1 .
b ) lewis structure of a single g4-tetrad with c@xmath4 symmetry containing a monovalent metal ion in its center .
guanines are bound together by 8 hydrogen bonds via hoogsteen pairing  @xcite .
the metal ion is coordinated either coplanar by 4 or cubic by 8 o6 oxygen atoms of the respective guanines which depends on the ionic radius of m@xmath5 . ]    depending on the helix length the g4 and dsdna molecules are solvated in a rectangular box with 4000 to 11000 water molecules using the tip3p model  @xcite and periodic boundary conditions ( pbc ) are applied . in order to neutralize the total charge of the system , due to the negatively charged backbones , an appropriate number of sodium counterions ( na+ ) are added . concerning the quadruplexes
, we carried out 4 simulations for ( tg@xmath1t)@xmath1 and ( g@xmath2)@xmath6 , respectively , firstly without any central monovalent alkali ions within the quadruplex , and then in presence of either lithium ( li+ ) , sodium ( na+ ) , or potassium ( k+ ) ions within the respective g4 molecules .
therefore , either 3 ( ( tg@xmath1t)@xmath1 ) or 11 ( ( g@xmath2)@xmath6 ) of these ions are subsequently placed in the center between two g4 tetrads coordinated by o6 oxygen atoms ( see also fig .
[ fig : meth - tetrad](b ) ) . nevertheless , the longer ( g@xmath3)@xmath1 quadruplex is simulated only twice in absence and presence of central sodium ions .
all simulations are carried out with the gromacs software package  @xcite using the amber parm99 forcefield  @xcite including the parmbsc0 extension  @xcite . after a standard heating - minimization protocol followed by a 1 ns equilibration phase , which is discarded afterwards , we performed 30 ns ( 50 ns for ( g@xmath3)@xmath1 ) md simulations with a time step of 2 fs .
snapshots of the molecular structures were saved every 1 ps , for which the charge transfer parameters were calculated with the scc - dftb - fo approach as described in sec .
[ subsec : meth - hop ] .      in this section
we will shortly describe how the electronic structure of g4-dna molecules is mapped to a coarse - grained transfer hamiltonian using the fragment orbital ( fo ) approach .
the method has already successfully been applied to dsdna molecules  @xcite .
a more detailed description of the methodology can be found in ref .
@xcite . to begin with
, the electronic structure of the dna system is written in an effective tight - binding basis ( fragment basis ) as : @xmath7 the onsite energies @xmath8 and the nearest - neighbor hopping integrals @xmath9 characterize , respectively , effective ionization energies and electronic couplings of the molecular fragments .
the evaluation of these parameters can be done very efficiently using the scc - dftb method  @xcite combined with a fragment orbital ( fo ) approach  @xcite : @xmath10 and @xmath11 the indexes @xmath12 and @xmath13 correspond to the fragments of the quadruplexes which , in this model , are constituted by single guanine bases .
accordingly , the molecular orbitals @xmath14 and @xmath15 are the respective highest - occupied molecular orbitals which are obtained by performing scc - dftb calculations for these isolated fragments , i.e. the individual guanine bases . using a linear combinations of atomic orbitals ansatz @xmath16
the coupling and overlap integrals can be efficiently evaluated as : @xmath17 and @xmath18 @xmath19 and @xmath20 are the hamilton and overlap matrices in the atomic basis set as evaluated with the scc - dftb method .
note , since @xmath21 is built from non - orthogonal orbitals @xmath14 and @xmath15 , which is for various problems not suitable , we apply the loewdin transformation  @xcite : @xmath22 the effect of the environment ,  the electrostatic field of the dna backbone , the water molecules and the counter - ions ( including the central alkali ions within the quadruplex ) , is taken into account through the following qm / mm hamiltonian : @xmath23 @xmath24 are the mulliken charges in the qm region and @xmath25 are the charges in the mm region ,  the dna backbone , counterions and water molecules . the coupling to the environment
is therefore explicitly described via the interactions with the @xmath25 charges . in the following
, we will denote the calculation set - up based on the complete expression in eq .
[ eq : qmmm ] as qm / mm ; neglecting the last term will be denoted as `` vacuo '' .
the electronic parameters @xmath8 and @xmath9 are evaluated for every snapshot of the simulations .
subsequently , the transfer hamiltonian is used to calculate the ct in the quadruplexes as described in sec .
[ subsec : meth - tc ] .
once the transfer hamiltonian has been constructed , the transport properties will be calculated using landauer theory for each snapshot along the md trajectories .
for this , we consider a two - terminal set up where the g4 oligomer is contacted to left ( l ) and right ( r ) metallic electrodes through the four terminal bases .
a central quantity to be computed is the molecular ( g4 ) green function , which can be obtained through a dyson equation : @xmath26 here @xmath27 and @xmath28 are self - energy matrices characterizing the coupling to the electrodes . to simplify the calculations
we do not take the full energy dependence of the self - energies into account , but rather use the so called wide band limit , where @xmath27 and @xmath28 are replaced by energy - independent parameters  @xcite : @xmath29 we set @xmath30 and @xmath31 to 1 mev . within the landauer approach , the transmission function @xmath32 for a given set of electronic parameters
is then obtained as : @xmath33\end{aligned}\ ] ] where @xmath34 and @xmath34 are the broadening matrices calculated as the anti - hermitian part of @xmath35 : @xmath36\end{aligned}\ ] ] using the former expressions , the conformational ( time ) dependent electrical current will be simply calculated by : @xmath37 the i - u characteristics presented in this work should be interpreted only qualitatively , for the fermi energy ( @xmath38 ) is artificially placed as average of the onsite energies for all guanine sites and for each snapshot , respectively .
the reader should note that the current could exhibit quite different shapes depending on where @xmath38 is located .
however , only the current - voltage gap is affected by @xmath38 , whereas the maximum current obtained at high voltages is not altered .    the transmission function and the current
are evaluated for every snapshot of the md simulation in order to obtain statistical average quantities .
however , we are aware of the limits of the coherent transport model which is valid only in the adiabatic regime .
hence , we assume that the timescales of the ct process is shorter than the fastest dominant structural fluctuations .
this issue has also been addressed in ref .  @xcite .
first of all , the ( tg@xmath1t)@xmath1 quadruplex from the crystal structure 244d has been simulated in order to validate further simulations with 12 and 30 tetrads . corresponding rmsd time series in absence and presence of central ions as well as molecular snapshots after 30 ns simulation time
can be found in sec .
si.a . in ref .
based on the crystal structure 244d we built up an artificial all - parallel stranded guanine quadruplex ( g@xmath2)@xmath6 without terminal thymine residues and 12 stacked guanine tetrads ( for details see  [ subsec : meth - md ] ) .
as for ( tg@xmath1t)@xmath1 , we carried out 4 simulations , one in absence and three in presence of 11 lithium , sodium and potassium ions initially placed in the center of two adjacent g4 tetrads . the rmsd time series and molecular snapshots after 30 ns md simulation for ( g@xmath2)@xmath6 and poly(g ) ( g@xmath2 ) are given in fig .
[ fig : rmsd2 ] and fig .
[ fig : strucs-12g4 ] , respectively .    )
@xmath1 in absence and presence of centrals ions li@xmath5 , na@xmath5 and k@xmath5 , b ) comparison with double - stranded g@xmath2 . ]    the rmsd values show that without central ions the quadruplex can not be considered as a g4 structure in equilibrium .
in fact , the rmsd increases for the whole simulation time from 1.5 to 3  indicating considerable structural changes .
in contrast , the rmsd for structures containing central ions reveal rather consistent values around 1 to 1.4  with only minor fluctuations suggesting rigid g4 structures with only slight structural deviations .
compared to the ( tg@xmath1t)@xmath1 simulations , the type of ions within the quadruplex has a lesser influence on the structural deviations , since the rmsd values are rather similar . in order to validate these small structural deviations and therefore the rigidity of g4-dna we carried out a rmsd calculation for double - stranded g@xmath2 ( fig .
[ fig : rmsd2](b ) ) .
generally , the rmsd values and their fluctuations for the dsdna molecule are substantially larger than those for the quadruplex structures .
this clearly supports the notion that guanine quadruplexes are much more rigid than corresponding dsdna and therefore g4-dna reveals less dynamical disorder which might promote ct .    )
@xmath1 after 30 ns md simulation in absence and presence of centrals ions li@xmath5 ( green spheres ) , na@xmath5 ( blue spheres ) and k@xmath5 ( purple spheres ) . ]
the snapshots in fig .
[ fig : strucs-12g4 ] exhibit that ( g@xmath2)@xmath6 in absence of central ions is not completely destabilized which has not been the case for ( tg@xmath1t)@xmath1 ( see sec .
si.a . in ref .
the inner tetrads roughly maintain a regular g4 structure , whereas the outer ones show a considerable disordered behavior with no planarity at all , even the g4 hydrogen bridge pattern has been entirely destroyed .
recently , cavallari _
et al . _ have shown in their md simulations that long g4 quadruplexes can be stable without central ions  @xcite .
this finding was supported by experimental results in which longer g4 quadruplexes could be synthesized and characterized in absence of monovalent ions  @xcite .    therefore , an additional 50 ns simulation of a quadruplex with 30 tetrads ( g@xmath3)@xmath1 in absence of central ions is carried out .
the simulations clearly exhibit that the major part of the quadruplex forms a quite regular structure with adequate parallel stacking .
therefore , it seems to confirm previous results that large g4 wires in absence of central ions become more stable , reasonably due to increased stacking interaction .
however , like for the dodecamer , towards the ends the quadruplex is slightly deformed , especially at the 3 terminus .
moreover , rmsd values reveal that the quadruplex is not completely structural balanced .
the major part of structural fluctuation originates from the 5 and 3 ends which is shown by four groups of large peaks in root mean square fluctuations per atom ( rmsf ) . the corresponding data is shown in fig .
s2 in ref .
@xcite .
the quadruplexes which contain central monovalent ions in fig .
[ fig : strucs-12g4 ] exhibit highly regular g4 structures indicating again the strong rigidity and stability of these molecules .
generally , the ions appear to prefer different coordinations within the quadruplex cavity , that is lithium favors the planar position within the tetrads while sodium and potassium are most likely to be found in the center of two g4 planes .
however , one sodium ion in the middle of the quadruplex is located in - plane and two others are rather closer to one of the corresponding tetrads than sandwiched between them .
also , lithium ions are sometimes found slightly apart from the in - plane position .
this might be an indicator for the longitudinal mobility of central sodium and lithium ions within guanine quadruplexes .
in contrast , potassium ions are rather fixed in their center position between two tetrads , most likely because of their larger ionic radii compared with lithium and sodium .
the resulting larger rigidity of g4 compared to dsdna plus the preferred locations of central ions within the quadruplex are in perfect agreement with findings obtained by pakov _
et al . _ as well as previous work by cavallari _ et al . _
@xcite .
notwithstanding , the ions are equally capable of moving out of the quadruplex through the minor groove and the 5 and 3 ends . in our simulations 3 sodium , 3 lithium and 2 potassium ions went out of the g4 dodecamer , respectively .
sometimes , these ions remain in the minor groove for several ps .
likewise , solvent molecules are able to infiltrate the quadruplex via the minor groove and the 5 and 3 ends provided that these locations are not already occupied by initially placed ions . usually , the infiltration of solvent molecules into the helix will not significantly destabilize the structure as long as there is at least a minimum amount of central monovalent cations left in the quadruplex . otherwise , the interior solvation might destroy the regular 4-stranded structure as can be seen for the structures without central ions in fig .
[ fig : strucs-12g4 ] .
these observations are consistent with results of previous works  @xcite .        before calculating ct characteristics
, we have to focus on the electronic parameters , especially the molecular orbitals involved in the hole transfer process . similarly to our previous work on ct parameters in dsdna  @xcite we carried out benchmark calculation of g4 tetrads with dftb and compared to dft ( pbe and b3lyp ) as well as hf calculation .
since we intend to deal with hole transfer we are most notably interested in the highest molecular orbitals which ought to have @xmath0-symmetry in order to have sufficient mo overlap along the quadruplex . tab .
[ tab : mo - ideal ] shows energies of the 5 highest occupied mos for a idealized g4 structure with c@xmath4-symmetry obtained after geometry optimization with b3lyp/6 - 31g(d , p ) .
corresponding visualization of mos can be found in tab .
s3 in ref .
@xcite .
.energies of highest occupied molecular orbitals for a idealized g4 tetrad with c@xmath4-symmetry , comparison between dftb , dft and hf , for the latter the 6 - 31g(d , p ) basis set is used , all values in ev . [ cols="<,^,^,^,^,^",options="header " , ]     as has been discussed above , the central ions stabilize the g4 structure significantly , the structural fluctuations of g4 with ions are much lower compared to g4 without ions or dsdna .
the effect of ions on the electronic parameters is twofold : first of all , the onsite energies are shifted down by about 0.5 ev but , surprisingly , the enhanced structural stability does not lead to smaller fluctuations of the onsite parameters .
they are in the order of 0.4 ev , similar to the situations in dsdna , as discussed recently  @xcite .
these large parameter variations are introduced by solvent fluctuations ( and not by fluctuations of the dna / g4 structure itself ) , therefore , they are not affected by the higher structural rigidity .
the t1 values , on the other hand , are even lower in the g4 structures with ions .
this would indicate , that these structures conduct even less , when compared to four strands of dsdna .
however , the t2 and t3 values in g4 are still large , indicating that interstrand transfer can occur quite frequently .
this opens a multitude of pathways for charge transport in g4 ( compared to dsdna ) , which will be the key to understand g4 conductivity , as discussed in more detail below .    as a first result
, we do not see any indication that the higher structural stability of g4 with central ions will lead to a higher conductivity due to a reduced dynamical onsite disorder or due to increased electronic couplings , i.e. because of somehow better stacking interactions owing to the more regular structure .
therefore , the higher conductivity must have different reasons .
the more stable structure leads to smaller couplings , in contrast to prior expectations @xcite .
a more detailed analysis of onsite energies and electronic couplings is given in tab .
s6-s11 in ref .
@xcite .      the parameters @xmath39 and t@xmath40 analyzed in the previous section are now used to evaluate the charge transport properties using the approach described in sec .
[ subsec : meth - tc ] .
first , we focus on static structures to get insight in particular ct features for g4 quadruplexes and also to have a reference for the following calculations which are applied to the whole md trajectory . later on , we are interested in ensemble averages , for single snapshots or static structure can not elucidate the ct process in dna .
subsequently , we will compare the ct properties of g4-dna with those of double - stranded dna and distinguish the differences in ct efficiency by means of conformational analysis .
additionally , we determine the effects of the mm environment and in particular of the central ions on ct in g4 quadruplexes .      primarily , we compare t(e ) and i(u ) for three static structures : i ) idealized g4-dimer as used in tab .
[ tab : hop - g4 ] , ii ) the two central tetrads of the crystal structure 244d and iii ) double - stranded idealized b - dna ( g@xmath41 ) , i.e. two stacked g - c base pairs ( 36 twist and 3.4  rise ) .
as can be seen from fig .
[ fig : static - strucs ] there are only slight structural differences between the idealized g4-dimer and the stacked tetrads of the crystal structure 244d . fig .
[ fig : g4-trans - ideal ] exhibits the corresponding transmission functions showing eight resonances due to the eight g sites of the transport model .
however , since the idealized g4-dimer has a very regular structure these eight peaks lie in a very narrow energy range of about 0.15 ev and have similar high transmission values of at least 0.9 . on the other hand ,
the eight peaks of the two central tetrads of 244d lie in a broader energy range of about 0.26 ev and have lower transmission values ranging from 0.4 to 0.8 .
this indicates , how small structural deviations lead to a broadening and a decrease of the transmission function . nevertheless , we know from the values in tab .  [ tab : hop - g4 ] that the t1 and t2 couplings are much larger for the crystal structure than for the idealized structure , thus the cause for decreased conductivity has to be the onsite energy disorder .    ) , i.e. two stacked g - c base pairs ( 36 twist and 3.4  rise ) .
]    , comparison of transmission function @xmath32 and current @xmath42 for idealized g4-dimer , the two central tetrads of crystal structure 244d and poly(g ) ( g@xmath41 ) .
note , only the stacked g bases are considered , whereas the cytosine bases and the backbone composed of sugar and phosphate groups are neglected . ]
[ fig : g4-trans - ideal ] also shows the corresponding current - voltage characteristics .
the current for the g4 structures is significantly larger than for g@xmath41 .
this is what we would expect , because we compare two idealized molecular units , one with two relevant transmissive states and the other with eight within the relevant energy window .
note that in all current calculations shown in this work , the fermi energy is artificially placed as the mean value of onsite energies for each snapshot respectively .
generally , the conductance for the idealized double - stranded poly(g ) dna would become twice as large once the cytosine states are included , though it is well established that the c states have only a minor impact in hole transfer due to their higher oxidation potential than g.      previous studies revealed that idealized static structures can not exhibit reasonable ct properties in dna , since dynamical as well as environmental effects were shown to be crucial  @xcite .
dramatic effects are induced by dynamic disorder , for it suppresses ct in homogeneous sequences on the one hand , but can also enhance ct in heterogeneous ( random ) sequences on the other .
moreover , only a minority of conformations appears to be ct - active as has been indicated in ref .
@xcite . therefore , neglecting these significant factors or assuming purely random distributions for dynamical disorder leads to a considerable loss of a vital part of the ct in dna . on this account
the ct properties are evaluated for every snapshot along classical md trajectories which then leads to ensemble averaged quantities , that is the average transmission function @xmath43 and the current @xmath44 .    in fig .
[ fig : compare - dsdna1 ] @xmath43 and @xmath44 are shown for both quadruplex molecules ( with central sodium ions ) , the two central tetrads of ( tg@xmath1t)@xmath1 ( a ) and ( g@xmath45)@xmath6 ( b ) as well as their corresponding double - stranded poly(g ) analogues g@xmath41 and g@xmath45@xcite .
because of the substantial onsite energy fluctuations of about 0.4 ev due to the dna environment and dynamics , the transmission function reveals large broadening for both dna species . the transmission maxima for the quadruplexes
are shifted to lower energies by about 0.3 ev due to the presence of the central sodium ions .
basically , the average transmission is strongly reduced compared to the idealized static structures in fig .
[ fig : g4-trans - ideal ] . nevertheless , the maximum of @xmath43 for the two central tetrads in ( tg@xmath1t)@xmath1 is almost five times larger compared to g@xmath41 .
accordingly , the maximum current is about 4.4 times larger in the quadruplex .        considering the octamers in fig .  [ fig : compare - dsdna1](b ) the conductance difference between g4 and poly(g ) even increases . here , the average transmission in the relevant energy range for the quadruplex is to a great extent two orders of magnitude larger than those for the poly(g ) sequence .
the poly(g ) spectra shows much larger spikes at certain energies .
these spikes indicate the strong impact of dynamics and may be explained by the existence of few charge transfer active conformations , which dominate the average transmission , i.e. few conformations , which feature a high transmission .
this may indicate that even longer sampling ( more than 30 ns ) would be required to converge the spectra .
notwithstanding , it also reveals that in dsdna the average current is dominated to a much larger degree on few non - equilibrium structures , as concluded from our earlier work  @xcite . as a result ,
the average current @xmath44 for ( g@xmath45)@xmath6 is almost two orders of magnitude larger than for g@xmath45 suggesting that the enhancement of ct in g4 with respect to poly(g ) might grow with increasing dna length .    thus , additional sets of ct calculations for longer dna species with 14 and 20 tetrads ( base pairs ) , respectively , are carried out .
the corresponding data was obtained from md simulations of a quadruplex ( g@xmath3)@xmath1 ( including central na@xmath5 ) , and double - stranded dnas g@xmath3 and a heterogeneous sequence containing 31 sites  @xcite recently used by scheer and co - workers in a ct - measurement  @xcite ( note , for ct calculations only the 14 ( 20 ) central sites are used ) .
as expected , the transmission strongly decreases for both dna species .
however , in g4 this effect is not as strong as in dsdna .
for instance the current for ( g@xmath46)@xmath1 is only about 1 order of magnitude lower than for ( g@xmath47)@xmath1 , although expanding to ( g@xmath48)@xmath1 the current drops significantly additional 3 orders of magnitude . in contrast , @xmath44 for poly(g ) and the `` scheer '' sequence decreases by more than 10 orders of magnitude by increasing the number of base pairs from 14 to 20 .
this indicate there is a much stronger distance dependence of ct in dsdna , hence the notion of coherent ct for longer molecular wires might be considerably more likely in g4 than in dsdna .
admittedly , the landauer formalism used in this work performs well for short dna species ( less than 10 sites ) , where the transport is assumed to be at least partially coherent .
on the other hand , it clearly fails for long dna sequences , hence the ct results for the longer molecules should be interpreted only qualitatively and with caution .
for instance , the currents obtained for the 14mer and 20mer of both dsdna molecules poly(g ) and the scheer sequence are orders of magnitude smaller than pa which is far beyond any measurable range .
the complete data is given in fig .
s3 in ref .
@xcite .
the significant conductance difference of g4 and dsdna may not be attributed to the fact that g4 is composed of four poly(g ) like wires . to analyze this further the ct in ( g@xmath47)@xmath1 with central sodium ions ( full md )
is compared to various models , in which i ) only intra- and interstrand@xcite couplings are non - zero ( intra+inter ) and ii ) the quadruplex is separated into its 4 single strands , for which the intrastrand transport is calculated independently and added up afterwards ( @xmath494s ) .
furthermore the results are analyzed with reference to g@xmath45 ( 4xpg ) .
note , for comparison the ct quantities of poly(g ) are multiplied by 4 .    as appears from fig .
[ fig : transport - fullvsintrainter ] , the average transmissions for the two models : full md and intra+inter in the spectral support region clearly reveal the largest plateaus , also showing similar peak structures .
however , both quadruplex models exhibit similar moderate fluctuations . by contrast , the sum of the four single g4 strands ( @xmath494s ) shows much larger fluctuations , comparable with those for g@xmath45 ( 4xpg ) .
moreover , the average transmission for @xmath494s is significantly reduced and is even slightly lower than 4xpg .
interestingly , in @xmath494s are barely ct - active conformations , i.e. single dominating peaks like in 4xpg which sometimes even outreach the maximum transmission of the quadruplex ( full md ) .
this might reflect the structural differences between quadruplexes and dsdna , since the four strands in the quadruplex are not as flexible as those in poly(g ) ( see also rmsd fluctuations in fig .
[ fig : rmsd2 ] ) , hence in poly(g ) the structural phase space is much larger and therefore , several high - transmissive structures can arise . on the other hand ,
the more rigid quadruplex covers only a smaller conformational phase space which indeed ensures a large number of structures showing moderate ct properties for each single strand , respectively .
this clearly indicates , that the most important factor for the enhanced conductance in g4 are the interstrand couplings between the 4 strands in the quadruplex .
thus , if there are conformations in which the 4 isolated channels are not transmissive , there is a considerable probability that ct might occur via coupling between the individual stands . as can be extracted from tab .
[ tab : hop - g4 ] those interstrand couplings are sufficiently large with about 0.01 - 0.02 ev .
as a consequence , at each snapshot there is a substantial amount of pathways over which the ct might occur through the quadruplex .
these findings are supported by the i - v characteristics in fig .
[ fig : transport - fullvsintrainter ] .
the current for the intra+inter model almost matches that of the full md with 2.0 and 2.3 na , respectively .
if we switch off the interstrand couplings in the quadruplex ( @xmath494s ) the current will drop down to 0.08 na which is more than order of magnitude smaller .
a slightly larger current of 0.11 na is obtained for 4xpg .        for further insights into the different ct properties between g4 and dsdna we make use of conformation analysis . for that purpose
, we investigate the amount of conformations which dominate the average ct .
additionally , the distribution of transmissions for the multitude of conformations is explored for both ( g@xmath47)@xmath1 and g@xmath45 .
the results are presented in fig .
[ fig : transport - conf - g4vspg ] .
panel a ) evidently indicates that there are substantially more conformations contributing to the average ct in g4 than in poly(g ) .
consequently , virtually every 10th g4 conformation is ct - active ( about 3000 out of 30,000 ) , whereas only 128 ( again out of 30,000 ) single non - equilibrium structures characterize the average ct in poly(g ) .    )
@xmath1 ( na@xmath5 ) and g@xmath47 ,
a ) : number of conformations that make up 90% of the average transmission maximum @xmath50 , b ) : probability distribution functions ( pdfs ) of @xmath51 for ( g@xmath47)@xmath1 ( na@xmath5 ) and the two models 4xpg and @xmath494s as used in fig .
[ fig : transport - fullvsintrainter ] . ]
second , as demonstrated in the transmission probability distribution functions ( pdfs ) in fig .
[ fig : transport - conf - g4vspg]b ) , the distribution width for the g4 quadruplex is considerably narrower compared to 4xpg and @xmath494s .
note the x - axis is scaled logarithmic .
in addition , the g4 pdf is significantly shifted to higher transmission .
this underscores that the majority of poly(g ) structures is not transmissive , yet a few single conformations are responsible for the average ct . on the other hand , for @xmath494s
the major part of conformations reveal higher transmission than in 4xpg , but the single dominating conformations are missing . clearly , this explains that these single high - transmissive poly(g ) conformations are the reason for a better conductivity in 4xpg compared to @xmath494s .
the higher average t1 values in dsdna are due to few highly conducting conformations .
the smaller average of the g4 t1 couplings therefore resembles a more stable structure , however , not leading to a higher conductance , as could be argued beforehand .
the advantage of g4 over dsdna is due to the existence of non - negligible interstrand couplings in g4 .
the amount of high - transmissive structures is remarkably increased compared to double - stranded poly(g ) dna .    [ [ contact - effects ] ] contact effects + + + + + + + + + + + + + + +    in most of the conductance experiments for dna the molecules are only connected with one strand to the respective left and right contacts .
for instance , in a very recent experiment by scheer _
a single stranded nucleotide with sequence 5-(t@xmath52g@xmath53[ttaggg]@xmath54t@xmath52)-3 ( t@xmath52 denotes modified thymine residues ) , which is known to form stable quadruplexes , has been attached between two contacts  @xcite .
therefore , the question arises whether the ct in an all - parallel stranded quadruplex differs if only 2 or even 1 strand of the quadruplex are coupled to the left and right contacts , respectively . for that purpose , the ct in ( g@xmath47)@xmath1 ( na+ ) is calculated for various contact models : i ) all 4 strands are connected to the left and right electrodes , respectively ( 1 - 4 ) , ii ) only one strand is connected to both contacts ( 1-intra ) and iii - v ) one strand is attached to the left electrode while one of the remaining is contacted to the right one ( 12-inter , etc . ) . as can be seen in fig .
[ fig : transport - contact - ana ] , reducing the number of connected strands from 4 to 1 leads to an decrease in @xmath43 for the relevant energy range of about 1 order of magnitude . however , the transmission for the 1-intra and inter models is very similar indicating a minor significance for which strand or strands are connected to the contacts .
this finding is also supported by the i - v characteristics , for the current ranges around 0.2 na for these models which is roughly one order of magnitude smaller as though all 4 strands are attached ( 1 - 4 model ) .    )
@xmath1 ( na@xmath5 ) : average transmission and current calculated for various contact models i ) all 4 strands are connected to the left and right electrodes , respectively ( 1 - 4 ) ii ) only one strand is connected to both contacts ( 1-intra ) , iii - v ) one strand is calculated to the left electrode while one of the remaining is contacted to the right one ( 12-inter , etc . ) , and vi ) comparison with poly(g ) . ]
interestingly , there seems to be an increase in the transmission fluctuations if the quadruplex is contacted through only one strand at each end .
notwithstanding , the transmission plateau for the 1-stranded ( intra and inter ) contacted models is still about 1.5 orders of magnitude larger compared to those for poly(g ) resulting in an average current which is again 1 order of magnitude smaller with 0.025 na . thus our results suggest that independent on the various contact linking schemes of g4 and dsdna one might expect a higher conductivity for the quadruplex .
certainly , the optimal conductance for all - parallel stranded g4-dna is ensured if all four strands are coupled to the contacts .
[ [ effect - of - dna - environment - on - transport ] ] effect of dna environment on transport + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    as discussed above , the major part of the dynamical disorder is induced by the qm / mm environment , i.e. the last term in eq .
[ eq : qmmm ] which is built of the mm charges of dna backbone , solvent and counterions .
previous results have indicated that the disorder due to the dna environment might not only suppress ct in homogeneous sequences like poly(g ) , rather it is able to enhance ct in random sequences like the heterogeneous dickerson dodecamer  @xcite .
very recent experiments by scheer _
_ confirmed this notion , since they found the current for a 31 base pair sequence  @xcite to be 2 orders of magnitude smaller in vacuo than in aqueous solution  @xcite .
therefore , it might be interesting to investigate , what happens if we switch of the environment for ct in g4-dna , also with respect to the effect in dsdna . as can be seen in fig .
[ fig : transport - envvsvacuo](a ) , the transmission maximum for ( g@xmath47)@xmath1 in vacuo is about 2 orders of magnitude larger than with the qm / mm environment . besides , its broadening and also the fluctuations are significantly reduced and the plateau is shifted to higher energies due to the neglect of the electrostatic interaction with the sodium ions within the quadruplex .
the transmission for poly(g ) in vacuo shows basically the same features , for the broadening is likewise strongly reduced and the maximum is located in the same energy range as for g4 in vacuo , although it is 2 order of magnitude smaller . in general , the transmission for g4 including the qm / mm environment clearly reveals the largest broadening which might indicate that the central sodium ions have an additional strong influence on the dynamic disorder due to longitudinal mobility within the quadruplex , which is not the case in dsdna . despite all that , the current at high voltages is larger for g4 with qm / mm than for poly(g ) in vacuo
.        basically , for heterogeneous sequences ( not shown here ) a reduced transmission in vacuo is found , which is caused by energy gaps between a and g states .
notwithstanding , for dna molecules with no static energy gaps like double - stranded poly(g ) and g4-dna ( both with uniform dna bases ) , the qm / mm environment is most likely to increase the dynamical disorder , thus will suppress ct compared to the vacuo model . as a result , there is no significant difference in the effect of the dna environment on ct for g4 and poly(g ) dna .
this is also underscored by conformation analysis given in panel b ) and c ) in fig .
[ fig : transport - envvsvacuo ] , which indicate that there are considerably more ct - active conformations in vacuo than with qm / mm environment for both dna species g4 and poly(g ) . moreover , for g4 in vacuo nearly every conformation appears to be high - transmissive , since the average maximum transmission @xmath50 increases almost linearly with the number of conformations .
thus , the ct in vacuo is only marginally affected by single non - equilibrium conformations , rather the whole ensemble of g4 conformations seems to be ct active which is demonstrated in the pdf of transmission .
one last issue remains when considering ct in g4-dna that is the effect of ions within the quadruplex channel .
the structural influence of these ions has already been addressed in detail in sec .
[ subsec : res - md ] . as could be seen in sec .
sii.b in ref .
@xcite we did not find the central ions to contribute states in the relevant energy range for hole transfer in g4 , therefore , the effect will only be investigated electrostatically . in fig .
[ fig : transport - ion - ch ] the average transmission and current is shown for the quadruplex simulations of ( g@xmath47)@xmath1 in absence and presence of either lithium , sodium and potassium ions .
obviously , the transmission maximum is only marginally affected by the presence of different types of central ions .
furthermore , as expected the transmission function is shifted to lower energies if central ions are present . however , for potassium the transmission is to be found slightly reduced ,
also the broadening is not as large as for the other species .
this is also reflected in the pdf of transmission maxima which can be found in fig .
s4 in ref .
@xcite .    )
@xmath1 in absence and presence of monovalent central ions li@xmath5 , na@xmath5 and k@xmath5 . ]    as a consequence , the average current for the simulation with central potassium ions is half as large as for those with lithium and sodium ions which might be attributed to the different mobilities of li@xmath5 and na@xmath5 compared to k@xmath5 . interestingly , there is no significant difference for ct in absence and presence of ions , although it is known from sec .
[ subsec : res - md ] that the g4 molecules without central ions exhibit significant destabilization .
once more this supports the notion that the enhanced conductance in g4 may not exclusively be explained in terms of higher structural rigidity , i.e. less dynamical disorder , rather it is the multitude of ct pathways via interstrand couplings that brings on an increased number of high - transmissive conformations
. apparently , those interstrand and in - plane couplings t2 and t3 are not altered by the presence of central ions , i.e. by a more rigid quadruplex structure .
in this work , we have investigated the conductivity of g4 with respect to dsdna using classical md simulations , combined qm / mm methods to compute ct parameters and landauer theory to compute transmission and i - u characteristics .
these properties have been evaluated for phase - space ensembles of g4 and dsdna structures .
the approach adopted here to calculate the transport characteristics at a given time ( landauer theory ) has clear limitations related to the underlying assumption of tunneling transport .
though this may be an efficient pathway over very short segments of the molecules under study , its validity becomes questionable with increasing length .
this becomes quite evident when investigating sequences as used in recent experiments e.g. by scheer and co - workers  @xcite , where our methodology would predict currents order of magnitude smaller than those found in the experiments .
we are currently developing a methodology in order to describe transport in complex materials in a more general way  @xcite .    nevertheless , the landauer calculations reported here give interesting insight into particular properties of g4 , especially when compared to dsdna .
first of all , g4 with central ions is structurally much more stable than dsdna and has more @xmath0-contacts along the chain @xcite . at first sight
, one may argue that this leads to improved electrical conduction , however , it turns out that the large fluctuations found for dsdna lead to highly conducting structures , which dominate the transport . therefore , a dsdna molecule conducts better than one strand of g4 , or equivalently , the higher conduction of g4 is not due an increased structural stability of its single strands .
the phase space of the single strands still contains a vast majority of conformations , which are not ct active .
it is the ability of g4 to allow for a large number of conformations due to the interstrand couplings t2 and t3 which makes this species better conducting .
if the pathway along one strand is blocked , e.g. when one t@xmath40 along the chain vanishes , many other conduction channels may be viable due to interstrand hopping . at the end , g4 has a much larger number of ct active conformations than four dsdna ( poly(g ) ) offers .
this is the basis of the higher conductivity of g4 .
this advantage is even maintained , when contacted only at one g - site , instead of four sites .    as a second point
, dsdna or g4 may be exerted to strain due to the contacting procedure . here , clearly the higher stability of g4 may help to maintain a conducting conformation , while the conduction in dsdna may be much more easily disrupted  @xcite .
in a recent work , we have predicted dsdna with non - uniform sequences ( i.e. not poly(g ) or poly(a ) ) to conduct better in solution than in gas phase .
this is due to the larger structural fluctuations in solvent , which introduce conducting structures  @xcite .
this finding has been confirmed by a recent experiment  @xcite . for homogeneous sequences like poly(g ) or g4 however
, a higher conductance in gas phase should be found  @xcite .
finally , our results suggest that the presence of central metal ions within the quadruplex has only a marginal impact on hole transfer in g4 , although they are vital for the stability and rigidity of g4-dna .
this work has been supported by the deutsche forschungsgemeinschaft ( dfg ) within the priority program 1243 `` quantum transport at the molecular scale '' under contract el 206/5 - 2 and cu 44/5 - 2 , by the volkswagen foundation grant nr
.  i/78 - 340 and by the european union under contract ist-029192 - 2 .
note , generally , the two last tetrads ( base pairs ) at the 5 and 3 end are not considered for ct calculations to avoid end effects , although the simulations were done with 12 tetrads and base pairs for g4 and dsdna , respectively .