2003
Доклады БГУИР
АПРЕЛЬ-ИЮНЬ
Том 1, № 2
УДК 621.38.049.77.002.3
ВОЗМОЖНОСТИ ПРИМЕНЕНИЯ МОЛЕКУЛЫ ДНК В КАЧЕСТВЕ ПЕРЕКЛЮЧАЮЩЕГО ЭЛЕМЕНТА
Н.В. НОВИК, Ю.А. БЕРАШЕВИЧ, В.Е. БОРИСЕНКО
Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь
Поступила в редакцию 29 мая 2003
Предложена электрическая эквивалентная схема молекулы ДНК, которая описывает перенос носителей заряда через молекулу и состоит из параллельно включенных туннельных и выпрямляющих диодов, проводимость которых определяется скоростью переноса носителей заряда через различные последовательности динуклеотидов посредством туннельного и термоактивированного механизмов переноса соответственно. Показана возможность создания приборных структур на основе молекулы ДНК, выполняющих логические операции инверсии, дизъюнкции и конъюнкции. Преимуществами таких приборов являются возможность регулирования проводимости молекулы с помощью изменения набора образующих ее нуклеотидных оснований, высокая скорость протекания процессов переноса, порядка 10-1410-12 с, широкий диапазон рабочих температур.
Ключевые слова: молекула ДНК, транспорт носителей заряда, переходные характеристики
Введение
Научно-техническое направление, связанное с исследованием и возможным применением молекул ДНК для реализации новых сверхбыстродействующих устройств электронной обработки информации, в последние годы интенсивно развивается и во многих отраслях современной науки в настоящее время занимает ключевые позиции [1, 2]. Это обусловлено большой перспективностью таких структур для изготовления уникальных наноэлектронных приборов и возможностью их использования в качестве модельных объектов для изучения поведения различных биологических систем.
Согласно экспериментальным данным [1-3], молекула ДНК любой конфигурации отличается возможностью совмещения с различными материалами благодаря имеющимся в ее полимерной цепи химически активным концам. В перспективе это позволит помещать ее как между металлическими контактами, так и между полупроводниковыми контактами. Для использования структуры молекулы ДНК в качестве активного переключающего элемента в твердотельных приборах возникает необходимость создания адекватного аппарата моделирования характеристик проектируемого устройства в составе интегральных схем. Одной из важнейших составляющих такого конструкторского инструментария является разработка эквивалентной схемы, построенной на базе молекулы ДНК. Кроме того, функция базовой приборной структуры, заключающаяся в переключении двух основных состояний, соответствующих высокому и низкому уровням сигнала, т. е. высокопроводящих участков молекулы на низкопроводящие, может обеспечить создание на ее основе логических переключающих элементов, однако эти возможности не изучены в полной мере.
Целью данной работы является разработка эквивалентной схемы молекулы ДНК и логических элементов, пригодных для использования в стандартных пакетах схемотехнического моделирования типа Р8рюе.
Эквивалентная электрическая схема молекулы ДНК
Полинуклеотидная последовательность в молекуле ДНК может быть представлена в виде периодической наноразмерной структуры, в которой за счет различного энергетического положения нуклеотидных оснований образуется система потенциальных барьеров и ям. При этом тимин-адениновые (А-Т) основания выполняют функцию потенциального барьера как для дырок, так и для электронов, а гуанин-цитозиновые (О-С) — потенциальной ямы (рис. 1) [4].
-Е,
Е*
[©рррр
С
С
0
С
С
С
С
Е*
Рис. 1. Энергетическая диаграмма и механизмы переноса дырок в 5'-О-(Л)4-ООО-3' молекуле ДНК: 1 — туннельный механизм переноса, 2 — термоактивированный перенос носителей заряда
Из результатов экспериментальных данных следует [3, 5, 6], что перенос дырок в молекуле ДНК является более эффективным по сравнению с электронным переносом, что связано с иерархией потенциалов окисления нуклеотидных оснований (О < Л < С, Т). Низкоэнергетическое расположение азотистого основания гуанина О обеспечивает наилучшие условия для захвата, последующего возбуждения и перемещения дырки от одного нуклеотидного основания к другому, что открывает возможность для наблюдения транспорта носителей вдоль молекулы ДНК, в то время как захват и возбуждение отрицательного носителя заряда в молекуле затруднены из-за подобия потенциалов окисления тимина и цитозина [7]. При этом при образовании дуплетов ОО и триплетов ООО в структуре молекулы ДНК увеличивается вероятность захвата дырки гуанином за счет изменения его потенциала окисления [6].
Рассмотрим возможные механизмы переноса носителей заряда в молекуле ДНК. В случае образования потенциального барьера, сформированного нуклеотидной А-Т парой, расположенной между двумя парами О-С динуклеотидов, возбужденная дырка может перескакивать или туннелировать через потенциальный барьер, как представлено на рис. 1. Увеличение числа Л-Т пар ведет к росту ширины потенциального барьера. В этом случае вероятность туннелиро-вания экспоненциально уменьшается, а термоактивированный перенос осуществляется в несколько этапов. Первый этап включает термоактивированный захват дырки первой Л-Т парой, последующий электронный суперобмен между двумя одинаковыми парами и переход дырки в О-С пару (рис. 1). При этом ввиду перекрывания п-электронных орбиталей нуклеотидных оснований, которое имеет наибольшее влияние в случае идентично расположенных динуклеоти-дов вдоль полимерной цепи ДНК, предполагается сохранение энергии и импульса переносимых таким образом носителей заряда [8]. Механизм такого транспорта носителей заряда был описан в рамках туннельной суперобменной теории электронного переноса [8].
Изменение пространственного расположения и числа нуклеотидных остатков в спирали молекулы ДНК приводит к варьированию условий, определяющих механизм переноса в различных молекулярных последовательностях, и влияет как на величину скорости протекания процессов переноса, так и на время перемещения носителей заряда по полимеру. Положение донорного участка молекулы на более высоком энергетическом уровне относительно уровня акцепторной области Еа обеспечивает наиболее выгодные условия для перемещения возбужденного носителя заряда через углеводно-фосфорный скелет ДНК, т. е. при выполнении условия Еа > Еа. Такое условие соблюдается в 5-О-(Л)п-ООО-3' конфигурации молекулы ДНК.
Увеличение числа последовательно соединенных пар гуанина за счет перекрывания электронных орбиталей ведет к появлению энергетического зазора между 5'-О участком и ООО-З'. Вследствие существования этого энергетического зазора между О-нуклеотидом и триплетом ООО, равного порядка Еа — Еооа =0,68 эВ [8], появляются идеальные условия для переноса
возбужденных дырок от донора, в качестве которого выступает 5-О, к акцептору. При этом перенос носителей заряда в последовательно соединенных идентичных участках молекулы осуществляется за счет делокализации п-электронов при перекрывании их электронных орби-талей, имеющих возможность свободно перемещаться вдоль колец нуклеотидов.
Для моделирования протекающих в структуре молекулы ДНК механизмов переноса примем следующие допущения. Перенос носителей заряда в направлении действия приложенного внешнего смещения значительно выше, чем в противоположном направлении, перенос носителей заряда через потенциальный барьер рассматривается как одномерный процесс, а составляющие импульса носителя, простирающиеся в других направлениях, перпендикулярных направлению переноса, задаются как фиксированные параметры.
Исходя из принятых предположений, темп носителей заряда посредством термоактивированного механизма переноса в приближении термодиффузионной теории определяется как
кьор = 4тп*к2вТ2 (
И3
1 — ехр
V аквТ У
(и
УУ
| ехр
(
ЛЕ — У д'УъгаХ8^™)-1
ЫаБ'
квТ
л ^
dx
(1)
У У
г-г *
где И — постоянная Планка; т е — эффективная масса электрона; кв — постоянная Больцма-на; Т — температура; а — число нуклеотидных оснований в молекуле ДНК; п — число нуклеотидных оснований, образующих потенциальный барьер; ЛЕ х квТ — высота потенциального барьера; е — диэлектрическая проницаемость барьера; УЬга!1 — внешнее приложенное напряжение.
Темп туннельного переноса можно описать в рамках аппроксимации Вентцеля-Крамерса-Бриллюэна [9]:
4*0 0
ки = Лт_[[) — р(е + дУы<аа/п)]Е [Б(Е
И3d
(2)
ЛЕ
где d = аЖ — длина молекулы; Р(Е) — функция распределения Ферми-Дирака; D(E) — вероятность туннелирования, которая равна:
Б = ехр
( г г
2т*
—- [ ЛЕ —
И *
х1 V - -
д3п¥ъ
Ьгаз
а
4п:
1/2
— Е + Е,
1/2
Л
dx
(3)
у
где Е — энергия туннелирования электронов; х1, х2 — классические точки поворота; Е, — энергия электрона, связанная с поперечной составляющей импульса.
Для учета взаимосвязи явлений на всех периодах структуры ДНК рассмотрим кинетику изменения заряда в г-й квантовой яме структуры, описываемую следующей системой уравнений:
—1
Е
dPi _ G -
- к^°Р — ktm
dt
dp- _ kh°p + k—l — kh°p — k;u", (4)
dt ' i i i
dpN _ k^°P + ktun _ k^°P _ ktu"
j — kN—l + kN—l kN kN ' dt
где ' — период структуры (i—1...N); Gp — скорость возбуждения дырки; pi — концентрация
l h°p i tun «
дырок в i-яме; ki , ki — темп термоактивированного и туннельного переноса носителем заряда через '-потенциальный барьер, сформированный набором А-Т пар. Для описания переноса электронов через полинуклеотидную структуру молекулы можно построить аналогичную систему кинетических уравнений.
Для моделирования переходных процессов, протекающих в периодической структуре молекулы ДНК, находим решение системы кинетических уравнений (4) в нестационарных условиях. Основным параметром при этом является время переключения приборной структуры (рис. 1) из одного состояния в другое, что позволяет определить быстродействие этого прибора и моделировать работу переключающего элемента в импульсном режиме. Так, за время распространения импульса напряжения tl захваченные А-Т парами дырки накапливаются в ямах, при этом уравнение, описывающее процесс накопления в нестационарных условиях, имеет вид
p,(ti) _ Pi—l(tl) 'т' (l — exp(—11 /T')) + p(t _ 0) ■ exp(—1, /zj, (5)
T—l
_ h°p _ tun
T • T
i i
T' _ T h°p + Ttun ' (6)
h°p tun „
где Tj ' Tj определяются из темпов переноса носителей заряда [7J. При этом первый член
уравнения описывает процесс накопления носителей заряда в -яме за счет прихода носителей из (i -1) ямы. Второй член уравнения учитывает рассасывание носителей заряда в i -яме и имеет место в случае ненулевой концентрации носителей заряда в ямах при подаче импульса напряжения.
За время задержки импульса t2 происходит рассасывание дырок, в течение которого импеданс переключающего элемента сохраняется низким и выходной ток определяется внешней цепью. При этом темп генерации дырок на инжектирующем контакте принимается равным нулю. В этом случае решение системы кинетических уравнений (4) для процесса рассасывания на всех периодах структуры молекулы ДНК в нестационарных условиях имеет вид
P(t2) _ PiМТ' (l — exp(—2 / Т )) + P(t _ 0)exp(—t2 / г,). (7)
T—l
Причем для первой ямы, так как генерация носителей заряда на инжектирующем контакте принимается равной нулю, это уравнение выглядит следующим образом:
Pl(t2) _ Pl(t _ 0)exp(—t2/rt). (8)
Предложенная кинетическая модель переноса носителей заряда через систему потенциальных ям и барьеров описывает транспорт носителей заряда в молекуле ДНК через темп перехода носителей в i-ю яму и темп ухода из этой ямы посредством упругого туннелирования и термоактивированного переноса заряда через потенциальные барьеры.
Помещение молекулы ДНК между двумя металлическими контактами и приложение к ее структуре внешнего электрического поля увеличивает смещение энергетических зон донор-ного и акцепторного уровней в контактных областях молекулы. В результате выполнения условия Ел > Еа возникает беспрепятственный переход носителей заряда от одной квантовой ямы к другой посредством двух упомянутых выше механизмов, что можно охарактеризовать как прямое включение молекулы ДНК.
Иначе, в случае обратного смещения происходит уменьшение энергетического зазора между донором и акцептором вплоть до появления условия, когда Ел < Еа. Соответственно при небольшом обратном смещении, когда условие Ел < Еа еще не достаточно, туннелирование носителей заряда в молекуле ДНК отсутствует, а высокий потенциальный барьер между Еа и валентной зоной А-Т пар уменьшает вероятность переноса носителей заряда по термоактивированному механизму. При достаточно большом внешнем напряжении (УЫаб. > Еа - Ел) появляется компонента туннельного переноса, и общий перенос носителей заряда через структуру экспоненциально возрастает с увеличением внешнего смещения. Однако скорость экспоненциального роста ниже, чем для прямого смещения, за счет разницы потенциальных барьеров в структуре для возбужденных носителей заряда в двух этих случаях. Таким образом, поведение молекулы ДНК при приложении к ней внешнего напряжения Уыа. определенной полярности и величины можно положить в основу работы переключающего элемента.
На основе описанных выше закономерностей электронного переноса носителей заряда в молекуле ДНК была предложена следующая приборная структура на основе 5'-0-(Л)3-000-3' последовательности молекулы ДНК (рис. 2). Схема переключающего элемента включает три параллельных диода: туннельный диод Б1, в основе работы которого лежит практически независимый от температуры туннельный перенос носителей заряда, и два выпрямляющих диода Б2 и Б3, характеризующихся экспоненциальной зависимостью тока от температуры. Второй выпрямляющий диод Б3 учитывает дополнительный процесс термоактивированного переноса носителей заряда, имеющий место в случае изменения конфигурации молекулы ДНК, например в 5'-0-(Л-Т-Л-А)-000-3' последовательности. Возникающий в этом случае дополнительный потенциальный барьер, сформированный тимином за счет различного энергетического положения нуклеотидных оснований, нарушает условия электронного переноса по суперобменному механизму и снижает общий перенос дырок в пределах А-Т-А участка. Однако появление дополнительной компоненты переноса носителей по термоактивированному механизму через сформированный тимином потенциальный барьер компенсирует снижение тока, протекающего через приборную структуру.
D1
i
Шп
bias
ч D3 У • hop i ■3 ,
R
N
Рис. 2. Эквивалентная схема переключающего элемента, разработанного на базе 5'-в-(Л)3-ввО-3' последовательности молекулы ДНК
Результаты расчета и их обсуждение
Согласно экспериментальным данным [3, 8, 10], скорость переноса носителей заряда в молекуле ДНК определяется энергетическим смещением донорного и акцепторного участков относительно друг друга при воздействии на систему внешнего электрического поля, временем накопления заряда в молекуле во время протекания прямого тока, а после переключения рассматриваемой системы - временем рассасывания заряда при протекании обратного тока, а также температурой. Поэтому моделирование эквивалентной схемы молекулы ДНК и ее переключающих характеристик проводили с учетом перечисленных факторов. Параметры, использованные в расчетах, приведены в таблице.
Параметры, используемые для моделирования эквивалентной схемы 5'-С-(Л)п-ООС-3' структуры молекулы ДНК
Параметр Значение
Расстояние между соседними нуклеотидными основаниями Высота барьера для дырок между О-Л нуклеотидами Высота барьера для дырок между (О)2-Л нуклеотидами Высота барьера для дырок между (О)3-Л нуклеотидами Число динуклеотидов в молекуле ДНК Число динуклеотидов, образующих потенциальный барьер Эффективная масса дырки Относительная диэлектрическая проницаемость молекулы Температура 0,34 нм [8] 0,5 эВ [4] 0,9 эВ [11] 1,1 эВ [11] 5-12 1-7 * 0,9 те 4 [5] 50-350 К
Результаты моделирования ВАХ переключающего элемента на базе однопериодных 5'-О-(Л)п-ООО-3' и 5'-ОО-(Л)п-ООО-3' молекул ДНК приведены на рис. 3. Очевидно, что увеличение числа А-Т динуклеотидов ведет к пропорциональному изменению как прямого, так и обратного тока, в то время как изменение конфигурации 5'-О-(Л)п-ООО-3' молекулы ДНК на 5'-ОО-(Л)п-ООО-3' уменьшает различие между этими токами на два порядка.
10"т
1010'
109-
108-
107-
10й-
105-
104-
103-
102-
Рис. 3. ВАХ переключающего элемента, разработанного на базе 5'-в-(Л)3-ООО-3' и 5'-00-(Л)з-00в-3' последовательностей молекулы ДНК
Моделирование переходных процессов в приборной структуре, построенной на базе двухбарьерной 5'-О-(Л)п-О-(Л)п-ООО-3' молекулы ДНК, показало, что основными параметрами, определяющими переключающие свойства подобного прибора, являются время накопления и время рассасывания заряда (дырок) на определенном периоде полинуклеотида. Результаты моделирования переходных процессов в двухбарьерной 5'-О-(Л)2-О-(Л)2-ООО-3' структуре молекулы ДНК представлены на рис. 4. Время накопления дырок в первой потенциальной яме, зависящее от положения первой О-С пары оснований относительно дна квантовой ямы и А-Т динуклеотида, лежит в фемтосекундном диапазоне и составляет 1 пс, причем во второй потенциальной яме процесс накопления носителей заряда происходит с некоторой задержкой (кривая 3, рис. 4). При этом рассасывание дырок во второй потенциальной яме, а также вблизи области акцепторного участка молекулы, сформированного триплетом гуанина ООО, происходит с учетом вклада концентрации дырок в предыдущих потенциальных ямах, о чем свидетельствуют кривые 2, 4 (рис. 4). Кроме того, уменьшение количества А-Т пар, приводящее к сужению потенциальных барьеров в структуре молекулы ДНК, уменьшает время переключения приборной структуры в два раза.
1:, пс
Рис. 4. Зависимость концентрации дырок от времени действия приложенного смещения в двухбарьерной 5'-О-(Л)2-в-(Л)2-ООО-3' молекуле ДНК: 1, 2 — накопление и рассасывание дырок в первой потенциальной яме соответственно; 3, 4 — накопление и рассасывание дырок в контактной области ООО триплета соответственно
Темп термоактивированного переноса носителей через потенциальный барьер экспоненциально снижается с уменьшением температуры, в то время как туннельный перенос носителей практически не зависит от температуры. Поэтому на графике зависимости плотности протекающего в 5'-О-(Л)3-ООО-3' молекуле ДНК тока от температуры (рис. 5) наблюдается участок перекрытия термоактивированного механизма с туннельным, который с увеличением напряжения сдвигается в область более высоких температур. Уменьшение числа динуклеоти-дов ведет к превалированию туннельного переноса носителей заряда в широком диапазоне внешних смещений и соответственно к стабильности тока при изменении температуры. В свою очередь, увеличение числа п приводит к росту вклада термоактивированного переноса и к его очевидному преобладанию над туннельным механизмом переноса при температуре свыше 200 К.
т, к
Рис. 5. Зависимость плотности протекающего в 5'-0-(Л)3-000-3' молекуле ДНК тока от температуры
Таким образом, перспектива создания переключающих элементов, построенных на базе молекулы ДНК, очевидна. Это обусловлено в первую очередь высоким быстродействием таких структур, которое создает предпосылки для разработки новых сверхбыстродействующих нано-электронных приборов. Кроме того, существует возможность уменьшения времени переключения приборных структур на основе молекулы ДНК до 0,1 пс за счет использования частичного насыщения и рассасывания заряда в них. При этом роль температуры в проектируемых устройствах на основе молекулы ДНК достаточно велика, так как диапазон ее значений определяет превалирование того или иного механизма переноса в полинуклеотиде молекулы и может обеспечить адекватный выбор конфигурации ДНК для получения наилучших переключающих свойств биоэлектронного устройства в целом.
На основе переключающих свойств однопериодной структуры 5'-0-(Л)п-000-3' молекулы ДНК с несимметричной вольт-амперной характеристикой, определяемой полярностью внешнего смещения (рис. 3), схемотехнически могут быть реализованы цифровые цепочки для интегральных схем. Так, предложенная приборная структура (рис. 2) с двумя устойчивыми состояниями может использоваться в логических схемах как элемент, выполняющий операцию инверсии, так как при подаче на его вход смещения высокого уровня У1=0,8^1,3 В, соответствующего прямому включению, на выходе прибора имеет место низкий уровень сигнала, а при подаче на этот же вход отрицательного смещения низкого уровня У0= -0,6^- -1 В, соответствующего обратному смещению — высокий уровень сигнала (рис. 3). В рамках стандартной полупроводниковой диодно-резисторной логики операция сложения на основе приборной структуры реализуется при параллельном включении двух структур, а при их антипараллельном включении реализуется операция умножения. Кроме того, выбором числа последовательностей в молекуле ДНК можно регулировать не только быстродействие переключения данных элементов, но и зависимость их электрических характеристик от температуры за счет различной зависимости темпа термоактивированного и туннельного переноса от ширины потенциального барьера.
Заключение
Разработана эквивалентная электрическая схема приборной структуры на базе молекулы ДНК, состоящая из параллельно включенных туннельных и выпрямляющих диодов, проводимость которых определяется скоростью переноса носителей заряда через различные последовательности динуклеотидов посредством туннельного и термоактивированного механизмов переноса соответственно.
Преимуществами приборов на основе молекулы ДНК являются возможность регулирования проводимости молекулы с помощью изменения набора образующих ее динуклеотидов, высокая скорость протекания процессов переноса, порядка 10-14-10-1 2 с, низкая потребляемая мощность, отсутствие реактивного сопротивления, широкий диапазон рабочих температур, а также простота разработки переключающих и логических элементов для наноэлектроники на ее основе.
APPLICATION OPPORTUNITIES OF DNA MOLECULE FOR THE SWITCHING ELEMENT
N.V. NOVIK, J.A. BERASHEVICH, V.E. BORISENKO Abstract
The electric equivalent circuit of molecule DNA, which describes charge carrier transport via molecule, is proposed. The circuit consists of parallel-connected tunnel and rectifying diodes. Their conductivity is determined by rate of charge carrier transport composed of activated and tunneling mechanism and depended on prevalence of these mechanisms in the case of various sequences of nu-cleotides. Based on DNA molecule, the opportunities of elaboration of logical devices performed the inversion, disjunction and conjunction operations, are considered. It is shown that advantages of such devices are possibility to regulate a conductivity of a molecule in dependence on molecular sequences, high rate of transport processes with the switching time of 10-14-10-12 s and operation in a wide range of temperatures.
Список литературы
1. Collier C. P., Jeppesen J. O., Luo Y., Perkins J., Wong E. W., Heath J. R., Stoddart J. F. // J. Am. Chem. Soc. 2001. Vol. 123. P. 12632-12641.
2. Storm A. J., Noort J., Vries S., Dekker C. // Appl. Phys. Lett. 2002. Vol. 79, № 23. P. 3881-3883.
3. Giese B., Amaudrut J., Kothrin A.-K., Spormann M., Wessely S. // Nature. 2001. Vol. 412. P. 318-320.
4. Berashevich J. A., Filonov A. B., Borisenko V. E. // Physics, Chemistry and Application of Nanostructures. Singapore: World Scientific. 2003. P. 235-238.
5. Bixon M., Giese B., Wessely S., Langenbacher T., Michel-Beyerle M. E., Jortner J. // PNAS. 1999. Vol. 96, № 21. P. 11713-11716.
6. Meggers E, Kusch D, Spichty M, Wille U, Giese B. // Angew. Chem. Ent. Ed. 1998. Vol. 37. P. 460-461.
7. Берашевич Ю. А. // Изв. Белорус. инж. акад. 2002. № 2 (14)/2. С.5-7.
8. Troissi A., Orlandi G. // Chem. Phys. Lett. 2001. № 344. P. 509-518.
9. Duke C. B. In: Tunneling in solids. Solid physics supplement // New York: Academic Press, 1969. Vol. 10.
10. Wan C, Fiebig T, Kelley S. O, Treadway C. R., Barton J. K., Zewail A. H. // Proc. Natl. Acad. Sci. USA. 1999. Vol. 96. P. 6014-6019.
11. Sugiyama H, Saito I. // J. Am. Chem. Soc. 1996. № 188. P. 7063-7068.