Độ dịch chuyển hóa học NMR là biến thiên tương đối giữa tần số cộng hưởng của một hạt nhân từ tính cụ thể so với tần số cộng hưởng của chất chuẩn đối chiếu quy ước khi đặt trong cùng một từ trường ngoài, phản ánh mật độ phân bố electron và môi trường liên kết vi mô bao quanh hạt nhân đó. Đại lượng không thứ nguyên này được biểu diễn bằng đơn vị phần triệu và là thông số quan trọng nhất trong phổ cộng hưởng từ hạt nhân giúp định danh nhóm chức, xác định cấu trúc phân tử, lập bản đồ tương tác phối tử và phân tích động học phân tử. Bài viết trình bày toàn diện về cơ sở vật lý lượng tử, biểu thức toán học, quy ước chất chuẩn quốc tế, các khoảng dịch chuyển đặc trưng, phương pháp mô phỏng tính toán cùng các ứng dụng kiểm nghiệm thực tế.
Cơ sở vật lý và cơ chế chắn từ hạt nhân
Hiện tượng cộng hưởng từ hạt nhân xuất hiện khi các hạt nhân mang spin từ tính khác không được đặt trong một từ trường tĩnh bên ngoài. Tuy nhiên, hạt nhân trong phân tử không chịu tác động trực tiếp của từ trường ngoài trần mà được bao bọc bởi một đám mây electron chuyển động liên tục. Dưới tác dụng của từ trường ngoài, các electron này sinh ra một dòng điện cảm ứng nhỏ, từ đó tạo ra một từ trường thứ cấp cục bộ có hướng chống lại hoặc tăng cường từ trường ngoài.
Tương tác che chở này được mô tả thông qua đại lượng hằng số chắn, ký hiệu là . Hằng số chắn là một đại lượng vô hướng trong dung dịch đẳng hướng thu được từ đường chéo của tensor chắn từ hạt nhân, như được IUPAC chuẩn hóa trong khuyến nghị bổ sung năm 2008 của Harris và các cộng sự. Từ trường hiệu dụng mà hạt nhân thực sự cảm nhận được, ký hiệu là , có mối liên hệ với từ trường ngoài, ký hiệu là , theo biểu thức:
Trong đó là cường độ từ trường hiệu dụng tại tâm hạt nhân, là cường độ từ trường ngoài do nam châm siêu dẫn của thiết bị tạo ra, và là hằng số chắn từ không thứ nguyên của môi trường electron bao quanh. Khi mật độ electron xung quanh hạt nhân tăng cao, giá trị tăng lên, hạt nhân được gọi là bị chắn mạnh hơn. Khi đó từ trường hiệu dụng giảm xuống, khiến tần số cộng hưởng của hạt nhân dịch chuyển về phía tần số thấp hơn, theo thuật ngữ phổ học hiện đại gọi là dịch chuyển về vùng trường mạnh (upfield). Ngược lại, khi các nhóm hút electron làm giảm mật độ điện tích quanh hạt nhân, giá trị suy giảm, hạt nhân bị giải chắn và cộng hưởng tại tần số cao hơn, tương ứng với vùng trường yếu (downfield).
Đóng góp nghịch từ và thuận từ vào hằng số chắn
Về mặt cơ học lượng tử, hằng số chắn toàn phần được phân tách thành hai thành phần chính theo lý thuyết Ramsey:
- Thành phần chắn nghịch từ: Ký hiệu là , sinh ra từ chuyển động quay cục bộ của các electron ở trạng thái cơ bản trong các orbital cầu đối xứng quanh hạt nhân. Thành phần này luôn mang giá trị dương, có xu hướng chống lại từ trường ngoài và đóng vai trò chi phối tuyệt đối đối với các hạt nhân nhẹ như proton.
- Thành phần chắn thuận từ: Ký hiệu là , xuất hiện do sự trộn lẫn của các trạng thái electron kích thích vào trạng thái cơ bản dưới tác động của từ trường khi các orbital electron không có tính đối xứng cầu, đặc biệt là các orbital p và d. Thành phần này thường mang giá trị âm, có tác dụng khử chắn và đóng vai trò quyết định độ dịch chuyển của các hạt nhân nặng hơn như carbon hay dị tố.
Định nghĩa toán học và thang đo ppm
Tần số cộng hưởng Larmor của một hạt nhân tỉ lệ thuận với cường độ từ trường ngoài đặt vào máy đo. Nếu báo cáo vị trí hấp thụ của hạt nhân bằng đơn vị tần số tuyệt đối là hertz, giá trị này sẽ thay đổi tùy thuộc vào cường độ từ trường của từng máy quang phổ cụ thể, gây khó khăn cho việc so sánh kết quả giữa các phòng thí nghiệm sử dụng thiết bị có công suất khác nhau. Để loại bỏ sự phụ thuộc vào cường độ nam châm, cộng đồng hóa học quốc tế đã thiết lập khái niệm độ dịch chuyển hóa học tương đối.
Theo bản khuyến nghị danh pháp chính thức của IUPAC công bố năm 2001 do Harris và các cộng sự chủ trì biên soạn, độ dịch chuyển hóa học, ký hiệu bằng chữ cái Hy Lạp , được xác định theo công thức sau:
Trong công thức trên, đại lượng là độ dịch chuyển hóa học tính bằng phần triệu; là tần số cộng hưởng quan sát được của hạt nhân khảo sát tính bằng hertz; là tần số cộng hưởng của chất chuẩn đối chiếu đo trong cùng điều kiện tính bằng hertz; và là tần số hoạt động danh định của máy đo phổ đối với hạt nhân đó tính bằng megahertz. Nhờ nhân với thừa số mười mũ sáu, giá trị trở thành một hằng số vật lý nội tại của phân tử, độc lập với tần số máy đo dù thực hiện trên thiết bị nam châm từ trường nhỏ hay các hệ phổ kế từ trường siêu dẫn cực mạnh.
Quy ước chất chuẩn đối chiếu theo IUPAC và thực tiễn kiểm nghiệm
Để đảm bảo tính nhất quán trên toàn thế giới, IUPAC đã tiêu chuẩn hóa việc lựa chọn chất chuẩn nội cho từng loại hạt nhân từ tính khác nhau trong khuyến nghị năm 2001.
Chất chuẩn trong dung môi hữu cơ
Đối với mẫu hòa tan trong các dung môi hữu cơ không phân cực hoặc phân cực trung bình, hợp chất tetramethylsilane, viết tắt là TMS, được lựa chọn làm chất chuẩn nội sơ cấp cho cả phổ proton và phổ carbon-13. TMS sở hữu cấu trúc đối xứng tứ diện cao với bốn nhóm metyl tương đương gắn vào nguyên tử silic trung tâm. Do silic có độ âm điện thấp hơn carbon, mật độ electron xung quanh các nguyên tử hydro và carbon trong phân tử TMS rất lớn, tạo nên hiệu ứng chắn nghịch từ cực mạnh. Do đó, tín hiệu của TMS xuất hiện ở vùng trường rất mạnh, tách biệt hoàn toàn khỏi hầu hết các tín hiệu của các hợp chất hữu cơ thông thường. IUPAC quy ước độ dịch chuyển hóa học của TMS bằng đúng 0,0 ppm.
Chất chuẩn trong dung dịch nước và mẫu sinh học
TMS không tan trong nước, vì vậy khi phân tích các mẫu sinh học, dược phẩm hoặc protein trong môi trường nước và nước nặng, IUPAC khuyến nghị sử dụng hai hợp chất dẫn xuất silan tan trong nước làm chất chuẩn thứ cấp:
- DSS: Muối natri của acid 2,2-dimethyl-2-silapentane-5-sulfonic, sở hữu nhóm trimetylsilyl tạo mũi cộng hưởng đơn sắc nét được quy ước tại 0,0 ppm.
- TSP: Muối natri của acid 3-(trimethylsilyl)propionic-d4, cũng có tín hiệu chuẩn quy ước tại 0,0 ppm và thường được dùng trong các nghiên cứu chuyển hóa sinh học.
Thực tiễn nghiên cứu và kiểm nghiệm tại Việt Nam
Tại các phòng thí nghiệm và cơ sở nghiên cứu phân tích cấu trúc hóa học tại Việt Nam, phương pháp phổ cộng hưởng từ hạt nhân là công cụ thiết yếu để giải mã cấu trúc phân tử các hợp chất tự nhiên phân lập từ dược liệu và tổng hợp hữu cơ. Trong các phép đo này, tetramethylsilane luôn được sử dụng làm chất chuẩn quy ước 0,0 ppm theo thông lệ quốc tế.
Các dải dịch chuyển hóa học đặc trưng trong thực nghiệm
Giá trị độ dịch chuyển hóa học phụ thuộc chặt chẽ vào độ âm điện của các nguyên tử lân cận, trạng thái lai hóa của orbital và hiệu ứng dị hướng từ trường của các liên kết bất bão hòa xung quanh.
Thang đo proton (1H-NMR)
Theo tài liệu tra cứu thực nghiệm chuẩn của Pretsch và các cộng sự (2009), thang dịch chuyển của proton thông thường nằm trong khoảng từ 0 ppm đến 14 ppm. Các vùng tín hiệu đặc trưng bao gồm:
| Nhóm chức hóa học | Khoảng dịch chuyển hóa học (ppm) | Đặc điểm cấu trúc và hiệu ứng chắn |
|---|---|---|
| Nhóm alkyl no (R-CH3, R2-CH2, R3-CH) | 0,5 đến 2,0 ppm | Mật độ electron cao, proton được che chở mạnh bởi các liên kết sigma. |
| Proton cạnh nhóm thế hút electron (CH-X) | Vùng trung gian (liền kề alkyl) | Hiệu ứng cảm ứng hút electron làm giảm mật độ điện tích quanh hạt nhân. |
| Proton liên kết alken (C=C-H) | Vùng trường yếu trung gian | Chịu ảnh hưởng của lai hóa sp2 và từ trường cảm ứng cục bộ của liên kết pi. |
| Proton vòng thơm (Ar-H) | 6,5 đến 8,5 ppm | Hiệu ứng dị hướng từ mạnh sinh ra từ dòng vòng electron pi của nhân benzen. |
| Proton nhóm aldehyde (R-CHO) | Vùng trường yếu liền kề nhân thơm | Chịu tác động đồng thời của liên kết đôi carbonyl phân cực và dị hướng từ. |
| Proton nhóm carboxylic acid (R-COOH) | 10,5 đến 13,0 ppm | Liên kết hydro liên phân tử rất mạnh làm suy giảm tối đa lớp vỏ chắn electron. |
Thang đo carbon-13 (13C-NMR)
Theo bảng dữ liệu thực nghiệm của Pretsch và các cộng sự (2009), do hằng số chắn thuận từ đóng vai trò chi phối và khoảng biến thiên năng lượng kích thích electron rộng hơn nhiều so với proton, thang đo phổ carbon-13 thực nghiệm trải rộng trong khoảng từ 0 ppm đến 220 ppm:
- Vùng carbon no aliphatic lai hóa sp3: Thường xuất hiện trong khoảng từ 0 ppm đến 50 ppm đối với các hydrocarbon mạch hở và mạch vòng no.
- Vùng carbon gắn dị tố (C-O, C-N, C-halogen): Dịch chuyển về phía trường thấp hơn, thường dao động trong khoảng từ 50 ppm đến 90 ppm.
- Vùng carbon lai hóa sp2 và nhân thơm: Nằm trong dải từ 100 ppm đến 160 ppm do tác động kết hợp của mật độ electron pi và dị hướng từ trường.
- Vùng carbon carbonyl: Xuất hiện ở vùng trường thấp nhất, trải từ 160 ppm đến 220 ppm, trong đó nhóm ester và amide thường nằm ở 160 đến 180 ppm, còn nhóm ketone và aldehyde có thể đạt tới 200 đến 220 ppm.
Các yếu tố chi phối độ dịch chuyển hóa học trong thực nghiệm
Trong thực hành phổ nghiệm, vị trí pic phổ không phải là một giá trị bất biến tuyệt đối mà chịu tác động rõ rệt bởi các điều kiện vật lý và hóa lý của môi trường đo:
Hiệu ứng dung môi và liên kết hydro
Dung môi sử dụng trong phép đo NMR không chỉ đóng vai trò hòa tan mẫu mà còn tương tác trực tiếp với các phân tử chất tan thông qua lực phân tán van der Waals, tương tác lưỡng cực và liên kết hydro. Sự tạo thành liên kết hydro làm kéo giãn mật độ electron xung quanh proton tham gia liên kết, dẫn đến sự dịch chuyển rõ rệt về phía trường thấp đối với các nhóm hydroxyl và amine. Ngoài ra, tính dị hướng từ của các phân tử dung môi có nhân thơm như pyridin-d5 hoặc benzen-d6 có thể làm thay đổi đáng kể độ dịch chuyển của các nhóm chức lân cận so với khi đo trong clorofom deuterat.
Năm 2010, Fulmer và các cộng sự đã công bố một nghiên cứu chuẩn hóa trên tạp chí Organometallics, tổng hợp chi tiết vị trí pic phổ của các dung môi hữu cơ dư và tạp chất thường gặp trong năm loại dung môi deuterat thông dụng. Công trình của nhóm Fulmer chỉ ra rằng vị trí hấp thụ của cùng một phân tử nước tạp chất có thể thay đổi rất lớn: trong clorofom deuterat tín hiệu nước xuất hiện tại 1,56 ppm, nhưng trong toluen deuterat tín hiệu này bị dịch chuyển mạnh về phía tần số thấp hơn và xuất hiện tại 0,43 ppm do hiệu ứng chắn từ trường của vòng thơm dung môi.
Nhiệt độ và nồng độ mẫu
Nhiệt độ đo làm thay đổi cân bằng dung môi hóa và phá vỡ liên kết hydro liên phân tử. Khi tăng nhiệt độ, các liên kết hydro yếu dần, khiến độ dịch chuyển của các proton linh động như OH và NH dịch dần về phía trường mạnh. Đối với các mẫu đo ở nồng độ cao, hiện tượng tự liên hợp phân tử có thể làm dịch chuyển pic và làm giãn rộng vạch phổ.
Phương pháp tính toán và dự đoán độ dịch chuyển hóa học
Cùng với sự phát triển của năng lực tính toán hiệu năng cao, việc kết hợp đo đạc thực nghiệm với mô phỏng lý thuyết đã trở thành quy trình tiêu chuẩn trong hóa học hiện đại để xác minh cấu trúc của các hợp chất tự nhiên phức tạp hoặc các sản phẩm tổng hợp hữu cơ mới.
Phương pháp hóa học lượng tử ab initio và DFT
Trong bài tổng quan hệ thống công bố năm 2011 trên tạp chí Chemical Reviews, Lodewyk, Siebert và Tantillo đã tổng kết phương pháp tính toán hóa học lượng tử dựa trên lý thuyết phiếm định mật độ để tính toán hằng số chắn từ và độ dịch chuyển hóa học 1H cũng như 13C. Bằng cách kết hợp phiếm hàm mật độ phù hợp cùng bộ cơ sở thích hợp và sử dụng các hệ số hiệu chỉnh tuyến tính tối ưu, các tác giả chứng minh rằng phương pháp DFT có thể đạt độ chính xác cao với sai số tuyệt đối trung bình dưới 0,15 ppm đối với proton và dưới 2,0 ppm đối với hạt nhân carbon-13. Độ chính xác này cho phép phân biệt rõ ràng giữa các đồng phân lập thể và cấu hình không gian đối quang khó phân lập bằng phương pháp tinh thể học.
Phương pháp học máy và mô hình dữ liệu
Mặc dù phương pháp DFT mang lại độ chính xác cao, thời gian tính toán cho các đại phân tử sinh học hoặc các thư viện hợp chất lớn thường kéo dài nhiều ngày. Nhằm giải quyết điểm nghẽn này, các phương pháp dự đoán dựa trên dữ liệu đã thu hút sự quan tâm mạnh mẽ của giới nghiên cứu. Trong bài tổng quan công bố năm 2021 trên tạp chí Magnetic Resonance in Chemistry, Jonas, Kuhn và Schlörer đã phân tích toàn diện các thuật toán học máy, mạng nơ-ron đồ thị và mô hình học sâu được huấn luyện trên các cơ sở dữ liệu phổ thực nghiệm công cộng. Các mô hình học máy hiện đại có khả năng dự đoán độ dịch chuyển hóa học của toàn bộ phân tử chỉ trong vài mili-giây với độ chính xác tiệm cận tính toán lượng tử lượng tử, mở ra triển vọng tự động hóa phân tích cấu trúc trong sàng lọc thuốc và hóa học tính toán.
Hạn chế và thách thức thực nghiệm
Mặc dù là công cụ đắc lực, việc giải đoán độ dịch chuyển hóa học trong thực tế vẫn đối mặt với một số ranh giới kỹ thuật quan trọng:
- Hiện tượng chồng lấn tín hiệu: Trong các phân tử có kích thước lớn như chuỗi polypeptide hay polysaccharide, nhiều hạt nhân nằm trong môi trường hóa học gần như tương đồng, dẫn đến sự chồng lấn nghiêm trọng của các pic phổ trong không gian một chiều, đòi hỏi phải triển khai các kỹ thuật phổ hai chiều và ba chiều phức tạp.
- Độ nhạy nội tại thấp: Hạt nhân carbon-13 chỉ chiếm tỷ lệ tự nhiên khoảng một phần trăm trong tổng lượng carbon, khiến cường độ tín hiệu rất yếu so với proton và đòi hỏi thời gian tích lũy tín hiệu kéo dài hoặc cần máy đo phổ trang bị đầu dò đông lạnh từ trường cao.
- Ảnh hưởng của trao đổi hóa học: Các proton linh động tham gia vào quá trình trao đổi nhanh với dung môi có thể tạo ra các pic phổ giãn rộng, làm mờ ranh giới pic hoặc biến mất hoàn toàn khỏi thang đo phổ nếu tốc độ trao đổi nằm trong dải thời gian trung gian của phép đo.